Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,794 (C)7→6 intergenic (‑18/‑124) yfjM ← / → rnlA CP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,7880C.99.3% 596.7 / ‑2.6 142intergenic (‑12/‑130)yfjM/rnlACP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (48/93);  total (49/93)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 163 total aligned reads displayed.

GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT  >  NC_000913/2765650‑2765919
                                                                                                                                             |                                                                                                                                    
ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                                     <  1:1675132/140‑1 (MQ=255)
ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                                     >  2:205589/1‑140 (MQ=255)
 gCTCTCTGCTCGTTCTGTTTCGACGAAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccc                                                                                                                                    >  2:1369594/1‑140 (MQ=255)
  ctctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                                   <  1:1555962/140‑1 (MQ=255)
   tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                  <  1:822942/140‑1 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGTACATCATCCCCCATTTTGTGGACTGccacccg                                                                                                                                >  2:1180571/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                >  1:96729/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                >  2:1640018/1‑139 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                                <  2:144398/140‑2 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                                   <  1:746225/137‑1 (MQ=255)
     tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                                   >  2:746225/1‑137 (MQ=255)
          tCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGa                                                                                                                          <  2:1221197/140‑1 (MQ=255)
          tCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGa                                                                                                                          <  2:139796/140‑1 (MQ=255)
          tCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGa                                                                                                                          <  1:337885/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:589289/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:1397069/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:1425908/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:615614/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  2:441502/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:1586175/140‑1 (MQ=255)
               cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑CGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATc                                                                                                                     <  1:1311270/140‑1 (MQ=255)
                 gTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCt                                                                                                                   >  1:1724353/1‑140 (MQ=255)
                 gTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCt                                                                                                                   >  2:423467/1‑140 (MQ=255)
                  tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTg                                                                                                                  <  1:181341/140‑1 (MQ=255)
                   ttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                  <  1:1171218/124‑1 (MQ=255)
                   ttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                                  >  2:1171218/1‑124 (MQ=255)
                    tCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGcc                                                                                                                <  2:1612565/140‑1 (MQ=255)
                     cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGccc                                                                                                               <  2:498995/140‑1 (MQ=255)
                     cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGccc                                                                                                               <  2:114043/140‑1 (MQ=255)
                     cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGccc                                                                                                               <  2:1176419/140‑1 (MQ=255)
                       aCGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGcc                                                                                                                <  2:1327311/140‑1 (MQ=255)
                                 ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                                   <  2:1773777/140‑1 (MQ=255)
                                 ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                                   <  2:1417295/140‑1 (MQ=255)
                                 ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                                   <  1:423467/140‑1 (MQ=255)
                                    tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                               <  1:1177276/140‑1 (MQ=255)
                                    tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                               <  1:265222/140‑1 (MQ=255)
                                     cTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                              >  1:748329/1‑140 (MQ=255)
                                           tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                    >  2:135708/1‑129 (MQ=255)
                                           tgctgGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAg                                                                                                    <  1:135708/129‑1 (MQ=255)
                                                  gAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTgg                                                                                 <  1:658618/140‑1 (MQ=255)
                                                   aataatTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGc                                                                                >  2:904872/1‑140 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  1:662319/140‑1 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:1724353/140‑1 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:748329/140‑1 (MQ=255)
                                                            tAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  1:969670/140‑1 (MQ=255)
                                                              aGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAGATTTAGCa                                                                     >  2:132108/1‑140 (MQ=255)
                                                               gTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCaa                                                                    <  1:1624543/140‑1 (MQ=255)
                                                                tATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCaaa                                                                   >  1:1659733/1‑140 (MQ=255)
                                                                  tCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                  >  2:981158/1‑140 (MQ=255)
                                                                  tCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCa                                                                     <  1:1310272/136‑1 (MQ=255)
                                                                   cGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCTCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                                                                >  2:596731/1‑140 (MQ=255)
                                                                   cGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAATGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                                                                >  1:41457/1‑140 (MQ=255)
                                                                     aGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                              >  2:111517/1‑140 (MQ=255)
                                                                      gTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                             <  1:1608072/140‑1 (MQ=255)
                                                                         gctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                          >  2:760826/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       >  1:1394106/1‑126 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       <  2:1394106/126‑1 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                                       >  1:785298/1‑126 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  1:1667452/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  2:1532561/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  1:521035/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  2:1741764/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  2:688308/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  2:1035614/1‑140 (MQ=255)
                                                                          ctgct‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTACCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                         >  2:822654/1‑140 (MQ=255)
                                                                              t‑‑‑tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                     <  1:205589/140‑1 (MQ=255)
                                                                              t‑‑‑tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                     <  1:46623/140‑1 (MQ=255)
                                                                              t‑‑‑tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                     <  2:611238/140‑1 (MQ=255)
                                                                              t‑‑‑tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                     <  2:301995/140‑1 (MQ=255)
                                                                              t‑‑‑tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                     <  2:1667452/140‑1 (MQ=255)
                                                                                   ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCTCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                   <  1:596731/140‑1 (MQ=255)
                                                                                      gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                <  2:473686/140‑1 (MQ=255)
                                                                                        gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                              <  2:1258373/140‑1 (MQ=255)
                                                                                        gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                              <  1:1369594/140‑1 (MQ=255)
                                                                                        gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                              <  2:941258/140‑1 (MQ=255)
                                                                                         aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat                                             >  1:1487717/1‑140 (MQ=255)
                                                                                               aaGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGt                                       <  1:599119/140‑1 (MQ=255)
                                                                                                  agGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                    <  1:228227/140‑1 (MQ=255)
                                                                                                    gTTTCCGAAACGATGGACATCATCCCCCATCTTGTGGACTGCCTCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                  <  2:781516/140‑1 (MQ=255)
                                                                                                    gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                  <  2:1487717/140‑1 (MQ=255)
                                                                                                    gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                  <  2:642526/140‑1 (MQ=255)
                                                                                                    gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                  >  2:704096/1‑140 (MQ=255)
                                                                                                    gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGtt                                 >  2:1810508/1‑140 (MQ=255)
                                                                                                      tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGttt                                >  2:1014371/1‑140 (MQ=255)
                                                                                                      tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGttt                                >  2:1367542/1‑140 (MQ=255)
                                                                                                      tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGttt                                >  1:1529043/1‑140 (MQ=255)
                                                                                                        ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                              >  2:1457518/1‑140 (MQ=255)
                                                                                                        ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                              <  1:1457518/140‑1 (MQ=255)
                                                                                                          gAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcccc                            <  1:1014371/140‑1 (MQ=255)
                                                                                                             aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata                         <  1:200199/140‑1 (MQ=255)
                                                                                                              cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCatat                        >  2:150003/1‑140 (MQ=255)
                                                                                                               gATGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATg                       <  1:1761305/140‑1 (MQ=255)
                                                                                                                 tGGACATCATCCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGtt                     >  1:1029671/1‑140 (MQ=255)
                                                                                                                       tcatcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg               <  1:730594/140‑1 (MQ=255)
                                                                                                                       tcatcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg               <  1:1367542/140‑1 (MQ=255)
                                                                                                                        catcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                          >  1:1243497/1‑128 (MQ=255)
                                                                                                                          tcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg               <  1:784975/137‑1 (MQ=255)
                                                                                                                          tcCCCCATTTTGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg               >  2:784975/1‑137 (MQ=255)
                                                                                                                                 ttttGTGGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGa     <  1:1463469/140‑1 (MQ=255)
                                                                                                                                    tgtgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa   >  2:692627/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                    
GGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCT‑‑‑TGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAT  >  NC_000913/2765650‑2765919

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: