Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,251,973 T→C 15.9% L303P (CTC→CCC)  dhaR → dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,251,9730TC15.9% 188.9 / 4.4 75L303P (CTC→CCC) dhaRdhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base T (48/14);  new base C (3/10);  total (51/24)
Fisher's exact test for biased strand distribution p-value = 3.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCG  >  NC_000913/1251741‑1252206
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
cAATTTATAAAGTCCCAGGGGGCGGGGGTTTGGCCCCTTTGCCGCAGGGCAAGTAGGGGGAGGGAAAACCACAGAACTTTAAATTTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATc                                                                                                                                                                                                                                      <  1:1261042/239‑1 (MQ=255)
 aaTTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCAt                                                                                                                                                                                                                                       >  1:1020593/1‑236 (MQ=255)
          aTGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTTTTTCCCCATCCTGGGGaaaaa                                                                                                                                                                                                                           >  2:566428/1‑236 (MQ=255)
               cAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCCCCATCCTGTGGaccgaatgc                                                                                                                                                                                                                       >  1:1061158/1‑230 (MQ=255)
                 ggcggcGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGc                                                                                                                                                                                                                    >  1:1325656/1‑239 (MQ=255)
                                      ttGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAAttt                                                                                                                                                                                                >  1:1069069/1‑238 (MQ=255)
                                       tGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAAATGCGGCCGTTGATGACCAGTCAATTAg                                                                                                                                                                                               >  2:58740/1‑239 (MQ=255)
                                                                   aacccgtGAACTTTTAACTTTCCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCata                                                                                                                                                                    <  1:566428/233‑1 (MQ=255)
                                                                           aaCTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGGAAAGTCAGCCATACCTTCGCTCATATGCCACAg                                                                                                                                                            >  1:766706/1‑239 (MQ=255)
                                                                           aaCTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATAc                                                                                                                                                                               >  1:500856/1‑220 (MQ=255)
                                                                                ttAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACg                                                                                                                                                        >  1:382257/1‑238 (MQ=255)
                                                                                ttAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAATAAACCAAACCTTCGCCCAAAAGCCAAAGGagg                                                                                                                                                        >  2:536964/1‑236 (MQ=255)
                                                                                  aaCGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCTTATGCACCAGGACGAt                                                                                                                                                      >  1:175526/1‑238 (MQ=255)
                                                                                       tACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGCAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCaa                                                                                                                                                >  1:414858/1‑239 (MQ=255)
                                                                                        aCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATGATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGAGCCGGCAGTTGTTGACCAGTCAATTAGGAAAAGTAAGCCATACCTTCCCCCATATGCACCAGGAAAAACCGCaa                                                                                                                                                >  2:37324/1‑238 (MQ=255)
                                                                                        aCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCaa                                                                                                                                                >  1:32458/1‑238 (MQ=255)
                                                                                        aCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCaa                                                                                                                                                >  1:770110/1‑238 (MQ=255)
                                                                                        aCCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTAATGCGGTGATAACCCTTAAACCGATAATAGAAACGACGGGAACCAGCTTTATTTTGTTGCTCCATCCTGGGGAACAGATGCGGCAGTGGCAGACCAGTCACTTCGGCAAAGACAGCCAAACCTTCGCTCATTTGCCAAAGGCCGAACCGCaa                                                                                                                                                >  2:87787/1‑238 (MQ=255)
                                                                                              cgTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGGAAACccgc                                                                                                                                          >  1:264375/1‑238 (MQ=255)
                                                                                                    ggaacaaGCAAAAAAACCGGCAAATCCGGTCAAACACGTAAAACCACCCTGTAAAAGCCAGCACCAGTTTATCGATCGGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCCCAATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGa                                                                                                                                   <  1:762216/237‑1 (MQ=255)
                                                                                                    gCAAAAAGAAATAAAACAGGCACACCCGCTAAAACACGAAAAACCAACCTGTAAAAGCAAGCACCAGTTAATGGAGGGGTTGATAACCTTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAATTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGa                                                                                                                                   <  2:1020593/239‑1 (MQ=255)
                                                                                                        caaGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCCCTTGATTCa                                                                                                                               >  1:1162494/1‑239 (MQ=255)
                                                                                                          aGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCTGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCAtt                                                                                                                             >  1:615370/1‑239 (MQ=255)
                                                                                                          aGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATCCGGCAGTTGATGACCAGTCAATTAGGAAAAGACAGCCAACCCTTGGTGCATATGCCCAAGGCCGATCCGACAACCCGCCGATTGTTTCttt                                                                                                                             >  1:1318887/1‑239 (MQ=255)
                                                                                                           gCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGGAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCtttt                                                                                                                            >  1:807561/1‑239 (MQ=255)
                                                                                                               tAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTGGCTCAAATGCCCCATGACAGACCCCAAACCCCCCCCTGGAATCATTTTGGt                                                                                                                        >  2:1051268/1‑239 (MQ=255)
                                                                                                                      ggCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAg                                                                                                                  >  1:680893/1‑238 (MQ=255)
                                                                                                                      ggCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATCCGGCCGTTGATGACCAGTCCATTAGGACAAGTCAGCCCTACCTTCGCCCATATGCCACAGAGCGATGCCCAAACCCCCCCCTTGATTCATTTTgaccgacggg                                                                                                                 >  2:715199/1‑230 (MQ=255)
                                                                                                                       gCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAggc                                                                                                                >  1:1127884/1‑239 (MQ=255)
                                                                                                                           aTCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCAAATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTctgcctcttata                                                                                                            >  1:78294/1‑228 (MQ=255)
                                                                                                                              cgttCAAACCCGTAGAAGCAACTTTTGAAAGCCAACACCAGTTTATTGATGCGGTGATAACCCTTAAACCGAAAATAGAAACGCAGGGAACCAGTTTTATTTTGTTGCTCCATCCTGTGGAACAGAGGCGGCAGTTGAGGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGgcgcgc                                                                                                         <  1:1296017/236‑1 (MQ=255)
                                                                                                                                   aaaCCCGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTTATACCCCTTAAACCGATAATAGAAACGCAGGGAACCAGTTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                                                                    <  1:916654/239‑1 (MQ=255)
                                                                                                                                   aaaCACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATACCCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAg                                                                                                    <  1:536964/239‑1 (MQ=255)
                                                                                                                                   aaaCACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAAACCGCCGCTTGATTCATTTTGGCCGCCAGGCGgcgcgcggaa                                                                                                     >  2:551658/1‑234 (MQ=255)
                                                                                                                                     acacGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGCTGCGGCAGCTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTGCCCCATATGCCACAGGCCGATCCGCAACCCCCCCCCTTGATTCTTTTTTGTCCCCAGGCTGCCCGCAgagcct                                                                                                  >  2:82456/1‑234 (MQ=255)
                                                                                                                                         gTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCa                                                                                                                                                                                                        >  2:1306127/1‑133 (MQ=255)
                                                                                                                                         gTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCa                                                                                                                                                                                                        <  1:1306127/133‑1 (MQ=255)
                                                                                                                                                 caCCTTTTAAAACCAGCCCCCAGTTTATTGATGCGGTGAAACCCCATAACCAGATAAAAGAAACGCAGGGAACCCGCTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATAAGGAAAAGGCAGCAATTCCTCCGCTCATGGGCCACAGGCCGATCCGAAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTtcctgtcctg                                                                                       <  2:32458/238‑1 (MQ=255)
                                                                                                                                                  cccTTTAAAAGCCAGCACCAGTTTATGGATGCGGTGATAACCCTTAACCCGATAATAGAAACGCAGGGAACCAGTTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                                                                     <  1:58740/238‑1 (MQ=255)
                                                                                                                                                  aCCTTTAAAACCCAGCACCAGTTTATTGATGCGTGTAAAACCCTTAAACCAATAATAGAACCGCGGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGAGGACCAGTCAATTAGGAAAAGTCAGCCATCCCTTCGCTCATATGCCGCAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                                                                     <  1:744113/239‑1 (MQ=255)
                                                                                                                                                       tGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGTAAAGTAAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCAGCCGCTTGATTCATTTTGGTCCCACGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTtttg                                                                                >  1:1331581/1‑239 (MQ=255)
                                                                                                                                                        gAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGTTCCATCCTGTGGAACAGATGCGGCAGTTGAGGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                                                               <  1:1302335/239‑1 (MQ=255)
                                                                                                                                                        gAAAACCGGCCCCAGTTTTTTGATTGGGTTATCACCCTAAAACCTATAATAAAAACGCAGGGAACCGGTTTTATTTTGTTGCTCACTCCTGTGAAACAGAGGCGGCGGTTGATGACCAGTCAATTAGGAAAAGTCAGCCAATCCTTCGCTCATATGCCCAAGGACGATCCGAAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTgg                                                                               <  2:1331581/239‑1 (MQ=255)
                                                                                                                                                              cAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCACCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGggaaga                                                                         >  1:141040/1‑239 (MQ=255)
                                                                                                                                                                    cccattttttGGTGCGTTGTATACCACTAAAACCGATAAGAGAAACGAAGAGAACCAGCTTTATTTTTTTGCTCCACCCTGTGAAACAGGCGCGGCAGTGGATGGCCAGTCATTTAGGAAAAGCCAGCCATACCTGCGCTCATATGCCACAGGACGTCCCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGc                                                                     <  2:1061158/235‑1 (MQ=255)
                                                                                                                                                                       gcttATTGAGGCACTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTCATGCTCCATCCGGTGGAACAGAGCCGGCAGTTGATGATCATTCAGTTAGGAAAAGGCAGCCCTACCTTGGACCATATGCCCCAGGACGCTCCGCAAACCCGCCACTTGACTCCTTTTGTACGCCAGTCGCCCAGCGTAACCATTCCTGTCCACCTTTGTGGGAAAGAgggggggg                                                                 >  2:1333924/3‑234 (MQ=255)
                                                                                                                                                                             tgagtgggTGATAACCCTTAACCCGATAATAGAAACGCAGGAAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACATATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAgg                                                           <  1:1012700/233‑1 (MQ=255)
                                                                                                                                                                               aTGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGc                                                          >  1:669065/1‑238 (MQ=255)
                                                                                                                                                                                 gCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGGGAAGAGGGCGTGGGAAAGGcccg                                                       >  1:676858/1‑236 (MQ=255)
                                                                                                                                                                                  cGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACt                                                       >  1:314783/1‑238 (MQ=255)
                                                                                                                                                                                  cGGCGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACCGGACGATCCGCAAACCCGCCGCCTGATTCATTTTGGTCCCCAGGCGGGCCGCCGTAGCGTTCCGTACCTACTTTGTGGAGAAAAGGGGGTGGGCCAGGAAcg                                                       >  2:467904/1‑237 (MQ=255)
                                                                                                                                                                                   ggTGAAAACCCTTAAACCGATAATAGAACCGCAGGGAACCAGCTTTATTTTGTGGCCCCATCCTGTGGACCAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTg                                                      <  1:754224/238‑1 (MQ=255)
                                                                                                                                                                                    gTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                                                            >  2:634043/1‑183 (MQ=255)
                                                                                                                                                                                     gttaTAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGgcggcg                                                                                                            <  1:634043/180‑1 (MQ=255)
                                                                                                                                                                                       ttattcccTTTAACCCAAAAAAGAAACCCGAGGAACCCATTTTTTTTTTGTTGCCCCATCCTGTGGAACAGAGTCGGCAGTTGATGCCCAGACAATTAGGAAAAGTCAGCCATCCCTTCGCTCATAGCCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCt                                                    <  2:770110/233‑1 (MQ=255)
                                                                                                                                                                                             ctttAACCCGTTAATAGAAACGGGGGGAACCATTTTTTTTTTTTTGCCCCATCCTGGTGACAGGGTGGGGGAGTTGATCACCAGTCAATTAGGAAAAGTCAGCCATCCCTTCGCTCATATGCCACGGGACGATCCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAgg                                            <  2:676858/237‑1 (MQ=255)
                                                                                                                                                                                             ccTTAACCCGATAATAGAAACGCAGGGAACCGGCTTTATTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGAGGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAgg                                            <  1:37324/239‑1 (MQ=255)
                                                                                                                                                                                                tAAACCGATAATAGAAACGCAGGGAACCAGTTTTATTTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                                                          <  1:951332/207‑1 (MQ=255)
                                                                                                                                                                                                tAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGagaag                                                                          >  2:951332/1‑207 (MQ=255)
                                                                                                                                                                                                    aagagaataGAACCGGAGGTATCCAGTTTTATTTTGTTGCCCTAGCGTGTGGAACAGATGCGGCATTTGAGCACCAGTAAATTGGGAAAGGTCACCCATCCCTTCGCTCATATGCCACGGGACGTTCCGCAACCCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCa                                     <  2:255543/234‑1 (MQ=255)
                                                                                                                                                                                                       ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGTTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTTGGGGAAAAAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCaaaa                                  >  2:977342/1‑239 (MQ=255)
                                                                                                                                                                                                       ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTTTGGAGAAGAGGGCGTGGGCAAGGCCCTGCTAAATCAGGCCATTCATaa                                  >  2:1092112/1‑239 (MQ=255)
                                                                                                                                                                                                       ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCCCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGGTGGCTGGGGCAAGGCACTGCTAAGTCAGGCAATTCCTaa                                  >  1:1210525/1‑239 (MQ=255)
                                                                                                                                                                                                       ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCCGCAGGCGCCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATaa                                  >  1:535483/1‑239 (MQ=255)
                                                                                                                                                                                                       ataataGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCAGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACCGGACGATCCGCAAACCCGCCGCTTGATTCATTTTTGTCGCCAGGTGGCGCGCGGTACCTTTCCTGTCCTGCTTTGTGGGGAAAGGGGCGGTGGCAGGGCACTGCTAAGTACGGCAATTctaaa                                  >  2:1074851/1‑235 (MQ=255)
                                                                                                                                                                                                                 cGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGCTTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGGGGGCAAGGCACTGCTAAGTCAGGCACTTAATAATGAAagcgag                        >  2:264658/1‑239 (MQ=255)
                                                                                                                                                                                                                 cGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAagg                           >  1:255543/1‑235 (MQ=255)
                                                                                                                                                                                                                      gggAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGGGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAaacgaggggg                    >  2:1085420/1‑235 (MQ=255)
                                                                                                                                                                                                                       ggAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCCCCGCTTCATTCATTTTTGTCGCCACGCGGGGCCCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGGCCCTCTAAGGCAGGCAATTCTTAATGAAAACGCGCGggc                   >  2:1332852/1‑239 (MQ=255)
                                                                                                                                                                                                                         aaCCAGTTTTAGTTTGTTGCTCCAACCGGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACg                                                                                                                                                        <  2:648262/104‑1 (MQ=255)
                                                                                                                                                                                                                          aCCGGCTTTTTTTTGTTGCCCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCATGTAGGAAAAGTCAGCCATCCCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGtgctgc                <  1:1051268/239‑1 (MQ=255)
                                                                                                                                                                                                                               cTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCt                                                                                     >  1:1011815/1‑165 (MQ=255)
                                                                                                                                                                                                                                   aTTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAAc                                                                                                                                              >  1:648262/1‑104 (MQ=255)
                                                                                                                                                                                                                                          ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCCAACCCGGCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGGCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGAGTGCTGCAGGTCCTTATATCg  >  1:305741/1‑238 (MQ=255)
                                                                                                                                                                                                                                          ttgCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCCAGGCCCTGCTAAGCCAGGCAATTCATAATGAAAACGAGCGGGCGGCAGGTCCTTATATcc  >  2:661254/1‑237 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
CAATTTATTAATGCCCAGGCGGCGCGGGTCTTGCGCCTTGACGCGACGGCAAGTCAGGGACGGGCAATCACTGAACTCTTAACGTTACCCGCCGTATTGCAACAAGCAATAAAACAGGCACATCCGCTCAAACACGTAGAAGCAACCTTTGAAAGCCAGCACCAGTTTATTGATGCGGTGATAACCCTTAAACCGATAATAGAAACGCAGGGAACCAGCTTTATTTTGTTGCTCCATCCTGTGGAACAGATGCGGCAGTTGATGACCAGTCAATTAGGAAAAGTCAGCCATACCTTCGCTCATATGCCACAGGACGATCCGCAAACCCGCCGCTTGATTCATTTTGGTCGCCAGGCGGCGCGCAGTAGCTTTCCTGTCCTGCTTTGTGGAGAAGAGGGCGTGGGCAAGGCACTGCTAAGTCAGGCAATTCATAATGAAAGCGAGCGTGCTGCAGGTCCTTATATCG  >  NC_000913/1251741‑1252206

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