Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,585 T→G 5.8% E199A (GAA→GCA)  malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,5850TG 5.8% 788.2 / 32.4 225E199A (GAA→GCA) malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base T (101/110);  new base G (5/8);  total (107/118)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 228 total aligned reads displayed.

CGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAGTGAAAACCGCTGTGCCAGC  >  NC_000913/1698316‑1698857
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
cGGTCAATGCCGCCCGCCGAACCTGTTGCGCAAAGTAACGATCTGTTCGGCTTTCCCCGCTTTGCCGCCCTGCTTTGGGCAACCCAAACCATCCCCCCATCGCATGAGTTTCTTTATGTCACACCACGTCACTGATGGTCGGTTTATGACTTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCa                                                                                                                                                                                                                                                                             <  4:106087/275‑1 (MQ=255)
      agcccgcaAGCGAAACCAGTCTCTCAACGTCACGCTGGACTCCCCTTTCCCCGCGTTGGCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGAAATCGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGt                                                                                                                                                                                                                                                                        <  4:275961/271‑1 (MQ=255)
                                                           ttttGCCTCCCCGATTTGAACAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTTATAGCTTCCTCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                                                                                                                                                                                                  <  1:38504/274‑1 (MQ=255)
                                                            tttGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCCGGCCAACTTTTAGTAAAGTTGCACAATAGCCCCCAACCCGTTCt                                                                                                                                                                                                                                                                                >  4:82653/1‑212 (MQ=255)
                                                                tcgccCTGCTTTGGGCAAACCGAACCAGGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGTACAATAGCCCCCAACCCTTTCTGCAc                                                                                                                                                                                                                                                                            <  2:82653/211‑1 (MQ=255)
                                                                                            gCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACTGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCa                                                                                                                                                                                  >  2:59582/1‑274 (MQ=255)
                                                                                            gCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTGTAGTAGAGTTTCACAATAGACACCAACCCGTCCTGCACGGGTTAATGAGGAACTTTGCCCTCCAAGCCAGGCGAGCCGCTGACGCCCATTGCGAATGAAAGGCCCCGTCAACAACTGTGCAGCCTTCa                                                                                                                                                                                  >  4:29725/1‑274 (MQ=255)
                                                                                                  aTCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTg                                                                                                                                                                                                                                                         >  2:250316/1‑197 (MQ=255)
                                                                                                                          aCCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTc                                                                                                                                                                                                                                                                                 <  8:748077/149‑1 (MQ=255)
                                                                                                                          aCCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTc                                                                                                                                                                                                                                                                                 <  7:1368072/149‑1 (MQ=255)
                                                                                                                            cACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTg                                                                                                                                                                                                                                                                               >  8:374501/1‑149 (MQ=255)
                                                                                                                                    tGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGtt                                                                                                                                                                                                                                                                       >  8:416580/1‑149 (MQ=255)
                                                                                                                                       tGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                                                                                                                                                                    <  8:1237352/149‑1 (MQ=255)
                                                                                                                                         gTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGa                                                                                                                                                                                                                                                                  <  7:1043996/149‑1 (MQ=255)
                                                                                                                                             ggTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGaa                                                                                                                                                                                                                                                              >  5:641061/1‑149 (MQ=255)
                                                                                                                                              gTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                                                                                                                                                                                                                             <  7:416580/149‑1 (MQ=255)
                                                                                                                                              gTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                                                                                                                                                                                                                             >  7:913096/1‑149 (MQ=255)
                                                                                                                                                tATGACGTAAAAGCGCCGTGATAGCTTGCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAAAACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAATCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                                                                                                                                                                                                       >  1:71755/1‑153 (MQ=255)
                                                                                                                                                   gACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                                                                                                                                                        <  8:1222608/149‑1 (MQ=255)
                                                                                                                                                   gACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGc                                                                                                                                                                                                                                                        <  6:14996/149‑1 (MQ=255)
                                                                                                                                                               ccGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                                                                                                                                                            <  7:305573/149‑1 (MQ=255)
                                                                                                                                                               ccGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                                                                                                                                                            >  5:348663/1‑149 (MQ=255)
                                                                                                                                                                   gATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGa                                                                                                                                                                                                                                        <  6:1188135/149‑1 (MQ=255)
                                                                                                                                                                    aTAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGACCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           <  2:310769/193‑1 (MQ=255)
                                                                                                                                                                     tAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATc                                                                                                                                                                                                                                      >  7:1006846/1‑149 (MQ=255)
                                                                                                                                                                        cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                                                                                                                                                   <  5:1093586/149‑1 (MQ=255)
                                                                                                                                                                        cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                                                                                                                                                   >  8:533703/1‑149 (MQ=255)
                                                                                                                                                                        cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                                                                                                                                                   >  6:382446/1‑149 (MQ=255)
                                                                                                                                                                        cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                                                                                                                                                                   >  6:480213/1‑149 (MQ=255)
                                                                                                                                                                         ttCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCt                                                                                                                                                                                                                                  >  5:106047/1‑149 (MQ=255)
                                                                                                                                                                             gcgGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATGTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATg                                                                                                                                                                                                                              >  8:242476/1‑149 (MQ=255)
                                                                                                                                                                                gcttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGccc                                                                                                                                                                                                                           >  5:1391655/1‑149 (MQ=255)
                                                                                                                                                                                gcttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCt                                                                                                                                                                                                          <  4:186277/166‑1 (MQ=255)
                                                                                                                                                                                 cttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCa                                                                                                                                                                                                                          <  7:276918/149‑1 (MQ=255)
                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAGCACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                                                                                                                                                                                                     >  5:900021/1‑149 (MQ=255)
                                                                                                                                                                                      ttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCg                                                                                                                                                                                                                     >  6:837208/1‑149 (MQ=255)
                                                                                                                                                                                          ggctggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATg                                                                                                                                                                                                                 >  6:308458/1‑149 (MQ=255)
                                                                                                                                                                                            ctggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgag                                                                                                                                                                                                               >  5:760994/1‑149 (MQ=255)
                                                                                                                                                                                             tggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgaga                                                                                                                                                                                                              >  7:375841/1‑149 (MQ=255)
                                                                                                                                                                                                 aGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCt                                                                                                                                                                                                          >  5:197296/1‑149 (MQ=255)
                                                                                                                                                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     >  5:307355/1‑149 (MQ=255)
                                                                                                                                                                                                      aCTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     >  8:1313672/1‑149 (MQ=255)
                                                                                                                                                                                                       cTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTca                                                                                                                                                                                                    >  5:578473/1‑149 (MQ=255)
                                                                                                                                                                                                        tCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaa                                                                                                                                                                                                   >  8:627319/1‑149 (MQ=255)
                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                                                                                                                                                                                                                          >  7:1147262/1‑92 (MQ=255)
                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACttt                                                                                                                                                                                                                                                          <  8:1147262/92‑1 (MQ=255)
                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTc                                                                                                                                                                                                         >  4:364566/1‑141 (MQ=255)
                                                                                                                                                                                                          cAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaaca                                                                                                                                                                                                 >  8:22768/1‑149 (MQ=255)
                                                                                                                                                                                                           aaCACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaacaa                                                                                                                                                                                                <  7:427641/149‑1 (MQ=255)
                                                                                                                                                                                                               cccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgt                                                                                                                                                                                            <  7:714203/149‑1 (MQ=255)
                                                                                                                                                                                                                ccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           <  6:1077158/149‑1 (MQ=255)
                                                                                                                                                                                                                ccAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                                                                                                                                                           >  8:182138/1‑149 (MQ=255)
                                                                                                                                                                                                                  aaTCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCa                                                                                                                                                                                         <  5:508418/149‑1 (MQ=255)
                                                                                                                                                                                                                  aaTCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTAGCCGGTCATGGCTACTTACTGCAGCTCCg                                                           >  4:145738/1‑275 (MQ=255)
                                                                                                                                                                                                                        tATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGc                                                                                                                                                                                   >  8:965996/1‑149 (MQ=255)
                                                                                                                                                                                                                          tGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 <  7:533703/149‑1 (MQ=255)
                                                                                                                                                                                                                            aaaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgt                                                                                                                                                                               >  6:1318137/1‑149 (MQ=255)
                                                                                                                                                                                                                             aaCGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAtgtt                                                                                                                                                                              <  6:900021/149‑1 (MQ=255)
                                                                                                                                                                                                                                ggCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTc                                                                                                                                                                           >  5:1203544/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCgg                                                                                                                                                                        <  7:159053/149‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCggg                                                                                                                                                                       >  8:937424/1‑149 (MQ=255)
                                                                                                                                                                                                                                     gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCa                                                                                                                                                                                  <  5:170663/137‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCa                                                                                                                                                                                  >  6:170663/1‑137 (MQ=255)
                                                                                                                                                                                                                                      ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga                                                                                                                                                                                                                <  7:213150/106‑1 (MQ=255)
                                                                                                                                                                                                                                      ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCg                                                                                                                                                                     >  7:698875/1‑149 (MQ=255)
                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                                                                                                                                                                 <  5:36781/136‑1 (MQ=255)
                                                                                                                                                                                                                                       cAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGa                                                                                                                                                                    >  5:862935/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aaaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATACGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGaa                                                                                                                                                                   <  8:101789/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATcc                                                                                                                                                                                                                                     <  8:1078600/82‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCg                                                                                                                                                                >  5:1040285/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGt                                                                                                                                                               <  7:1224498/149‑1 (MQ=255)
                                                                                                                                                                                                                                             ttAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTa                                                                                                                                                              <  5:1318137/149‑1 (MQ=255)
                                                                                                                                                                                                                                               agtagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                                                                                                                                                            <  6:1040285/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCaa                                                                                                                                                          <  6:641061/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCaa                                                                                                                                                          <  6:760994/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCaa                                                                                                                                                          <  7:1179112/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTGAATGAGGAACTTTGCCCTCCCCGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTa                         <  2:293810/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATg                                                                                                                                                                                                            >  6:41362/1‑94 (MQ=255)
                                                                                                                                                                                                                                                      ttGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATg                                                                                                                                                                                                            <  5:41362/94‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc                                                                                                                                                   <  6:348663/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGACCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACAtcatca                                                                                                                                                  >  6:1256765/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTc                                                                                                                                                                                                     <  2:45454/97‑1 (MQ=255)
                                                                                                                                                                                                                                                          acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCaa                                                                                                                                                 <  7:1313672/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGa                                                                                                                                               >  8:1181063/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCt                                                              <  2:169144/230‑1 (MQ=255)
                                                                                                                                                                                                                                                             aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGAt                                                                                                                                              <  8:499656/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGAt                                                                                                                                              <  6:1203544/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              tAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATACCTGACACGGGAAGCGCAAATCACCGGCATAGCTTCTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGCCAATTCCCACACCATCGACCCCCTGATTCAGTCACAGTGa                  >  4:69696/1‑272 (MQ=255)
                                                                                                                                                                                                                                                               aGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATaa                                                                                                                                            <  7:987487/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATaa                                                                                                                                            <  7:1108934/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  >  6:143116/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTcaac                                                                                                                                                                                                  <  5:143116/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc                                                                                                                                                   <  6:127932/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc                                                                                                                                                   >  5:127932/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                    ttaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGc                                                                                                                                       >  8:810022/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                        cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGc                                                                                                                                                                                          <  8:1291942/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGc                                                                                                                                                                                          >  7:1291942/1‑94 (MQ=255)
                                                                                                                                                                                                                                                                         ccGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACggg                                                                                                                                  >  5:1134776/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGCTCCGCTGATGTCCATCGCGACTGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCCACATCATCAAGATAACTTGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCCAAGGTCATCGTTACGTCCTGCAGCTCCGGCAATTACCCCCCCATCCACACCCTGATTCAGTAACAGTGCAAACCGCCTTGCCCgc  >  4:59848/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            ttCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAg                                                                                                                               <  5:849147/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
CGGTAAATGCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAGTGAAAACCGCTGTGCCAGC  >  NC_000913/1698316‑1698857

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: