Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,161,5650AC34.3% 95.1 / 12.9 67W369G (TGG→GGG) ftsPcell division protein FtsP
Reads supporting (aligned to +/- strand):  ref base A (12/32);  new base C (20/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 5.19e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATCGACACTTCATCGCCGTTGCTC  >  NC_000913/3161314‑3161815
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                 
cGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCT‑GCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCca                                                                                                                                                                                                                                                                   <  1:158034/251‑1 (MQ=255)
          gtGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTcc                                                                                                                                                                                                                                                         <  1:62608/250‑1 (MQ=255)
                tAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGGCATCAAAACGGTTGACGTCCCACC‑GCTGTCCATTAAt                                                                                                                                                                                                                                                   >  1:265176/1‑251 (MQ=255)
                   aaCGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCCCGGGGACAACAATACGGTTGGCGTCCCCCA‑GCTGTCCCTTAATAcc                                                                                                                                                                                                                                                >  2:289169/1‑251 (MQ=255)
                      gggAAGTGCGCCCAGGAAGGCTTACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCATACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCgg                                                                                                                                                                                                                                             <  1:202835/251‑1 (MQ=255)
                         aagggcgcCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGtc                                                                                                                                                                                                                                          <  2:94955/247‑1 (MQ=255)
                                      aaGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGGCCCGCGTTCCCCCGTTCCTTGTT‑GCGCGGGGACATCAATACGGGTGGCGTCCCCCC‑GCTGGCCATTAATTCCCggggactccccccgactg                                                                                                                                                                                                                             >  2:108873/1‑236 (MQ=255)
                                      aaGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCGTCCAGCGTTCCCCCGTTCCTTGCT‑GCGCGGGGACATCAATACGGTTGACGTACCCCA‑GCTGGCCCTTAATACCCCGGTCAT‑CACCCAGACTg                                                                                                                                                                                                                            >  1:64899/1‑251 (MQ=255)
                                         gCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGata                                                                                                                                                                                                                         <  2:142759/251‑1 (MQ=255)
                                          atacgctaaaTAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACCTTCCTTGCT‑GCGCGGGGACATCAAAACGGGTGACGTCCCACA‑GATGTACATTAATACCCCGTGCAT‑CACCCAGAATGacat                                                                                                                                                                                                                        >  1:217193/8‑248 (MQ=255)
                                            gacggAAATAAACAAGCACCACCTATGGTCCGACAACCCAAACGTGATCTATCCGGCCTCTGTCTTCCGGGAACGGCATCGCGACATTCACGGTACGGATCTGGAACATTACGCCTTCAATGTGGAACGCTTGTGGCTCGGCTGCTCGTAGATTCCAGCGGTCCTACGTTCCGTGCT‑GCGCGGTGTCATCAATACGGGTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATcg                                                                                                                                                                                                                      <  1:211980/247‑1 (MQ=255)
                                                gAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCCGCGTTCCCCCGTTTCTTGCT‑GCGCGGTGACATCAATACGGGTTACGTCCCACA‑GCTGTCCATTAATACCCGGGGCAT‑CACCCCGACTGATATCGCGAc                                                                                                                                                                                                                  >  1:53904/1‑251 (MQ=255)
                                                          agcagcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCCTCCAGCGTTCCCCCGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCCCC‑GCTGGCCCTTAATACCCGGGGCAT‑CACCCCGACTGATATCGCGACTGCGA‑ATTgg                                                                                                                                                                                                       >  1:384692/1‑251 (MQ=255)
                                                            cagcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCTCGCGTACCGTCCCGCGTTCCCACGCTCCCTGCT‑GCGCGGGGACATCCATACGGTTGACGTCCCCAA‑GCAGTCCATTAATACCCGGGGCATTCACCC‑GACTGGAATTGCGAATGCGG‑ATTgggg                                                                                                                                                                                                     >  2:99830/1‑249 (MQ=255)
                                                            cagcTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCGGCCCGCGGTTCCACCGTCCCTGCC‑GCGCGGTGACATCCATAAGGGTGACGTCCCCCA‑GCTGGCCATTAATAACCGGGTCCT‑CTACCAGACCGGTATCGCGGCCGCGA‑ACTgggg                                                                                                                                                                                                     >  1:85755/1‑249 (MQ=255)
                                                                  cACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCCTCCCCGCGTACCGTCCCGCGTTTCCCCCTTCCTTTCTTGC‑CGGGGAAATCACAACGGTTGACGTCCCCCC‑GCTGTCCCTTTATACCCGGGTACC‑CACCCCGACCGATATCGCGACCGAGA‑ATTGGGGAA‑cccg                                                                                                                                                                                              >  2:337845/1‑249 (MQ=255)
                                                                       gTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCAt                                                                                                                                                                                          <  1:125316/250‑1 (MQ=255)
                                                                         ccGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACCTTACGCCTTCAATACGGCACGCTTGCGGCTCCTACGCACGTACCGCCCACCGTTCCCACGTTCCTTGAT‑GCGCGGTGACATCAATCCCGCTGACGTCCCCCA‑GATGTCCATTAATAACCCGGTCAT‑CCCACCGACTGATATCCCGACTGCGA‑AATGGCGAA‑CCAGCCATgtt                                                                                                                                                                                       >  1:224289/1‑249 (MQ=255)
                                                                           gTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCGTACCGCCCACCCTTCCCCCCTTCCTTGCC‑GCGCCGGGACATCCATACGGGTGACGGCCCCCC‑GCCGTCCCTTTATAACCGGGGCAT‑CCCCCCACCTTATATCCCGACTGCGACATTGG‑GAA‑CCCGccctcattt                                                                                                                                                                                     >  2:164612/1‑244 (MQ=255)
                                                                             caacccaaccGGTATCTTTCCAGCCCCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCa                                                                                                                                                                                   <  2:193961/251‑1 (MQ=255)
                                                                             cAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCa                                                                                                                                                                                   <  1:190415/251‑1 (MQ=255)
                                                                                cccAAACGGTAGCTTTCCATCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTTCGACTCTGGAACATTACGCCTTCAATTTGGCACTCTTGCGGCTCGTCCGCGCGTACCGTCCAGCTTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTGGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCTCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑Gtt                                                                                                                                                                               <  2:108571/251‑1 (MQ=255)
                                                                                       ggTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCGTCCAGCGTTCCCACGTTCCTTTCT‑GCGCGGTGACATAAATACGGATGACGGCCCCCA‑GCTGTACATTAATACCCGGGGCAT‑CACCCCGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGggagc                                                                                                                                                                         >  1:461971/1‑250 (MQ=255)
                                                                                       ggTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCGGCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGGCATCAATACGGCTGGCGTCCCCCA‑GCTGTCCATTCATACCCCGGGCAG‑CACCCCGACCGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCACGATTTCA‑GGTGgcagac                                                                                                                                                                        >  2:331841/1‑249 (MQ=255)
                                                                                        gTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCg                                                                                                                                                                                                                                              <  1:341964/184‑1 (MQ=255)
                                                                                        gTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGCACCGTCCCCCGTTCCCACCCTCCTTGCT‑GCGCGGGGACATCCATACCGTTGACGTCCCCCC‑GCTGTCCCTTTATACCCCGGTCCT‑CCCCCACCCTGATATCCCCAGCGCGA‑ATTTGCGAA‑CCCGCCATTATTTCA‑CTTgggcgca                                                                                                                                                                        >  2:97786/1‑245 (MQ=255)
                                                                                          aTCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCCCA‑CCTGTCCCTTTATACCCCGGGCAT‑CACCCAGACCGATATCCCGACTGCGA‑ATTGGCGAA‑ACAGCCCTGCTTTTA‑GTTGgccgccagc                                                                                                                                                                     >  2:98678/1‑244 (MQ=255)
                                                                                                             cccGGGAACGCCCTCGCGCCATTCACGTGACGGATCTGGACCATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCTGCGTCCCCCCGTACCTTGCT‑GAGCGGTGACAGCAATACGGTTGACGTCCCCCG‑GCTGCCCATTAAGACCCGGGTCAT‑CGCCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGTCAGCAAGCGCATCGGAAGActgtctgt                                                                                                                                                  <  2:269483/250‑1 (MQ=255)
                                                                                                                gggAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCCCGTACCGCCCAACGGTCCCACCTTCCTTGCT‑GCGCCGCGACCTCCCTCCGGTTGACGTCCCCCA‑GCTGTCCCTTAATACCCGGGTCAT‑CACCCCGACTGATATCGCGGCTGCGA‑ATTTGCGAA‑CCAGCCCAGATTTCA‑GTTGGCAGCAAGCGCAACGGAAGACTGTCGGTG‑Ac                                                                                                                                              >  2:319727/1‑251 (MQ=255)
                                                                                                                gggAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCGTACCGTCCAGCGTTCCCCCGGTCCTTTCT‑GCGCGGGGACATCAATACGGGTGACGGCCCCCC‑GCCGTCCCTTTATACCCCGGTCCT‑CACCCCGACTGATATCGCGACCTCGG‑ATTTGCGAA‑CCAGCCATGATTTCATGTGGGC‑GCAAGGGCCTCGGAAGActgtcttgtgc                                                                                                                                               >  2:380179/1‑246 (MQ=255)
                                                                                                                   aaCGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGAGTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GGTGGCAGCAAGCGCATCGGAATACTGTCTGTG‑ACCAg                                                                                                                                           <  2:242340/251‑1 (MQ=255)
                                                                                                                    ccGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCCCA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACGGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGc                                                                                                                                          <  2:367406/250‑1 (MQ=255)
                                                                                                                       gCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGGGACATCAATACGGGTGACGTCCCACA‑GCTGTCCCTTAATACCCGGGGCAT‑CACCCCGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGc                                                                                                                                       >  1:122094/1‑251 (MQ=255)
                                                                                                                         aTCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGGTGACGTCCCCCA‑GCTGGCCATTAATACCCGGGTCAT‑CAACCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTTC‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGc                                                                                                                                       >  1:15806/1‑249 (MQ=255)
                                                                                                                              gcCATTCACGTTAAGGATCTGGAACATCACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCCAACCGTCCCGCGCTCCCACCTTCCTTGCC‑GCGCGGGGACATCAATCCCGTTGCCGGTCCACAAGCGGTCC‑TTCATACACGGGTCAT‑CACCCCGACTGATATCGCTCCTGCGA‑ATTGGCGGA‑ACAGGCATGATTTTA‑GTTGGGAGCCAGCGCATCGGGAGAATGTCTGTG‑AGCAGCGGccggaagg                                                                                                                                >  2:439315/1‑244 (MQ=255)
                                                                                                                                   ccACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTGGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGtt                                                                                                                           <  2:265176/250‑1 (MQ=255)
                                                                                                                                    cACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGATCCTCCAGGGTTACCACGATCCTTTAT‑GCGCGGTGACCTCAATACGGTTGACCTACCCCA‑GCTGTCCATTAATACCCGGGGCAT‑CACCCAGACTGATATCGCAACTGCGG‑CTTAGCGGA‑CCAGCCATGATTTCA‑GTTGGCCGCAAGCGCATCGGAAGACTGTCTGTG‑GCCCGCGGCAGAAGGCCGGtgg                                                                                                                          >  1:45719/1‑249 (MQ=255)
                                                                                                                                       gTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCCGCGTTCCCACGGTCCTTGCT‑GCGCGGGGACATCAATACGGGTGACGTCCCCCA‑GCTGTCCCTTAATACCCGGGTCAT‑CACCCCGACTGATATTCCGACTGCGG‑ATTGGGGAAACCA‑CCCTTATTTCC‑GTTGGCAGCAAGCGCACCGGAAGACTGTGTGTGGACCA‑CGGCAGAAGGGCGGGTgggg                                                                                                                       >  2:264179/1‑250 (MQ=255)
                                                                                                                                        ttACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGGCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTTGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGgcg                                                                                                                      >  1:370469/1‑251 (MQ=255)
                                                                                                                                           cGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCCGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTTA‑GTTGGCAGCAAGCGCATCGGAAGACAGTCTGTG‑ACCCGCGGCAGAAGGCCGGGTGGgcgc                                                                                                                     >  1:139205/1‑249 (MQ=255)
                                                                                                                                                 tGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCCCGCGTACCCTCCCGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCCCA‑GCTGTCCCTTAATACCCGGGGCCT‑CACCCAGAATGATATCGCCACAGCCA‑ATTGGCGAA‑CCAGCCCTGATTTCA‑GTTTGCAGCAAGCGCATCGGAAGACCGTCTGTG‑GCCCGCGGCAGAAGGGCGGTTGGGCGCcgcggttg                                                                                                             >  2:199022/1‑247 (MQ=255)
                                                                                                                                                    aaCATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGGACCGTCCAGCGTTCCCACGTTCCCTGCT‑GCGCGGGGGCATCAATACGGGTGACGGCCCCCA‑GCTGTCCCTTAATACCCGGGGCAT‑CACCCCGACCGATATCGCGACCGCGA‑ATTGGCGAA‑CCAGCCATGGTTTCA‑GTTGGCAGCAAGCGCATCGGGAGACCGGCCGTG‑ACCAGCGGCAGGAGGCCGGTTGGGCGCAGCGGTAGCAc                                                                                                          >  2:88945/1‑251 (MQ=255)
                                                                                                                                                     aCATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCAcc                                                                                                         <  1:316334/251‑1 (MQ=255)
                                                                                                                                                     aCATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCAcc                                                                                                         <  1:331841/251‑1 (MQ=255)
                                                                                                                                                      cATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCa                                                                                                        <  1:410734/251‑1 (MQ=255)
                                                                                                                                                                               tGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGAt                                                                                <  1:368171/249‑1 (MQ=255)
                                                                                                                                                                                gCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGc                                                                             <  1:269257/251‑1 (MQ=255)
                                                                                                                                                                                  ggCTCGTCCGCGCGTCCCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACC‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGCAGAAACCAGAATACTGGATGGCTc                                                                           >  1:110963/1‑251 (MQ=255)
                                                                                                                                                                                          cgcgcgTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CGCCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGc                                                                   <  1:328596/251‑1 (MQ=255)
                                                                                                                                                                                            cgcgTACCGTCCAGCGTTCCCACGTGCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCa                                                                 <  2:370469/251‑1 (MQ=255)
                                                                                                                                                                                                  ccGTCCAGCGTTCCCTCGTTCCTTTCT‑GCGCGGTGACATCAATTCGGTTGACGTCCCACA‑GCTGTCCATTAAGACCCGGGTCAG‑CACCCAGACTGATAGCGCGACATCTA‑ATTGGCGAA‑CCAGCCATTATTTGA‑TTTGGCAGCAAGCGCATCGGAGTAGGGTCTGTG‑GCCAGCGGCAGAATGCGGGTTCGGCGCAGCGTCAGCACCATTGGT‑GAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAAt                                                            <  1:430998/250‑1 (MQ=255)
                                                                                                                                                                                                      ccAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGa                                                        <  1:199022/250‑1 (MQ=255)
                                                                                                                                                                                                            ttCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATcaaca                                                  >  1:36946/1‑250 (MQ=255)
                                                                                                                                                                                                             tCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCGGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCGGAATACTGGATGGCTCAAAGAAGCCACGAATACGATcaacaa                                                 <  1:200942/250‑1 (MQ=255)
                                                                                                                                                                                                             ccccACGTTCCTTGCA‑GCGCGGTGACATCAATACGGTGGACGTCCCACA‑GCTGGCCAGTAATCCCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAAGACGATCAACAAt                                                <  2:468151/250‑1 (MQ=255)
                                                                                                                                                                                                                cACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCt                                             <  1:98678/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTTCCGTGCT‑GCGCGGTGAAATCAATACGTTTGACGTCCCACG‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATACCGCGACGGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTcgc                                          <  2:142268/251‑1 (MQ=255)
                                                                                                                                                                                                                     gccTTGCT‑GCGCGGTGACATCAATACGGTGGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAAGACTGGATGTGTCAAAGAAGCCACGAATACGAGCAACAATGCTcgcc                                         <  2:224289/249‑1 (MQ=255)
                                                                                                                                                                                                                      ccTTGCT‑GCGCGGTGACAGCAATACGGTTGACGCCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTcgccg                                        <  2:139205/250‑1 (MQ=255)
                                                                                                                                                                                                                         tgct‑gcGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACGGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCAGCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCtt                                     <  2:231065/250‑1 (MQ=255)
                                                                                                                                                                                                                            t‑gcGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGcc                                 >  2:195021/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gggACATCAATCCGGTTAACGCCCCACA‑GCTGTCCAGTAATACCCGGGTCAT‑CACCCAGACTGATATCGCTACTGCGA‑ATTGGCGAA‑CCTGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATTCGATCAACAATGCTCGCCGCTTCGCCACAGGt                           <  1:86438/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cATCAATCCGGTTGACGCCCCCCA‑GCTGGCCATTAAGACCCGGGTCCT‑CCCCCAGACTGATATCGCGACTCCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAAGACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATc                       <  2:167342/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATc                       <  1:337845/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ggTTGACGTCCCACA‑GCGGTCCATTAATCCCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATCGACACTTCa              <  2:36946/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATCGACACTTCAt             >  1:183318/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATCGACACTTCAt             >  1:320343/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        ccccACG‑GCTGTCCATTAATCCCCGGGTCCT‑CACCCAGACTGATATCGCGACGGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCGTCGCCACAGGTGATCGACACTTCATCGCCGtt      <  2:31841/250‑1 (MQ=255)
                                                                                                                                                                                                                                                            aca‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGAGCGACACTTCATCGCCGTTGCTc  <  1:336626/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                 
CGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGTCCGTCAACCCAAACGGTATCTTTCCAGCCTCTGTCTTCCGGGAACGGCATCGCGCCATTCACGTTACGGATCTGGAACATTACGCCTTCAATATGGAACGCTTGCGGCTCGTCCGCGCGTACCGTCCAGCGTTCCCACGTTCCTTGCT‑GCGCGGTGACATCAATACGGTTGACGTCCCACA‑GCTGTCCATTAATACCCGGGTCAT‑CACCCAGACTGATATCGCGACTGCGA‑ATTGGCGAA‑CCAGCCATGATTTCA‑GTTGGCAGCAAGCGCATCGGAAGACTGTCTGTG‑ACCAGCGGCAGAAGGCCGGTTGGGCGCAGCGTTAGCACCAG‑GGTAGAAACCAGAATACTGGATGGCTCAAAGAAGCCACGAATACGATCAACAATGCTCGCCGCTTCGCCACAGGTGATCGACACTTCATCGCCGTTGCTC  >  NC_000913/3161314‑3161815

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