Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3690CT24.0% 111.5 / 15.0 50E197K (GAA→AAA) yghSputative ATP‑binding protein YghS
Reads supporting (aligned to +/- strand):  ref base C (4/34);  new base T (12/0);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 1.50e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCTT  >  NC_000913/3133122‑3133608
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
cATCACCCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                                                                                                                                                                                  <  2:91693/251‑1 (MQ=255)
   cacaCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTaa                                                                                                                                                                                                                                              <  1:453727/251‑1 (MQ=255)
              accGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCCAGGATCTTTTCGCCCTTTAATGTCAGTTGCGGCGGAACGGGGATTTTTTTCTGTAACGCTGCCAGtt                                                                                                                                                                                                                                   >  1:58668/1‑251 (MQ=255)
                    agcAGCCCCGGAATGGAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTggtggt                                                                                                                                                                                                                            <  2:498182/249‑1 (MQ=255)
                      cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTc                                                                                                                                                                                                                           <  2:232647/250‑1 (MQ=255)
                      cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGGGCAGCCATTCGTTAAATCAGCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTc                                                                                                                                                                                                                           <  2:242938/250‑1 (MQ=255)
                                                   cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcatc                                                                                                                                                                                            <  2:58668/251‑1 (MQ=255)
                                                   cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTCAAGGATCTTTGGGCCCTTGTATTTCCGTTGGGGCGTAACGGCGGTTTTTTTCTGTTACGCCTCCAGTTTGTGGGTGGGTTTACGCGCAAACCCGGtgtttttatt                                                                                                                                                                                            >  2:4809/1‑242 (MQ=255)
                                                          aCTATTAACAGTTTGAGTTTGTAGTGGATCCATTCGTTAAATCATCCCCTATCCATTTTGGTTATCAGGAAAGCGCGGCGTGAATTCCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTTGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCa                                                                                                                                                                                     <  2:119514/251‑1 (MQ=255)
                                                           cgtttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCaa                                                                                                                                                                                    <  2:35786/248‑1 (MQ=255)
                                                             aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGATTTTTTTCTGTAAAGCTGCCCGTTGGTGGTCAGGGTTACGCGCAAACGAGGGTTTTTCATCAATGCCa                                                                                                                                                                                     >  1:316896/1‑248 (MQ=255)
                                                                aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAGCTTTTCGCCTTTGAATGTCAGTTGCGTCGTAACGGCGATTTTTGCCTGTAACGCTTCCAGTTGGTGGTCAGTTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACg                                                                                                                                                                                <  1:596656/250‑1 (MQ=255)
                                                                    ggttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCTTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCTGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTTGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCATGACGAATca                                                                                                                                                                          <  1:379007/249‑1 (MQ=255)
                                                                       tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATACAGTTCCAGGAGTTTTTCGGCATTTAAATTCAGGTGCGGCGTAAAGGCGATTTTTTTCTTTAAAGCTGGCCGTTGGGGGGTAGGTGTACGCGCCAAAACGGTTTTGTCATTAATGCCAAGACGAaacaac                                                                                                                                                                        >  1:130809/1‑250 (MQ=255)
                                                                              gTAGTGCTTCCATTCGTTAAATCATCCCCGATCCAGTTTGGTGATCAGGAAAGCGCGGCGTGTGTTGCGCGTAGTGACGCTTGCAGAATTTCAGCGGCGGGGTGTCGCCCATCCAGTTCAAGGAGTTTTGCGCCAGTGAATGTCAGTTGCGGCGTAACGGCGATTGTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATTCggg                                                                                                                                                                <  2:427068/251‑1 (MQ=255)
                                                                                   gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCaaa                                                                                                                                                            <  1:575347/250‑1 (MQ=255)
                                                                                         aTTCGTCAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCCTTGAATGTCAGATGCGGCGGAACGGCGATTTTTTTCTTTAACGCTGTCAGTTTGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCCCGACGAATCAACAATAAGGGGAAATAAAATg                                                                                                                                                     >  1:154438/1‑251 (MQ=255)
                                                                                            cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGAACTTTGCGCCATTTGATGTCAGTTGCGGCGGAAAGGGGAATTTTTTCTGTAAAGCTGGCCGTTGGGGGGTCGGTTTAAGCGCCAAAGCGGGTTGTTCATTAATTCCAAGACGAAATAACAAAACAGGCAAATAAGATGcc                                                                                                                                                   >  1:550542/1‑250 (MQ=255)
                                                                                             gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGGCCATCCAGTTCAAGGATCTTTGCGCCCTTTGATGTCAGTTGCGGCGTAACGGCGATTTTTTTCTTTAACCCTGCCCGGTGGGGGTCAGGTTTTCGCGCAAACCCGGGTTGTTTTTCCATTCCAAGACGAAATCAAAAAAACGGGaaaaaaaaagcca                                                                                                                                                   >  1:402061/1‑241 (MQ=255)
                                                                                               tAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGccatc                                                                                                                                                <  1:597797/250‑1 (MQ=255)
                                                                                                 aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGccatccat                                                                                                                                             <  1:224822/251‑1 (MQ=255)
                                                                                                  atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTAGCCCATCCCGTTCAAGGGAATTTTCGCCCCTTAATGTTAGTTGCGGCGTAAAGGCGATTTTTTTCTGTAACGCTGCCAGTTGGGGGGCAGGGTTATGCGCAAACGCGGGTTTTTTATCAATGCCCAGACGAAACAAAAAAACGGGCAAATAAAATGCCACCCAtt                                                                                                                                            >  1:585825/1‑251 (MQ=255)
                                                                                                       ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTAc                                                                                                                                       <  1:51872/251‑1 (MQ=255)
                                                                                                                aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAAGATCTTTGCGCCATTGATTGGCAGTTTCGGGGTAACGGGGCTTTTTTTCTGTAACGCTTCCAGTTGGGGGGTAGGTTTACGCGCCAACGCGGGTTGTTCATCAATGCCCAGAAGAAACAACAAAAAGGGGAAAAAAAAAGCCATCCATTGGTaaacctttaac                                                                                                                              >  1:105964/1‑241 (MQ=255)
                                                                                                                  tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGGTCTTTGCGGCATTGAATGTCAGGTGCGGCGGAACGGCGGTTTTTTTCTGTAACGCTGCCCGGTGGTGGTGAGGTTTAAGCGCAAAAGCGGGTTGTTTATTAATGGCAAAAAGGATCAACAAAAAGGGGAAATAAAATGGCATCCATTGGTACAACtttcgctc                                                                                                                            >  1:107704/1‑245 (MQ=255)
                                                                                                                   ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTc                                                                                                                            <  1:402422/250‑1 (MQ=255)
                                                                                                                   ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGCTCAAGAATTTTTTCGCCCCTGAATGGTCGTTGCGGCGGAACGGCGGTTTTTTTCTGTCAACCTGGCCGCTGGTGGTCAGGGTTCCGCGCAAAACCGGGTTGTTTATCCAAGCACCGACACAACTCCAATACCGGCAAATAAGAAGCCCTCCCTTGGGACAGCTTCAGCTcc                                                                                                                           >  1:587326/1‑250 (MQ=255)
                                                                                                                        gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTg                                                                                                                      <  1:407943/251‑1 (MQ=255)
                                                                                                                           cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCTGCGTAACGGCGATTTTTTCCTGTAACGCGGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGGTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGTTACAGCTGCAGCTCGCGCTGCCt                                                                                                                   <  2:591420/251‑1 (MQ=255)
                                                                                                                           cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCCCCCCGGTCCAGGGTCTTTTCGCCATTGTATGTCAGTTTCGGCGTAACGGCGGTTTTTTTCTTTTACGCTGCCAGTTGGGGGGCGGGGTTTCGCGCCAAAGCGGGTTTTTTATTCAATCCAAGACGGATCAACAATACGGGCAAATAAGATGCCATCCCTTTGTACAGCTTcctctcccgccgcct                                                                                                                   >  2:248612/1‑237 (MQ=255)
                                                                                                                              gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTaa                                                                                                                <  1:228921/251‑1 (MQ=255)
                                                                                                                                            aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAAc                                                                                                  <  1:367201/251‑1 (MQ=255)
                                                                                                                                             aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAAcc                                                                                                 <  1:408878/251‑1 (MQ=255)
                                                                                                                                                gtgcggaggggcgctatcAGAATTGCATCGGCGGGGTGTCGCCCATCCAGTTCAATGATGTTTGCGCCATTGAATGTTAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTTGTGTTGCGTGTTACGCGCAATCGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGTTGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAGCATTTTTATCCAACCGt                                                                                               <  2:453556/233‑1 (MQ=255)
                                                                                                                                                gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGtt                                                                                              <  1:105678/251‑1 (MQ=255)
                                                                                                                                                   cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTAcc                                                                                           <  1:527705/251‑1 (MQ=255)
                                                                                                                                                       gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGccc                                                                                       <  1:248612/251‑1 (MQ=255)
                                                                                                                                                                             ggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGc                                                                 <  2:414839/251‑1 (MQ=255)
                                                                                                                                                                              gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGcc                                                                <  1:42970/251‑1 (MQ=255)
                                                                                                                                                                                                           cTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTCATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCAATTCAACCt                                   >  2:259759/1‑251 (MQ=255)
                                                                                                                                                                                                             ttGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATGTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCaa                                      <  1:126453/246‑1 (MQ=255)
                                                                                                                                                                                                              tGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCACCAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGACTTAACATTTTTATCCAACCGTTACCACCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCg                                >  1:490307/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGc                          >  1:453822/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGtgatga              <  1:14210/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcggtgTAACGGCGATTTTTGCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTGTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGc            <  2:130809/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cggcgTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCa           <  1:4809/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGGCGATTTGTTCCTGTAACGCTGCCAGGTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTGCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATGTTTATCCAACCGTTCCCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAg      <  2:556504/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTCATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  <  1:259759/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  >  1:436031/1‑248 (MQ=255)
                                                                                                                                                                                                                                                gATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  <  2:436031/248‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
CATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCTT  >  NC_000913/3133122‑3133608

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