Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 6.6% intergenic (+28/+75) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3150CG 6.6% 413.7 / 14.1 141intergenic (+28/+75)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (62/68);  new base G (3/7);  total (65/76)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 143 total aligned reads displayed.

TACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGA  >  NC_000913/2286166‑2286464
                                                                                                                                                     |                                                                                                                                                      
tACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                        <  1:1796272/150‑1 (MQ=255)
    tGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCTCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                                    >  2:1485110/1‑150 (MQ=255)
     gCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                                                                                   >  1:554285/1‑150 (MQ=255)
      cAATCGGTGTGGAAATCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                                                                                  <  2:1867380/150‑1 (MQ=255)
      cAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                                                                                  >  1:854387/1‑150 (MQ=255)
        aTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggtt                                                                                                                                                <  1:1055598/150‑1 (MQ=255)
          gggTGGGGAAAAGGGTAGAATTAGCACCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAGGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGACGACCAAAAATCCCAAGAAAAAACCAACAGTAAGTGTTGGtttt                                                                                                                                             <  2:364397/149‑1 (MQ=255)
           ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCTAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                            >  2:1784666/1‑150 (MQ=255)
           ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                            >  2:579332/1‑150 (MQ=255)
           ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                            >  2:1017339/1‑150 (MQ=255)
            gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCATGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                           >  2:655525/1‑150 (MQ=255)
            gtgtGGAAAACGGTAGTATTAGCAGCCACGAGGCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                           <  1:490497/150‑1 (MQ=255)
             tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACTCTTACGGTTGGttttttt                                                                                                                                          <  2:882374/150‑1 (MQ=255)
              gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                                                                         <  2:1077228/150‑1 (MQ=255)
               tgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttctt                                                                                                                                        <  1:587758/150‑4 (MQ=255)
                ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                       >  1:896093/1‑150 (MQ=255)
                ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACAAACCATTACGGTTGGTTTTTTtata                                                                                                                                       >  2:1873267/1‑150 (MQ=255)
                     aCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                                                                                  <  1:1073678/150‑1 (MQ=255)
                      cGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                                                                                 <  1:1446408/150‑1 (MQ=255)
                          agtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACAAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                                                                             <  2:1941405/150‑1 (MQ=255)
                               tAGCAGCCACGAGTCGGCACGTAGCGGAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                        <  2:910384/150‑1 (MQ=255)
                               tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                        >  2:217812/1‑150 (MQ=255)
                               tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                        <  2:1470900/150‑1 (MQ=255)
                               tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                        <  1:1702891/150‑1 (MQ=255)
                                 gcagcCCCGCGTCGGCACGTATCGCAGCCTGGTAGCGCACCGTCATGGGGGGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGGTGGTTTTTTTATATCTGCAATTAATTCGAt                                                                                                                      <  1:1270712/150‑1 (MQ=255)
                                 gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCGTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCGTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                                                                                      >  1:1966635/1‑150 (MQ=255)
                                 gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAGGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAATAACCAACCCTTACGGTTGGTTGTTTTATATCTGCAATTAATTCGAt                                                                                                                      >  2:1721560/1‑150 (MQ=255)
                                 gaagcCACGAGTAGGCAGGTAGCGCAGCCTGGTAGCGCACCTTCATGGGTTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                                                                                      <  1:1532708/148‑1 (MQ=255)
                                  cagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                                                                                     <  2:1980809/150‑1 (MQ=255)
                                  cagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                                                                                     <  2:870140/150‑1 (MQ=255)
                                  cagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGGCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCTACCCTTACGGTGGGTTTTTTTATATCTGCAATTAATTCGATa                                                                                                                     >  2:1350409/1‑150 (MQ=255)
                                   agcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaa                                                                                                                    <  1:638496/150‑1 (MQ=255)
                                    gcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                   <  1:1540675/150‑1 (MQ=255)
                                    gcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                   >  2:400651/1‑150 (MQ=255)
                                     ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                                  >  2:1281412/1‑150 (MQ=255)
                                     ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                                  <  1:1399848/150‑1 (MQ=255)
                                      cACGAGTCGGCACGTAGCGCATCCTGGTATCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                                                                                 <  1:1951246/150‑1 (MQ=255)
                                      cACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                                                                                 <  2:556546/150‑1 (MQ=255)
                                       aCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAg                                                                                                                >  2:191488/1‑150 (MQ=255)
                                            tCGGCAAGTAGCGCAGCCTGGTAGAGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCTACCAAAAATCACAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                                                                           <  1:1319751/150‑1 (MQ=255)
                                              ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTGACGGTTGGTTTTTTTATATCAGCAATTAATTCGATAAACAGACCgaga                                                                                                         >  2:357450/1‑150 (MQ=255)
                                                  cGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                                                                                     >  1:1081108/1‑150 (MQ=255)
                                                    tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTATTTTGATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                                                                    >  1:1389290/1‑150 (MQ=255)
                                                         cAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                                                                              <  1:1853823/150‑1 (MQ=255)
                                                             cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGGCACATCACAGCCTGt                                                                                          >  2:1974359/1‑150 (MQ=255)
                                                              tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                                                                         >  2:527984/1‑150 (MQ=255)
                                                               ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttt                                                                                        <  2:403688/150‑1 (MQ=255)
                                                                   gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                                                    >  2:1413745/1‑150 (MQ=255)
                                                                       aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                                                                                <  2:1473125/150‑1 (MQ=255)
                                                                        ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                               <  2:497787/150‑1 (MQ=255)
                                                                            cATGGGGTGTCGGTGGTCTGAGGTTCAAATCCTCTCGTGACGACCAAAAATCCCAAGAAAAATCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCAATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                           <  1:1041963/150‑1 (MQ=255)
                                                                            cATGGGGGGTCGGGGGTCTGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGTTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                           <  2:1223460/150‑1 (MQ=255)
                                                                             aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                                                                          <  1:413561/150‑1 (MQ=255)
                                                                             aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                                                                          >  2:636924/1‑150 (MQ=255)
                                                                                gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                                                                       >  1:741137/1‑150 (MQ=255)
                                                                                   tgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGt                                                                    >  1:1596805/1‑150 (MQ=255)
                                                                                      cGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa                                                                 <  2:994068/150‑1 (MQ=255)
                                                                                           gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcaa                                                            >  1:727989/1‑148 (MQ=255)
                                                                                           gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca                                                            <  2:98146/150‑1 (MQ=255)
                                                                                            tCGGAGGTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGCCACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac                                                           <  1:125663/150‑1 (MQ=255)
                                                                                              ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCacac                                                         >  2:2015952/1‑150 (MQ=255)
                                                                                                 ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                      >  1:421419/1‑150 (MQ=255)
                                                                                                   ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                                    >  1:1669974/1‑150 (MQ=255)
                                                                                                   ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCAAGACCACACCCGcc                                                    >  1:621084/1‑150 (MQ=255)
                                                                                                    tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCACCCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGCCCGTGACACATCACAGCGTGTTTATTTTCTGTTATCAGAACGTCCAGCCCCCACCCGccg                                                   >  1:1091751/1‑149 (MQ=255)
                                                                                                    tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt                                                   >  1:1809130/1‑150 (MQ=255)
                                                                                                      aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCCCACCCGCCTGa                                                 >  1:333674/1‑150 (MQ=255)
                                                                                                       aaTCCTCTCGTGCCTACCAAAAATCCGAAGAAAAAACCAAACATTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGTCCACACCCGCCTGAg                                                <  2:1328964/150‑1 (MQ=255)
                                                                                                       aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                                <  2:896093/150‑1 (MQ=255)
                                                                                                       aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                                <  1:1440374/150‑1 (MQ=255)
                                                                                                       aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                                <  2:1596805/150‑1 (MQ=255)
                                                                                                       aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTAGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                                >  2:871711/1‑150 (MQ=255)
                                                                                                           ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTGCGGTTGGTTTTTTTAGATCTGCATTTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGATCGTCCAGACCACACCCGCCTGAGTAtt                                            >  2:1435252/1‑150 (MQ=255)
                                                                                                           ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCATTACGGTTGGTTTTTTTATATCTGCAATTAAATCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTAtt                                            >  1:1607385/1‑150 (MQ=255)
                                                                                                                 tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACgg                                      <  1:1041544/150‑1 (MQ=255)
                                                                                                                       ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGCCACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACGCCACCCGCCTGAGTATTCCACGGCGATTc                                >  1:68121/1‑150 (MQ=255)
                                                                                                                        cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                               >  2:236268/1‑150 (MQ=255)
                                                                                                                        cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                               >  2:1563883/1‑150 (MQ=255)
                                                                                                                         aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCACCGGCGATTCTa                              >  1:1020993/1‑150 (MQ=255)
                                                                                                                          aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAc                             <  2:1063518/150‑1 (MQ=255)
                                                                                                                          aaaaTACCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGACTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAc                             <  2:237086/150‑1 (MQ=255)
                                                                                                                           aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAcc                            <  2:1139811/150‑1 (MQ=255)
                                                                                                                                ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                       >  1:1514290/1‑150 (MQ=255)
                                                                                                                                ccAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTTTTTTCTGTTATCAGAACGTCAAGACCACACCCGCCTGAGTATTCGACGGCGATTCTACCCCTgc                       >  2:1247746/1‑150 (MQ=255)
                                                                                                                                 cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                      <  2:931065/150‑1 (MQ=255)
                                                                                                                                  aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcgc                     <  1:579332/150‑1 (MQ=255)
                                                                                                                                    gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                   <  2:1929054/150‑1 (MQ=255)
                                                                                                                                    aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                   >  2:1677112/2‑150 (MQ=255)
                                                                                                                                     aaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGAATTCCACGGCGATTCTAGACCTGCGCCAt                  >  1:439275/1‑150 (MQ=255)
                                                                                                                                      aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt                 >  1:880721/1‑150 (MQ=255)
                                                                                                                                       aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTg                <  2:672535/150‑1 (MQ=255)
                                                                                                                                        aaaccaaaccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGATAGTCCATACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGc               <  2:1691461/150‑1 (MQ=255)
                                                                                                                                        aaaCCAACCCTTTAGGTTGGTTTTTTTATATCTGCAATTAATTCGATAACCAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACCCCCGCCTGTGTATTCCACGGCGATTCTACCCCTGCGCCATTGc               <  1:776901/150‑1 (MQ=255)
                                                                                                                                        aaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGc               <  2:554445/150‑1 (MQ=255)
                                                                                                                                         aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc              <  2:610450/150‑1 (MQ=255)
                                                                                                                                           ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt            >  2:932004/1‑150 (MQ=255)
                                                                                                                                             aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTACTGTTATCAGAACGTCCAGACCACACCCGCCGAAGAATTCCACGGCAATTCAACCCCTCCCCAATTGCCaaaa          >  1:779424/1‑147 (MQ=255)
                                                                                                                                                ccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGGCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctg       <  2:334314/150‑1 (MQ=255)
                                                                                                                                                    aCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctgactg   <  2:1729898/150‑1 (MQ=255)
                                                                                                                                                     cGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctgactga  >  1:276499/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACTGA  >  NC_000913/2286166‑2286464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: