Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913726,3440GA28.3% 83.3 / 14.7 46V231V (GTC→GTT) kdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (6/27);  new base A (13/0);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 2.67e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTCCCCGGTGATGGCGCTTTCA  >  NC_000913/726121‑726575
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
gTAGCAGAACGTAAACGTCCCCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCa                                                                                                                                                                                                                             <  1:923575/239‑1 (MQ=255)
 tAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGTCCGTAAAAAGAAGCAAGTCAGGGGaa                                                                                                                                                                                                                            >  1:484665/1‑239 (MQ=255)
  aGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAAt                                                                                                                                                                                                                           <  1:399414/239‑1 (MQ=255)
      gAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGATTGCGGGACAAAAAAGGATAGCAAGGGCagcagc                                                                                                                                                                                                                        >  2:399414/1‑234 (MQ=255)
         cgtcaacgtcaCCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGGAAAAAAACGATAGTCAGGGGAATCAGCAGa                                                                                                                                                                                                                    >  1:12302/1‑239 (MQ=255)
                                         ccGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCg                                                                                                                                                                                   <  1:541286/239‑1 (MQ=255)
                                                    aaTCACGTTGCCGCCTAGCATCCGTCGCAGCCCGGCGCCCCGGGTCGCTAACGACAGGCCCCCAATAGTCATTGGGATCAAACAGACCAGCAGCCCCACCAGTACCGTTACGCTGACTGCATTTCCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑GTGCCAGTAAAAAGACGATAGTCGGGGCAATCAGCAGGATGGTCAGGGC‑AAAGTCGTTCGGCGTTAGGCGTCGGGGTgcg                                                                                                                                                                         <  2:435543/238‑1 (MQ=255)
                                                       cacaTGCGCGCCTAGCGTCGGGCTCATCCCGCCGACGCCGATCGCTGACAACGTGCCGCCAAAGGTGGTGGGGATCAGACAGCCCAACAGCGCCACCAGTACCGTTACGCTGACTGCATTGCCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTgcg                                                                                                                                                                         <  2:331051/235‑1 (MQ=255)
                                                         caTTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAAACGACAGTCAGGGCAATCAGCAGAATGGTCCGGGC‑AAATTCGTTCGGCGTTTTGCGTCGCCGTGCGCCTTcc                                                                                                                                                                   >  1:347221/1‑239 (MQ=255)
                                                                       aTCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGGGTGGATCAGACAGCACAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATttcgttcg                                                                                                                                                                                             <  2:837509/199‑1 (MQ=255)
                                                                                 tCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCCCGCGGAAACCGGCCACAGCGTGGCGG‑TTGCCCGTAAAAAAACGAAAGTCAGGGCAATCAGCAGAATGGGCAGGGG‑AAATTCGTTCGGGGTTTTGCg                                                                                                                                                                                   >  1:837509/1‑199 (MQ=255)
                                                                                             cGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGcc                                                                                                                              <  1:234213/238‑1 (MQ=255)
                                                                                                    tGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTa                                                                                                                       <  1:786424/238‑1 (MQ=255)
                                                                                                    tGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTa                                                                                                                       <  1:684585/238‑1 (MQ=255)
                                                                                                     gACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAACCGCTAGTGAGGGCAAACCGCACAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTTTGCGCCTTCTACCAT‑CGCGATCAGACG‑ATCCCGACAAGTCTTGCCGGGGTTAAc                                                                                                                     >  2:241221/1‑239 (MQ=255)
                                                                                                        aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGcc                                                                                                                                                                       <  1:117435/188‑1 (MQ=255)
                                                                                                                gccAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGTCTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑GTGGCAGGAAAAAAACGACAAGCAGGGGAAAACACAGAAAGGTAGCGGC‑AAATTTGTTTGGGGGTTTTGGTGGCCGTTGGCCTTCCCCCCT‑CGCGATCATCCG‑CGCCAGAAAAGGTTCGCCCGGGTGGAAGCTACACTCa                                                                                                           >  2:860820/1‑238 (MQ=255)
                                                                                                                  cAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑GTGCCAGTAAAAAAACGAAAGTCAGGGCAAACAGCAGAATGGGCAGGGG‑AAATTCGTTCGGCGGTTTTCGTCGCTGTGCGCCTTCCACCAT‑c                                                                                                                                                            >  2:117435/1‑188 (MQ=255)
                                                                                                                        tGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCCGGAAAAAAACGAAAGGCAGGGGAATCAGCAGGAAGGGCAGGGCAAATTTGTTTCGGGGTTTT‑CGTCGCTGGGGGGCTTCCCCCCT‑CCCGGTCATCCGAATCCAGAAATGCCGCGCCGGGGTTAACGCTACACCCAAACAcca                                                                                                    >  2:134790/1‑238 (MQ=255)
                                                                                                                            tGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATttcgttcg                                                                                                                                                                                             >  1:905569/1‑146 (MQ=255)
                                                                                                                            tGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATttcgttcg                                                                                                                                                                                             <  2:905569/146‑1 (MQ=255)
                                                                                                                               gATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCtgtg                                                                                                                                                                           <  1:869802/161‑1 (MQ=255)
                                                                                                                                         acCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAATCGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAAACGATAGTCAGGGGAATCAGCAGAATGGGCAGGGG‑CATTTTGTTTGGGGGTTTGCGTTGCTGTTGGCCTTtcaa                                                                                                                                                                 >  2:869802/1‑160 (MQ=255)
                                                                                                                                          ccAGCAGCGCCACCAGTACCGTTACGCTGACTGGATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCaga                                                                                          <  1:646905/229‑1 (MQ=255)
                                                                                                                                              cagcGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCAccagc                                                                                                  <  1:230041/217‑1 (MQ=255)
                                                                                                                                               agcGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAAt                                                                                                                                                                                                                           <  1:416797/98‑1 (MQ=255)
                                                                                                                                                acgcCACCAGTACCGTTACGTGGACTGCATTACCGCCCCACGCGGAAACCGGCCACAGCGTGGCGG‑TTGACAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGGCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTgccg                                                                          <  2:525581/238‑1 (MQ=255)
                                                                                                                                                 cgcCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AAtt                                                                                                                                                                                                    <  2:228898/118‑1 (MQ=255)
                                                                                                                                                 cgcCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AAtt                                                                                                                                                                                                    >  1:228898/1‑118 (MQ=255)
                                                                                                                                                   ccaccaGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCa                                                                                                                                                                                                                                  <  2:731452/87‑1 (MQ=255)
                                                                                                                                                   ccaccaGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCa                                                                                                                                                                                                                                  >  1:731452/1‑87 (MQ=255)
                                                                                                                                                     accaGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGt                                                                                                                                    <  1:167049/176‑1 (MQ=255)
                                                                                                                                                         gTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCCGTAAAAAAACGAAAGTCAGGGGAAACAGCAGAATg                                                                                                                                                                                                                 >  2:416797/1‑98 (MQ=255)
                                                                                                                                                               ttACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAAAAAAAAAAAGAAAGACAGGGGAAACAGCAGAAAGGGCAGGGC‑AAATTCGTTCGGGGATTTTCGTTCCTTTGTGCCTTCCACCATACCCGA‑CAACCC‑ATCCCGAAATGTCTCGCCcggg                                                                                                                          >  2:167049/1‑176 (MQ=255)
                                                                                                                                                                    cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAAACGAAAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCC‑CTCCAGAAAAGTCTCGCCCGGGGTAACGCTACACCCCATCACCAGCCCGACAGAAAGAAAACGCGCGCCGCACGTGACGGAGGCAAAATc                                                      >  2:318286/1‑239 (MQ=255)
                                                                                                                                                                         tGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTgcc                                                                           <  1:7740/213‑1 (MQ=255)
                                                                                                                                                                            aTTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCaga                                                                                          <  1:266264/195‑1 (MQ=255)
                                                                                                                                                                              tACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTa                                                                                                                 <  1:299709/170‑1 (MQ=255)
                                                                                                                                                                                        cgcgGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAAACGATAGGCAGGGCAATCAGCAGAATGGGCAGGGC‑AATTTCGTTTGGCGTTTTTCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCCGAAATGTCTCGCCGGGGGTTACGCCACAatcaatca                                                                                                       >  2:299709/1‑170 (MQ=255)
                                                                                                                                                                                                     cacaGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTc                     <  1:855537/239‑1 (MQ=255)
                                                                                                                                                                                                          cGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCg                                                                                                                                                 <  1:845952/111‑1 (MQ=255)
                                                                                                                                                                                                           gTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTccc                   <  1:134790/235‑1 (MQ=255)
                                                                                                                                                                                                                    tGCCAGTAAAAAGCCGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACCCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTCCCCGGTGATGgcg        <  1:318286/238‑1 (MQ=255)
                                                                                                                                                                                                                    ggCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTCCCCGGTGATGgcg        <  2:347221/237‑1 (MQ=255)
                                                                                                                                                                                                                          tAAAAAGACGATAGTCAGGGCAATCAGCCGAATGGGCAGGGC‑AATTTCGTTCGGCGGTTTTCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGGTAAACCTACACTCCATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCCAAATCGCCCCCCGGTTCACGGAGCACCGGTGCCCATTCCCCGGGGATGGCGCATTCa  >  2:507463/1‑238 (MQ=255)
                                                                                                                                                                                                                          tAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTCCCCGGTGATGGCGCTTTCa  >  1:249750/1‑238 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
GTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG‑TTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGC‑AATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAT‑CGCGATCATCCG‑ATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAATCGCCGCCGGATTCACGGATCACCGGTGCCGATTCCCCGGTGATGGCGCTTTCA  >  NC_000913/726121‑726575

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