Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,037,4960AT41.7% 45.8 / 88.3 72intergenic (+161/‑23)alaT/rrlAtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (24/18);  new base T (16/14);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGG  >  NC_000913/4037248‑4037641
                                                                                                                                                                                                                                                        |                                                                                                                                                 
aCTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATATCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgac                                                                                                                                                  >  1:101023/1‑249 (MQ=25)
    cATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgtaagaa                                                                                                                                              >  1:68242/1‑250 (MQ=25)
      ttttATGGGGCTATAACTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGCTTAAATACAACATACTTCAGAGTGTACCTGCACAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGCAACACTGAACAACGAAAGTTTTTCGTGAGTCTCTCAACTTTTCGCAACACGATGATGAATCGAAAGAAAcc                                                                                                                                           >  2:377819/3‑250 (MQ=11)
        ttATGGGGCTATAGCTCAGCTTGGAGAGCGCCTGCGTTGCACGCAGGAGTTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTTAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACAt                                                                                                                                          <  2:610614/250‑1 (MQ=255)
           tGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATATTc                                                                                                                                      >  1:294934/1‑251 (MQ=37)
           tGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTAACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAACCATCTTc                                                                                                                                      >  1:70944/1‑251 (MQ=21)
             ggtcTATAGCGCAGCTGGGAGAGCGCCTGCTTTGAACGCAGGAGGACTGCGGTTCGATCCCGCATAGCTCCTCCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTAACTGCTAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAATCTGAAAATTGAAACACTGAACAACGAAAGTTGTTGGTGAGTCTCTCAAATTTTCGCAACACTAGGATGAATCTAAATAAACATCTTCgt                                                                                                                                    <  2:216361/248‑2 (MQ=11)
                        cAGCTGGTAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCTGTTCGATCCCGCATAGCTCCCCCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAgg                                                                                                                           <  1:374382/249‑1 (MQ=255)
                         aGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTGCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTa                                                                                                                        <  1:357371/251‑1 (MQ=21)
                            tGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGGTGTGAGGTTAAGc                                                                                                                     >  2:250760/1‑251 (MQ=37)
                            tGGGAGAGCGCCTGCTTTGAACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACTTGCAAAGGTTCACTGCGATGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAGAATTGAAAAACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATGGAATGAAACAACTTCGGGTAGGGAGGTTAAGc                                                                                                                     >  2:359657/1‑251 (MQ=11)
                                agagCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                                                                                                                 >  1:632513/1‑251 (MQ=21)
                                          ccttGCGCGCAGGAGGTCTGCGGTTCGAGCCCGCATAGCTCCACCATCTCTGGAGTGATTAAATAAAAAATGCTTCAGAGTGTACCGGCAAAGGTTCACTGCGAAGGTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGATCAACGAAAGTTGTTCGTGAGTCTCTAAAATTTTCGCAACACGATGATGAATCGAAAGAAACAGCTTCGGGTTGGGAGGTTAAGCGACTAAGTGTacac                                                                                                       <  2:474826/249‑1 (MQ=11)
                                             tGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt                                                                                                    >  1:586257/1‑251 (MQ=21)
                                                  gCAGGAGGTCTGCGGTTGGATCCCGCATAGATCCACCATCTATGTAGAGATTAAATAAAAAATACTTCAGAGTGTAGCTGCAAAGGATCAATGAGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCAGAAAATTGTACCACTGCACAACGAAATATGTGCGTGAGTCTCTCAAGTTTTCGAAACTCGAGGATGAAGCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGGACACGGTGGAtt                                                                                               >  2:513623/1‑250 (MQ=11)
                                                         gTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGAACACGGTGGATGCCCTGGc                                                                                        >  1:297462/1‑251 (MQ=37)
                                                               ggttcgaaaacgaaTAGCGCCACCATCACTGTAGTGATTAAATAAAAAATACTTCAGAGTGAACCTGCAAAGGTTCACTGCGAAGTTTTGCTAATTAAACATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTATCTCAAATTTTCGCCACACGATGATGAATCGAAAGAAACATATACGGGATGTGGGGGTAAGCGACTAAGCGTACAAGGTGGATGCCCTGGCAGTc                                                                                    >  2:128710/14‑249 (MQ=11)
                                                               ggTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCag                                                                                  <  1:611101/251‑1 (MQ=21)
                                                               ggTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTc                                                                                    >  2:251782/1‑249 (MQ=37)
                                                                  tCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCagag                                                                                <  1:400139/250‑1 (MQ=25)
                                                                   cGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGc                                                                              >  1:113600/1‑251 (MQ=21)
                                                                      tCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTTTTCGTGAGTCTCTCAAATTTCCGCAACACGATGTGGAATCGAAAGAACAATCTTCGGGTTGTGAGGTTAAGCACCTAGCCGTACACGGTGGATGCCCTGGCGGTcggggggg                                                                             >  2:134691/1‑242 (MQ=11)
                                                                      tCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGAt                                                                           >  1:94804/1‑251 (MQ=37)
                                                                        ccGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGa                                                                         <  1:251782/251‑1 (MQ=37)
                                                                         cGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGCTGa                                                                         >  1:125726/1‑250 (MQ=37)
                                                                            aTAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGCGGCGATGAAgg                                                                      >  2:568936/1‑250 (MQ=37)
                                                                            aTAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGa                                                                     <  1:364328/251‑1 (MQ=37)
                                                                             tAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGCTGATGAATCGAAAGAAACATCTTCGGGTTGTGAGATTAAGCGACTAAGCGTACACGGAGGATGCACTGGCAGTCAGAGGCGATGAAAGAc                                                                    >  2:563619/1‑251 (MQ=255)
                                                                                   caccaTCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCt                                                               >  2:654359/1‑250 (MQ=37)
                                                                                      caTCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATc                                                           <  1:539377/251‑1 (MQ=21)
                                                                                          tctGTAGTGATAACATAAAAACTACTTCAGAGTGTCCCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGAGGATGCCCTGGCAGTCAGAGGCGCTGAAGGACGTGCTAATCTGCg                                                       >  2:292353/1‑251 (MQ=11)
                                                                                            tGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGAt                                                     >  1:422797/1‑251 (MQ=37)
                                                                                             gTAGTGATTAAATAAAAAATACTTCAGAGGGTACCTGCAAATGTTCACTGCGAAGGTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAAGTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCGCTCAAATTTTCGCAACACGATGATGAATCGAAATAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCGTGGCAGTCAGAGGCGATGAAGGACGTGCTAACCTGCGATa                                                    <  2:129948/251‑1 (MQ=11)
                                                                                               aGTGATTAAATAAAAAATACTTCAGAGTGTTCCCGCAAAGGTTCACTGCGAAGGTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTGTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAg                                                  <  1:41135/251‑1 (MQ=255)
                                                                                               aGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGAGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATa                                                    >  1:299257/1‑249 (MQ=25)
                                                                                               aGTGATTAAATAAAAAATACTGCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAg                                                  <  2:209448/251‑1 (MQ=21)
                                                                                                   aTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTc                                              >  2:213302/1‑251 (MQ=35)
                                                                                                   aTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTc                                              <  1:509664/251‑1 (MQ=35)
                                                                                                    ttAAATAGAAAATACTTCAGAGTGTACCTGCAAAGCTTCACTGTGAAGTTTTGATATTTACCAATCTGGATCAAGCTGAAAATTGAAACACTCAACACCGAAAGTTATTCGTGAGTCTCTCAAATTTTCGCGACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACGCGGTGTATGCCCTGGCAGTCAGAGGCAATGAAGGACGTGCTAATCTGCGATAAGCGTCg                                             >  2:353549/1‑251 (MQ=11)
                                                                                                      aaataaaAAATACTTCAGTGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGATATTTAAAAATCTGGATCAAGATGAAAATTGCAACACTGAAAAACGCAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTAGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGAGAATGCCCTGGCAGTCAGAGGCGATGAAGGAGGTGCTAATATGCGAGAAGCGTCggg                                           >  2:628451/1‑250 (MQ=11)
                                                                                                      aaataaaAAATACTTCAGAGTGTACCTGCAAAGGTTAACTGCGAAGTTTTGCTCTGTAAAAAGCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGt                                           <  2:198149/251‑1 (MQ=16)
                                                                                                       aataaaACATACTTCATAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGACCAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCGTCGGGTTGTGAGGTTAAGCGACTAAGCGTACAGGGTGGATGCCCTGGCTGTCAGAGGCGATGAAGGACGGGCTAATCTGCAATAAGCGCCGGt                                           >  2:39319/1‑250 (MQ=11)
                                                                                                          aaaaaaTACTGCAGAGTGTACCTGCAAAGGTTTACTGCGAAGTTTTGCTCTTTAAAAATCTGGATTAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCATCACGATGATGAATCGAAAGAACCATCTTCGTGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAgg                                       <  1:633777/251‑1 (MQ=255)
                                                                                                            aaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAGGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTACTCTGCGATAAGCGTCGGTAAGGTg                                     >  1:602829/1‑251 (MQ=18)
                                                                                                             aaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTg                                     >  1:574542/1‑250 (MQ=35)
                                                                                                              aaTACTTCAGAGTGTACCTGCAACGGGTCACTGCGAAGTTTTGCTCTTTCAAACTATGGATCAAGCTGAAAATTGAAACACTGAACAACGTACGTTGTTCGTGAGTCTCTCAGATTTTCGCAACCCGATGATGAATCGAAGGAACCATCTTCGGGTTGTGAGGTTAAGCGACTAGGCGAACACGGTGGGTGCCCTGCCAGTCCGGGGCGATGACAGACGTGCTAATCGGCGATAAGCGTCGGTAAGGGGat                                   >  2:67153/1‑251 (MQ=11)
                                                                                                               aTACTTCAGAGTTTACCGGCAAAGGTTCACTGCGAAGTTTAGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCTATGAAGGACGTGCGAATCTGCGATAAGCGTCGTTAAGGTGata                                  <  2:83703/251‑1 (MQ=255)
                                                                                                               aTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGCTCAAGCTGAAAATTGCAACACTGAACAACGACAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCATGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGat                                   >  2:106996/1‑250 (MQ=255)
                                                                                                                 aCTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATg                                <  2:629948/251‑1 (MQ=18)
                                                                                                                    tCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGaa                              <  1:588576/250‑1 (MQ=35)
                                                                                                                    tCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGaa                              <  1:160025/250‑1 (MQ=35)
                                                                                                                      agagTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGa                                                      >  2:119478/1‑224 (MQ=35)
                                                                                                                      agagTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGa                                                      <  1:119478/224‑1 (MQ=35)
                                                                                                                       gagTGTACCTGCAAAGGTTCACTACGAAGTTTTGCTCTTTCAAAATATGGATAAAGCTGAAAATTAAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTGAAATTTTCGCAACACGCTGATGAATCGTAAGAAACATCTTCGGGTGGTGAGGTTAAGCGCCTAAGCGTACACGGGGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAAGATGCGATCAGCGTCGGTAAGGGGATATGAACCGt                          >  2:483016/1‑251 (MQ=11)
                                                                                                                        agggtACCTGCAAAGGTTCACTGCGAAGTTTTGCGCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCgga                         <  2:297462/248‑3 (MQ=21)
                                                                                                                        agTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGtt                         <  2:94804/251‑1 (MQ=34)
                                                                                                                         gtgtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGtt                         >  1:81080/1‑250 (MQ=18)
                                                                                                                           gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtat                       <  1:45936/250‑1 (MQ=35)
                                                                                                                             aCCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGCCGGTAAGGTGATATGAACCGTTATAAc                    >  1:11208/1‑251 (MQ=21)
                                                                                                                               cTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCg                  >  1:443838/1‑251 (MQ=34)
                                                                                                                                 gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACCCGATGATGAATCGAAAGAAACATCTTCGGGTTGTTAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGc                <  1:377819/251‑1 (MQ=21)
                                                                                                                                 gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGc                >  1:577923/1‑251 (MQ=18)
                                                                                                                                 gCAAAGGTTCACTGCGAAGGTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGc                <  2:109508/251‑1 (MQ=18)
                                                                                                                                  cAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCg               <  1:329116/251‑1 (MQ=18)
                                                                                                                                   aaaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTACGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGGGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGGGATATGAACCGTTATAACCGGCgc              >  2:18345/1‑250 (MQ=16)
                                                                                                                                      ggTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGAttt           <  1:134691/251‑1 (MQ=18)
                                                                                                                                       gTTCACTGCGAAGTTTTGCTCTTGAAAAATCTGGATCAAGCTGAAAAGTGAAACACTGAACAACGAAAGTTGCTCGCGAGTCTCTCAAATTGTCGCAACACGATGATGAATCTTAAGAAACATCTTCGTGTTGTGAGGTTAAGCGACGAAGCGTACACGGTGGAGGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGTTAAGCGTCGGTAAGGTGATATGAACCGTAATAACCGGCGATTTc          <  2:70944/251‑1 (MQ=11)
                                                                                                                                         tCACTGCGACGTTTTGCTCTTTAAAAATCTGGATCACGCAGACCATTGAAACACTGCACAACGAAAGTTGTTCGTGAGTCTCTCACATTTTCGCAACACGATGCGGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGGGGATGCCCTGGCAGCCAGAGGCGATGAAGGACGTGCTAATATGCGATACGCGTCGGTAAGGGGATATGAACCGTTATAAccaccgctttc          >  2:662001/1‑240 (MQ=11)
                                                                                                                                           aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt                                                                                           <  1:54081/166‑1 (MQ=18)
                                                                                                                                           aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt                                                                                           >  2:54081/1‑166 (MQ=18)
                                                                                                                                              gCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGTTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGCTGAAGGCCGGGCTAATCTGCGATAATCATCGGTAAGGTGATATGAACCGTAAGAACCGACGATTTCCAAATgg   >  2:546665/1‑251 (MQ=11)
                                                                                                                                               cGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATggg  <  1:213302/251‑1 (MQ=32)
                                                                                                                                                   gTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAGATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTa                                          <  2:109481/207‑1 (MQ=11)
                                                                                                                                                   gTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAGATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTa                                          >  1:109481/1‑207 (MQ=11)
                                                                                                                                                                                                                                                        |                                                                                                                                                 
ACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGG  >  NC_000913/4037248‑4037641

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