Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1320GA25.4% 126.0 / 14.9 64A16V (GCT→GTT) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (12/35);  new base A (16/0);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 8.30e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCG  >  NC_000913/4003885‑4004370
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
atGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCg                                                                                                                                                                                                                                                 <  2:69000/250‑1 (MQ=255)
                        aaaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGTCAAGAGCGAGTAATACGCCCTGc                                                                                                                                                                                                                         <  1:466200/250‑1 (MQ=255)
                           cAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTGAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGAGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgcg                                                                                                                                                                                                                    <  1:346473/251‑1 (MQ=255)
                            aGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAAGGCTCATCAGCACCACCATAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACGCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                                                                                                                                                                                                                   <  1:273781/251‑1 (MQ=255)
                                           ccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                                                                                                                                                                                                    <  2:96947/251‑1 (MQ=255)
                                              aTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg                                                                                                                                                                                                 <  1:273234/251‑1 (MQ=255)
                                                  gCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                                                                                                                                                                                                                   <  1:437405/214‑1 (MQ=255)
                                                        cAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTGAAGTAAGCTGGCGCCCTTCCCCCCAAAAAAAAAACGGCAAGAGCGGGGAATAAGCCCCGCCGCGGTTGTTTTTCCTCCATCGGGGATACTcct                                                                                                                                                                                       >  2:590979/1‑249 (MQ=255)
                                                        cAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCTCGTAAATTAACTTGAAGTACACCTGCGCTATACCCCCAAAAAAAAAAACGGCCAGAGAGAGTAAAAAGCCCTGCCGCGTTTGTTTTTcatccacccgggaaactcct                                                                                                                                                                                       >  1:372247/1‑237 (MQ=255)
                                                         agagACAGCCCGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAAGAAATCAACTTGAAGTACGCTGGCGCTATTCCCCCCAAAAAATAAACGGCAAGAGGGAGGAATACGCCCCGCCGCGTTTGTTTTGATTCCATCCGGGATACTcctt                                                                                                                                                                                      >  1:148856/1‑248 (MQ=255)
                                                             aCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCCGGCGCTATACCCCAAATAAAAAAAACGGCAAAAGCGAGGAATACGCCCTGCCGCGGTTGTTTTGTTTCCATCCGGAAAACTCCtttttt                                                                                                                                                                                  >  1:133986/1‑250 (MQ=255)
                                                                              tttCTGTGGCGGCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACATCAGtt                                                                                                                                                                 <  2:30081/250‑1 (MQ=255)
                                                                                catgtgGCGTCGGAATCAGCGTTGTTAAATAGGACCACGGGCGTCAAATGCGCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTACATCGGCGTGCACGTAGTAAATCAGCTTGAACTACGCTGGCGCGATACCCCAAATAAAATAAGCGGCAAGTTCGAGGAATACGCCCTGCAGCGTTTGTTTTGGATCCATCCCGAATACGTATTTTTAATTTTTAACAGCATTTTAc                                                                                                                                                              <  2:255394/249‑1 (MQ=255)
                                                                                   gtgGCGGCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCt                                                                                                                                                            <  2:148856/250‑1 (MQ=255)
                                                                                    tgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGACTTCATCGGCGGGCACGTAGTAAATTAACTTGAAGTACGCTGGCGCTCTACCCCAAATAAAATAAACGGCAAGAGCGAGGAATACGCCCTGCCGCGGTTGTTTTTCATCCATCGGGAATATTCATTTTTTATTGGGAACAGCAGTTTATCTGc                                                                                                                                                          >  1:460839/1‑251 (MQ=255)
                                                                                     ggCGTCTGAAACAGCGTGTTTAAATAGGACCACTGGCGGCAAATGCGCATCAGCACCACCATAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTGCAGCGGCGGGCACGTAGTAAATCAACGTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCTCGTTTGTTTTGCATCCAGCGGGAAGACTCATGTTTAATTGGTAACAGCAGGTTATCTGCt                                                                                                                                                         <  2:35625/251‑1 (MQ=255)
                                                                                        gTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtttt                                                                                                                                                                                    <  2:16816/222‑1 (MQ=255)
                                                                                        gTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAGACTCATTTTTAATTGGTAACAGCAGTTTACCTGCtttt                                                                                                                                                      <  2:506698/251‑1 (MQ=255)
                                                                                        gTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACCACTTGAAAGACGCTGGCGCTTTACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTTTATCCCTCGGGGATAATCATTTTTGATTggaagaggacccgtccgaaccccc                                                                                                                                                      >  1:16816/1‑229 (MQ=255)
                                                                                            gAATCAGCGTTTTTAAATAGGACCACAGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTAt                                                                                                                                                    <  1:256006/249‑1 (MQ=255)
                                                                                             aaTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTc                                                                                                                                                 <  2:68617/251‑1 (MQ=255)
                                                                                                cAGCGGTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTGTATttcttc                                                                                                                                              <  2:262661/251‑1 (MQ=255)
                                                                                                  gCGTTTATAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTTCCTGCTTTTATGTCTTCaa                                                                                                                                            <  1:459336/251‑1 (MQ=255)
                                                                                                                  ccccTGGCGGCAAATGCGCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCAGCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAAGACGCCCTGCCGCGTTTGTTGTGGATCCATCGGGAATACTCATTTTTAAGTAGTAACAGCATTTTACCTGCTTTTATGTCTTCAACCTACCATATAggtgg                                                                                                                           <  1:595024/248‑1 (MQ=255)
                                                                                                                       ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGTGAATACTCAGTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAggtggcggt                                                                                                                       <  1:282505/250‑1 (MQ=255)
                                                                                                                          ggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCa                                                                                                                   <  2:109167/251‑1 (MQ=255)
                                                                                                                          ggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAAACAAATTTAAGTAAGCTGGCGGCATACCCCAAATAAAAAAAACGGCAAGAGCGAGGAATAAACCCTGCCCCGTTTGTTTTGCAACCAACGGGGATACTTATTTTTTATTTGTAACAGCAGGTTACCCGCTTTTATTTTTTCAACCTACCATATAggtggggggggg                                                                                                                    >  1:11538/1‑243 (MQ=255)
                                                                                                                               aTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGc                                                                                                              <  1:448102/251‑1 (MQ=255)
                                                                                                                                 gCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGCCGAAGGAAAACAACTTGAAGAAAGGGGGGGCCCGACCCCAAAAAAAATAAACGGAAAGAGGGGGAGAAAAGCCCCGACGCGGTTGTTTTTCATCCATCGGGGATAATTAATTTTAATTGGTACCAGCAGTTTACATGCTTTTATGTCTTCAAACTACCaacaaggggggggggggcagggca                                                                                                              >  2:163613/1‑227 (MQ=255)
                                                                                                                                       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATGTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATaaa                                                                                                      <  1:517866/251‑1 (MQ=255)
                                                                                                                                            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGCAGGAAAACAACTTGAAGTACGGCGGCGCCATACCCCCCATAAAAAAAACGGCCAGAGCGAGTAATACGCCCCGCCGCGGTTGTTTTTCATCCCTCGGGGAAACTACTTTTTTATTTGTTACCGCCGTTTTCCTGCCTTTATTTCTTCAAACTACCATAATagtggcggtgggccttgcactataaacatgt                                                                                                 >  2:421111/1‑231 (MQ=255)
                                                                                                                                              accaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGCGTAAAACAAAATGAAGAACGCTGGGGCTCTATCCCAAATAAAATAAACGGGCAGAGAGGGTAAAACGCCCTGCCCCGGTTGTTTTTCATCCATCGGGGAAACTCATTTTTTATTGGTTACAGCAGGTTATCTGCTTTTTTTTCTTCAACCTACCATATAGGTGGGGGGGGCACTGGCTATATAAAGCTGttt                                                                                               >  2:33329/1‑250 (MQ=255)
                                                                                                                                                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGttctt                                                                                             <  1:281206/251‑1 (MQ=255)
                                                                                                                                                caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATAAACTTGAAGTACGCTGGGGGTATAACCCAAATAAAATAAGAGGGAAAAGAGAGGAATACACCCTGCCGCGTTTTTTTTTGATCCCTCGGGGATACTCATTTTTTATTGGTAAAAGCAGTTTAACTGCTTTTTTTTTTTCAACCTAACCTATAGGGGGGGGGGGCACTGGCAAAATaaaaaggttttt                                                                                             >  1:428833/1‑242 (MQ=255)
                                                                                                                                                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATtgt                                                                                       <  1:326610/251‑1 (MQ=255)
                                                                                                                                                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGAAAATCAACTTGAAGTAAGCTGGCGCTATTCCCCAAATAAAATAAACGGCCAGAGCGGGGAATACGCCCTGCCGCGTTTGTTTTTcatcca                                                                                                                                                                                                     >  1:376310/1‑141 (MQ=255)
                                                                                                                                                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                                                                                                                                                                                                     <  2:376310/141‑1 (MQ=255)
                                                                                                                                                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGGAAAACAACTTGAAGTAAGCCGGCGCTTTACCCCCAATTAAATAAGCGGCAAGAGCGAGGAATACGCCCCGCCGCGGTTGTTTTTCTTCCCTCCGGAAGACTCCTTTTTTATTGGTAACAGAAGTTTACCTGCTTTTTTTTCTTCTACCCACCATATATGTGGTGGGGGGGCTGGCAATATAAAACTTTTCTTTCTTgtgtt                                                                                    >  2:440825/1‑250 (MQ=255)
                                                                                                                                                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGCAGGAAAACAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAAAAAACGGCAAGAGCGAGGAAAACGCCCTGCCCCGGTTTTTTTTCATCCATCGGGGATACTCATTTTTTATTTGGTACAGTAATTTTACTGCTTTTTTGTCTTCAACCTACCATATATGTGGGGGTGGCGCTGGCACTATAAAGCTGTTTTTCATTttgtg                                                                                    >  2:10601/1‑250 (MQ=255)
                                                                                                                                                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACCCAGAAAATCAACTTGAACAACGCAGCCGCTAAACCCCAAAATAAAAAACCGGCAAGAGAGGGTAAAAAAACCTGCACAGGTTTTTTTTAATCCATCAGGAATAATACTTTTTTATTTGTAAAAACAGTTTAACCGCTTTTATGTATTCCAACTACCCTATTTGTGGCGGTGGCCCTGGCAATATTAAACTGTTCTTCATTTTGTa                                                                                    >  1:274408/1‑250 (MQ=255)
                                                                                                                                                                cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAAc                                                                             <  1:419771/251‑1 (MQ=255)
                                                                                                                                                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAAACAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGTATCCATCGGGAATACTCATTTTTTATTGGTAACAGCAGTTTACCTGCTTTTTTGTCTTCAACCTACCCTATAGGTGGCGGTGGGACTGGCAATATAAAACTGTTCTTCATTGTTTAATTCAAcc                                                                            >  1:168520/1‑250 (MQ=255)
                                                                                                                                                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAAAAAAAAAAACGGCAAGAGCGAGGAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAAAAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGGGGCACTGGCAATATAAAACTGTTCTTCATTGTGTAATTCAAcc                                                                            >  1:235555/1‑250 (MQ=255)
                                                                                                                                                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGc                                                                                                    <  2:517660/214‑1 (MQ=255)
                                                                                                                                                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAGTCAACGCGggcgacggcga                                                               <  1:443028/251‑1 (MQ=255)
                                                                                                                                                                              cGTCAGGATTTCATCGGCGGACCCGTAGTAAATCAAATTGAAGGAAGGTTGCGCTTTACCCCCAATAAAAAAAAAGGCAAGAGCGAGGAAAACGCCCTGCCGCGTTTGTTTTGTATCCATCGGGGATACTCATTTTTAATTTGTAACAGCAGTTTACCTGCTTTTATGTCTTCAAACTACCATATTCGTGGGGGTGGGCCTGGGAATATAAAGc                                                                                                    >  1:517660/1‑214 (MQ=255)
                                                                                                                                                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCCAATAAAATAAACGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTTTGTCTTCAACCTACCATATAGGTGGCGGTGGCACTTGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggggacggcgac                                                              >  1:113589/1‑251 (MQ=255)
                                                                                                                                                                                            cggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGc                                                  <  2:74159/250‑1 (MQ=255)
                                                                                                                                                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGCGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCAc                                                >  2:398382/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAAc                                           >  1:207078/1‑249 (MQ=255)
                                                                                                                                                                                                             aTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACtagt                                >  2:217329/1‑251 (MQ=255)
                                                                                                                                                                                                               cAACTTGAAGTACGCTGGCGCTATAACCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACtagt                                >  1:477946/1‑249 (MQ=255)
                                                                                                                                                                                                                 ccgtGAAGGACGCTGTGGCTATACCCCAAATAAAATAAGCGGCAAGAGCTAGTAATTCGCCCTGCCGCGTTTGTTTTGCGTCCATCGGGAATACGCAGTTTTAATTTGTAACAGCAGTTTACGTGCTTTTATGTCTTCAACCTACCATATATGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAGTTCCTGCACGCAACAGGCTACTAGTAGCa                            <  2:197545/248‑1 (MQ=255)
                                                                                                                                                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAgg                                        >  2:238391/1‑235 (MQ=255)
                                                                                                                                                                                                                     gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAgg                                        <  1:238391/235‑1 (MQ=255)
                                                                                                                                                                                                                          aCGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCCTTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATACAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACCGGTTACTCGTAGCAGTAAAACGc                   >  2:576567/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTc                 <  2:460839/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcgcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGccc             >  2:152317/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccccAAATAAAATACGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCaa                                                                                                            >  2:521057/1‑148 (MQ=255)
                                                                                                                                                                                                                                        ccccAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCaa                                                                                                            <  1:521057/148‑1 (MQ=255)
                                                                                                                                                                                                                                         cccAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTggg    <  1:431114/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcg  >  2:316039/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aaGCTGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAt                                                         >  1:323932/1‑186 (MQ=255)
                                                                                                                                                                                                                                                      aaGCTGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAt                                                         <  2:323932/186‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
ATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCG  >  NC_000913/4003885‑4004370

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