Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913106,1220AG25.0% 80.2 / 10.2 44D273G (GAT→GGT) ftsZcell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base A (8/25);  new base G (11/0);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 9.86e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCAGCAGCCAGTGATGGATCGCTACCAGCAGCATGGGATGGCTCCGCTGACCCAGGAG  >  NC_000913/105876‑106354
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
gcgcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACAGCACGTCGGGCTTCGACCTGCGTCTGGATGAg                                                                                                                                                                                                                                        <  2:170661/251‑1 (MQ=255)
  gcTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGtt                                                                                                                                                                                                                                      <  1:100146/251‑1 (MQ=255)
     gtgCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTc                                                                                                                                                                                                                                     <  1:131730/249‑1 (MQ=255)
        cAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCCGCGTCTGGATGAGTTCGAAAc                                                                                                                                                                                                                                >  1:330776/1‑251 (MQ=255)
           ggTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACggt                                                                                                                                                                                                                             <  1:301795/251‑1 (MQ=255)
              aTCGCTGAACTGAGTACGCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGTAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTTGAAGATATCGACCGGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACggtagg                                                                                                                                                                                                                          <  1:414466/251‑1 (MQ=255)
                   tGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCCTCTACCTTCGTCTGGAGGAGTTCGTAACgggtgggtaa                                                                                                                                                                                                                      >  2:323125/1‑242 (MQ=255)
                    gAACTGATTACTCGGCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGAGGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGTCCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAacac                                                                                                                                                                                                                   <  2:226650/251‑1 (MQ=255)
                                      ggtgtgatgaACTTGTACTTTGCAGACGTACTCACCGTAATGGCTGAGATGGGCGACGCAAGGATGGGTTCTGGCGTGGCTAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATTGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCAGGCTTCGACCTGCGTCTGTATGAGGTCGAAACGGTAGG‑TAACTCCATCCGTG‑CATTTGCTTc                                                                                                                                                                                                <  2:177339/247‑1 (MQ=255)
                                      ggTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGTTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTc                                                                                                                                                                                                <  1:184223/251‑1 (MQ=255)
                                       gTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCCCCGTCTGGGTGGGTTTGAAACGGGAGG‑GAACACCCTCCCTG‑GTTTTGCTTcc                                                                                                                                                                                               >  2:102985/1‑251 (MQ=255)
                                              gaACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAAcg                                                                                                                                                                                        <  1:352233/251‑1 (MQ=255)
                                                   tGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGGTTAGGTTGAAACGGGAGG‑TAACACCCTCCGGG‑CATTTTCTTCCGACAAcgcgcct                                                                                                                                                                                   >  1:298346/1‑248 (MQ=255)
                                                     gACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGGCTCGACCTGCGTCTGGGTGGGTTCGGAAAGGGAGG‑GAACACCCTCCCGG‑CATTTGCTTCCGACAACGCGACttg                                                                                                                                                                                 >  2:84715/1‑249 (MQ=255)
                                                      aCTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGCCTTCGCCCTGCGGCTGGGTGGGTGTGAAAAGGTAGG‑GAAAAACATTCCTG‑CATTTTCTTCCCACAACGCGTCtgtg                                                                                                                                                                                >  2:305553/1‑251 (MQ=255)
                                                       cTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATTAGTTCGAAACGGTAGG‑TAAAACCATCCGTG‑CCTTTTCTTCCGACaaagccaatgtgg                                                                                                                                                                               >  1:159403/1‑240 (MQ=255)
                                                                 gTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGGAGGGTTCGAAAAGGGAGG‑TAACACCATCCGGG‑CATTTGCTTCCGACAACGCGACTGGGGTGATTGGTAc                                                                                                                                                                     >  1:165669/1‑251 (MQ=255)
                                                                       aCCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTctct                                                                                                                                                               <  2:113458/251‑1 (MQ=255)
                                                                        ccGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCtt                                                                                                                                                              <  1:112109/251‑1 (MQ=255)
                                                                               gTCTGAGATGGGCTACGCAAGGATTGGTTCTTGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAATATATCGACCTGTCTGGCGCGCGCGGAGTGCTGGTTAACAGCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCTAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGGCAACGCGACTGTGGTTAGCGGTACTTCTCTTGACCCgg                                                                                                                                                       <  2:278706/251‑1 (MQ=255)
                                                                                      aTGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑GAAAACCATCCGTG‑CATTTTCTTCCGACAAAGCGACTTGGGTTTTCGGTACTTCTCTTTACCCGGATatga                                                                                                                                                 >  1:319833/1‑250 (MQ=255)
                                                                                                  atgatgGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGACCTCGACCTGCGGCTGGGTGAGGTTGAAACGGGAGG‑TAACAACCTCCGTGGCATT‑GCTTCCGACCAAGCGACTGTGGTTATCGGTACTTTTCTTTACCCGGGTATGAATGACGAGCcgcg                                                                                                                                    >  1:282412/1‑251 (MQ=255)
                                                                                                       gggTTCTGGCGAGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCCCGGCGGGCTTCGCCCCCCGTCTGGATGAGTTCGAAACGGGAGGGTAA‑ACCCTCCGCG‑CCTTTTCTTCCCACAACCCGACTGTGGTTATCGGGACTTTTCTTTACCCGGGTATGGATGACCAGGTGcgcggt                                                                                                                                >  2:410805/1‑248 (MQ=255)
                                                                                                                 gtggcgAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGGTGAGGTTGCAACGGTAGG‑GAACACCATCCGTG‑CCTTTGCTTCCCACAAACCGGCTGTGGTTATCGGTACTTCTCTTTACCCCGGTATGAATGACGAGCTGCGCGGAACCGtttttgt                                                                                                                     >  2:368475/1‑250 (MQ=255)
                                                                                                                    gcgAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTTGATGAGTTCGAAACGGGAGG‑TAAAACCATCCGTG‑CATTTTCTTCCGCCAACGCGACTGTGGTTATTGGTAATTCTCTTTACCCCGAAATGAATTACGAGCTGCGCGTAACCGTTGTTGCGAc                                                                                                                  >  2:311139/1‑251 (MQ=255)
                                                                                                                          ggTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGGCTGGGAGGGTTTGGAACGGGAGG‑GAAAACCATCCCTG‑CCTTTTCTTCCGACCACGCGACTGTGGGTATCGGTACTTCTCTTTACCCCGATATGAATGACGAGCTGTGCGTAACCGTTGTTTCGACAgggat                                                                                                            >  2:198529/1‑248 (MQ=255)
                                                                                                                                   cGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGGCCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGAGATGAATGACGAGCTGCGCGTAACCGTTGTTGCGAGAGGTATCGACATGGa                                                                                                   <  2:272079/251‑1 (MQ=255)
                                                                                                                                        ggAAGAAGCTGCTGAAATGGCTATCTCTTCGCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAAc                                                                                              <  2:176756/251‑1 (MQ=255)
                                                                                                                                              agCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTg                                                                                        <  2:28880/251‑1 (MQ=255)
                                                                                                                                                   ctgAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACa                                                                                                                 <  2:349117/221‑1 (MQ=255)
                                                                                                                                                   ctgAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGGAGG‑GAAAACCATCCGTG‑CATTTTCTTCCGACAACGCGGATGTGGGTATCGGTACTTCTCTTGCCCCGGATATGAATGACGAGCTGCGCGTAACCGGTGTTGCGACa                                                                                                                 >  1:349117/1‑221 (MQ=255)
                                                                                                                                                              aTCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACGCTGGTGa                                                                         <  2:319833/250‑1 (MQ=255)
                                                                                                                                                                       ccGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATaa                                                                  <  1:343176/248‑1 (MQ=255)
                                                                                                                                                                       ccGCTGCTGGAAGATAGCGACCGGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAg                                                                 <  2:276570/249‑1 (MQ=255)
                                                                                                                                                                              tgGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGGAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGGTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGGAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGGGGCCAATAAACAGGTTc                                                          >  1:43468/1‑249 (MQ=255)
                                                                                                                                                                               ggAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTcag                                                        <  2:298346/250‑1 (MQ=255)
                                                                                                                                                                                  aGATATCGACCTGTCTGGGGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACTCCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCGCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCagcagc                                                    <  1:232229/251‑1 (MQ=255)
                                                                                                                                                                                     tatCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCTACCTGCGTCTGGGTGAGTTCGGAACGGGAGG‑TAAACCCCTCCGTG‑CCTTTGCTTCCGACAAACCGACTGTGGTTTTCGGTACTTCTCTTGACCCCGAAATGGATGACGAGGTGCGCGTAACCGTTGTTGGGACAGGGATTGGCATGGACAAAACTCCTGAAATCACTCCTGTGGCCAAAAAGCAGGTTTAGAAGcc                                                   >  2:234608/1‑249 (MQ=255)
                                                                                                                                                                                      atCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCAGCAGCCAGt                                                <  1:183602/251‑1 (MQ=255)
                                                                                                                                                                                      atCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCAGCAGCCAGt                                                <  1:311139/251‑1 (MQ=255)
                                                                                                                                                                                         gACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGGGTGGGTGAGTTCGAAACGGTAGG‑TAAAACCCTACGTG‑CATTTGCTTCCGACAAAGCGACTGTGGTTTTCGGTCCTTCTCTTGAACCGGATATGAATTACGAAGTGCGCGTAAACGGTGTTTCGACAAGGATAGGCAAGGGCAAACGGCCTGAAATCAATATGGGGACCCATACGCAGGTTtagcacgccgaggt                                             >  1:109643/1‑242 (MQ=255)
                                                                                                                                                                                                              gTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCAGGATATGAATGACGAGCTGCGCGTATCCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCATGTGGAGCAGCTAGTGATGGATCGGTACCAGTGGCATgg                        <  2:274394/251‑1 (MQ=255)
                                                                                                                                                                                                                gCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGGTGAGTTCGAAACCGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGGTATCGGTACTTCTCTTGACCCGGGGATGAATGACGCGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGGCAAAACTCCCGAAATCACTCTGGTGGCCAAAAAACCGGGTTAGCAGCCCGTGATGGATCGCCACCCCCAGCATgggg                      >  2:55612/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCAGCGGCCAGTGATGGATCGCTACCAGCAGCATGGGATGGCTCCGCTGACCCAGGAg  <  1:221852/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
GCGCTGTGCAAGGTATCGCTGAACTGATTACTCGTCCGGGTTTGATGAACGTGGACTTTGCAGACGTACGCACCGTAATGTCTGAGATGGGCTACGCAATGATGGGTTCTGGCGTGGCGAGCGGTGAAGACCGTGCGGAAGAAGCTGCTGAAATGGCTATCTCTTCTCCGCTGCTGGAAGATATCGACCTGTCTGGCGCGCGCGGCGTGCTGGTTAACATCACGGCGGGCTTCGACCTGCGTCTGGATGAGTTCGAAACGGTAGG‑TAACACCATCCGTG‑CATTTGCTTCCGACAACGCGACTGTGGTTATCGGTACTTCTCTTGACCCGGATATGAATGACGAGCTGCGCGTAACCGTTGTTGCGACAGGTATCGGCATGGACAAACGTCCTGAAATCACTCTGGTGACCAATAAGCAGGTTCAGCAGCCAGTGATGGATCGCTACCAGCAGCATGGGATGGCTCCGCTGACCCAGGAG  >  NC_000913/105876‑106354

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: ‑