Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,832 C→G 7.1% intergenic (+31/+108) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8320CG 7.1% 287.8 / 10.2 84intergenic (+31/+108)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (42/36);  new base G (2/4);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGT  >  NC_000913/2282690‑2282974
                                                                                                                                              |                                                                                                                                               
cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTTGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGGTCCG‑GCggg                                                                                                                                          <  2:757154/149‑1 (MQ=255)
cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                           <  3:607684/148‑1 (MQ=255)
cGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCggg                                                                                                                                          <  2:675488/149‑1 (MQ=255)
 gAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCgggg                                                                                                                                         >  1:624524/1‑149 (MQ=255)
     ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtttt                                                                                                                                     <  1:1068598/149‑1 (MQ=255)
     ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttt                                                                                                                                     <  2:595544/149‑1 (MQ=255)
      tCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGtctt                                                                                                                                     >  4:139141/1‑148 (MQ=255)
      tCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttt                                                                                                                                     <  3:263980/148‑1 (MQ=255)
       cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttttt                                                                                                                                    <  3:492415/148‑1 (MQ=255)
       cTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGttt                                                                                                                                      <  1:873904/146‑1 (MQ=255)
             aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCt                                                                                                                             <  2:26457/149‑1 (MQ=255)
             aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCt                                                                                                                             <  1:532780/149‑1 (MQ=255)
             aTCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCt                                                                                                                             >  1:614124/1‑149 (MQ=255)
                     tgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTa                                                                                                                     <  1:1003827/149‑1 (MQ=255)
                       aCGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTAt                                                                                                                    <  3:468727/148‑1 (MQ=255)
                         ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGcc                                                                                                                 <  1:904198/149‑1 (MQ=255)
                         ggTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGc                                                                                                                  <  4:390985/148‑1 (MQ=255)
                          gTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGcc                                                                                                                 >  3:799269/1‑148 (MQ=255)
                              aGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGAtt                                                                                                            <  1:104050/149‑1 (MQ=255)
                               gAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGAtt                                                                                                            >  4:351326/1‑148 (MQ=255)
                                   tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGt                                                                                                       >  1:492637/1‑149 (MQ=255)
                                     aaaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                     <  2:27967/149‑1 (MQ=255)
                                      aaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                     >  3:500392/1‑148 (MQ=255)
                                        gTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttt                                                                                                                                      <  3:822986/113‑1 (MQ=255)
                                        gTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttt                                                                                                                                      >  4:822986/1‑113 (MQ=255)
                                                     aCGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                      >  3:488914/1‑148 (MQ=255)
                                                       gTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGc                                                                                    <  4:317277/148‑1 (MQ=255)
                                                         aCAGCGTCACTCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttttta                                                                                                                                  <  2:647033/100‑2 (MQ=38)
                                                         aCAGCGTCACTCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGC‑CGGAGCGGGGtttttta                                                                                                                                  >  1:647033/1‑99 (MQ=38)
                                                          cAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAg                                                                                 >  4:421681/1‑148 (MQ=255)
                                                            gCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATa                                                                              >  2:344739/1‑149 (MQ=255)
                                                             cGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                             <  2:492637/149‑1 (MQ=255)
                                                             cGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATa                                                                              >  4:385133/1‑148 (MQ=255)
                                                                cACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                           <  3:385133/148‑1 (MQ=255)
                                                                 aCCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTCTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGc                                                                          <  3:139141/148‑1 (MQ=255)
                                                                    cGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                       >  3:426395/1‑148 (MQ=255)
                                                                      tACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAc                                                                    >  2:235081/1‑149 (MQ=255)
                                                                       aCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCAGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAc                                                                    >  4:695989/1‑148 (MQ=255)
                                                                           aCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgccgc                                                               <  2:638879/149‑1 (MQ=255)
                                                                           aCCGAAGCCGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgccg                                                                <  3:705792/148‑1 (MQ=255)
                                                                            ccGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCa                                                              >  1:1054994/1‑149 (MQ=255)
                                                                             cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGt                                                                        <  2:635825/138‑1 (MQ=255)
                                                                             cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGt                                                                        >  1:635825/1‑138 (MQ=255)
                                                                              gAAGCTGCAGCACATCGACTTTGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATc                                                            >  1:774102/1‑149 (MQ=255)
                                                                              gAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATc                                                            >  1:501776/1‑149 (MQ=255)
                                                                                aGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAAcgcctt                                                               >  4:489030/1‑142 (MQ=255)
                                                                                  cTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTgg                                                         <  4:939171/148‑1 (MQ=255)
                                                                                     aggcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAt                                                      >  4:657065/3‑148 (MQ=255)
                                                                                       gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                   <  2:429118/149‑1 (MQ=255)
                                                                                        cacaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                   <  3:351326/148‑1 (MQ=255)
                                                                                          caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGGTAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                 >  3:375473/1‑148 (MQ=255)
                                                                                          caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                 <  3:222116/148‑1 (MQ=255)
                                                                                           aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                               >  2:976651/1‑149 (MQ=255)
                                                                                           aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                               >  1:487612/1‑149 (MQ=255)
                                                                                                cTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaa                                           >  4:376569/1‑148 (MQ=255)
                                                                                                 ttcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaaa                                          <  4:799269/148‑1 (MQ=255)
                                                                                                  tcgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAc                                        >  2:374903/1‑149 (MQ=255)
                                                                                                    gttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCg                                      >  1:807781/1‑149 (MQ=255)
                                                                                                     ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCg                                      <  4:488914/148‑1 (MQ=255)
                                                                                                     ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCg                                      <  4:311886/148‑1 (MQ=255)
                                                                                                      tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                           >  3:182412/1‑110 (MQ=255)
                                                                                                      tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                           <  4:182412/110‑1 (MQ=255)
                                                                                                           cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAtt                               <  2:624524/149‑1 (MQ=255)
                                                                                                           cTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAAt                                >  3:541102/1‑148 (MQ=255)
                                                                                                              aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGCGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgcc                            >  1:880696/1‑149 (MQ=255)
                                                                                                                     gAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAcc                                       <  2:349552/131‑1 (MQ=255)
                                                                                                                     gAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAAcc                                       >  1:349552/1‑131 (MQ=255)
                                                                                                                      aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                    >  1:971857/1‑149 (MQ=255)
                                                                                                                        tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                  <  2:774102/149‑1 (MQ=255)
                                                                                                                            ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc              >  1:866195/1‑149 (MQ=255)
                                                                                                                             tacaaaaaaCCCCGCCGGA‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgcg             >  1:731713/4‑149 (MQ=39)
                                                                                                                               tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt           >  1:701382/1‑149 (MQ=255)
                                                                                                                                 aaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc         >  2:850341/1‑149 (MQ=255)
                                                                                                                                  aaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGAATACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAAGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa        >  2:1005407/1‑149 (MQ=255)
                                                                                                                                    aaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAt      >  2:328751/1‑149 (MQ=255)
                                                                                                                                    aaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAAt      <  1:344739/149‑1 (MQ=255)
                                                                                                                                     aCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTa                                  <  3:989321/120‑1 (MQ=255)
                                                                                                                                     aCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTa                                  >  4:989321/1‑120 (MQ=255)
                                                                                                                                       cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg   >  1:338100/1‑149 (MQ=255)
                                                                                                                                       cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg    >  3:57184/1‑148 (MQ=255)
                                                                                                                                       cccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg    <  3:376569/148‑1 (MQ=255)
                                                                                                                                        ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGGTAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg   <  4:375473/148‑1 (MQ=255)
                                                                                                                                        ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg   >  3:855364/1‑148 (MQ=255)
                                                                                                                                        ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg   >  3:466576/1‑148 (MQ=255)
                                                                                                                                         cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt  <  4:500392/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGT  >  NC_000913/2282690‑2282974

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑