Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,099,855 C→T G224D (GGC→GAC)  rhaT ← rhamnose/lyxose:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,099,8550CT98.5% 207.6 / NA 66G224D (GGC→GAC) rhaTrhamnose/lyxose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/28);  minor base A (1/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGCCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCA  >  NC_000913/4099608‑4100098
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
gacgaTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGCGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATAATCGcgtcgc                                                                                                                                                                                                                                                  >  2:511551/1‑251 (MQ=255)
    aTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCgtcgccgcc                                                                                                                                                                                                                                              >  1:568587/1‑251 (MQ=255)
    aTACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCgtcgccgcc                                                                                                                                                                                                                                              <  1:161409/251‑1 (MQ=255)
    aTACCGCCGAACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCgtcgccgcc                                                                                                                                                                                                                                              >  2:144118/1‑251 (MQ=255)
           cgcACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATTCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATACAACAGAAACCGAGGTTAATGATCGcgtcgcc                                                                                                                                                                                                                                                 >  2:148171/1‑241 (MQ=255)
           cgcACAATACATAGAAACGCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACATAAACCGATGGTAATGATCGcgtcgcc                                                                                                                                                                                                                                                 <  1:148171/241‑1 (MQ=255)
                aaTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGacaac                                                                                                                                                                                                                                  >  1:427128/1‑251 (MQ=255)
                 aTACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGacaaca                                                                                                                                                                                                                                 <  1:307504/251‑1 (MQ=255)
                     aTAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGc                                                                                                                                                                                                                             <  2:164062/251‑1 (MQ=255)
                          aaCTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGc                                                                                                                                                                                                                        <  1:532621/251‑1 (MQ=255)
                            cTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGc                                                                                                                                                                                                                        >  2:479696/1‑249 (MQ=255)
                                 atGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCg                                                                                                                                                                                                                  <  1:420610/250‑1 (MQ=255)
                                   gcagcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGAc                                                                                                                                                                                                                >  1:439971/1‑250 (MQ=255)
                                    cagcaTCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACAt                                                                                                                                                                                                              >  1:379917/1‑251 (MQ=255)
                                             aCTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGAt                                                                                                                                                                                                     >  1:540424/1‑251 (MQ=255)
                                               tGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGAATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCg                                                                                                                                                                                                   >  2:74579/1‑251 (MQ=255)
                                                 aTGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGAc                                                                                                                                                                                                 >  1:18628/1‑251 (MQ=255)
                                                        cATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCAGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAg                                                                                                                                                                                           <  1:467133/250‑1 (MQ=255)
                                                            cTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCgg                                                                                                                                                                                      >  2:611219/1‑251 (MQ=255)
                                                             tGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGc                                                                                                                                                                                     >  1:283069/1‑251 (MQ=255)
                                                                   ggAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGc                                                                                                                                                                               >  1:442422/1‑251 (MQ=255)
                                                                        gCGGGCGTGTCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGGCATTCAGTGGATCGAGGCCAAGGGCGGCAGCGGCTTCa                                                                                                                                                                           <  2:71436/250‑1 (MQ=255)
                                                                           ggCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCa                                                                                                                                                                       <  2:84462/251‑1 (MQ=255)
                                                                             cGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATc                                                                                                                                                                     <  2:568587/251‑1 (MQ=255)
                                                                                               aaaaGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATc                                                                                                                                                   >  1:613190/1‑251 (MQ=255)
                                                                                                  aGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAAAACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGAGCGACGCCAAGTGCGGCAGCGGCTTCCTGCATCGGTTTTGCGGCGTTCATCGc                                                                                                                                                 >  1:5657/1‑250 (MQ=255)
                                                                                                    aaTTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCaa                                                                                                                                               >  1:184274/1‑250 (MQ=255)
                                                                                                     aTTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAg                                                                                                                                             <  1:229124/251‑1 (MQ=255)
                                                                                                        gCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACTAATAAAACAGAAACCGAGGTTAATGAGCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATCCAGTGGATCGAAGCCAAGGGCAGTAGCGGCTTCATGCATCGGTTTGGCGTCGTTCATCGCAAAGGa                                                                                                                                           <  2:392762/250‑1 (MQ=255)
                                                                                                              aCCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCg                                                                                                                                    <  2:18628/251‑1 (MQ=255)
                                                                                                                cacaTCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGc                                                                                                                                  >  1:155016/1‑251 (MQ=255)
                                                                                                                      aaCCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGTCAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAgaga                                                                                                                             >  2:23840/1‑250 (MQ=255)
                                                                                                                      aaCCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTgg                                                                                                                                                                                                                         >  1:346412/1‑158 (MQ=255)
                                                                                                                      aaCCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTgg                                                                                                                                                                                                                         <  2:346412/158‑1 (MQ=255)
                                                                                                                         ccGCCCAGTGTCGAGAGTAACCCATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAAt                                                                                                                         <  1:275037/251‑1 (MQ=255)
                                                                                                                           gCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCAGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGc                                                                                                                       >  2:244907/1‑251 (MQ=255)
                                                                                                                            cccAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGCCAAATCCTTCCCTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGcc                                                                                                                      >  1:80247/1‑251 (MQ=255)
                                                                                                                                  gtCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACAt                                                                                                                <  1:224204/251‑1 (MQ=255)
                                                                                                                                     gagaGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCAc                                                                                                             >  2:583518/1‑251 (MQ=255)
                                                                                                                                       gagTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGAGCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGGTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCg                                                                                                           >  1:324579/1‑251 (MQ=255)
                                                                                                                                          gaacgaatagtgAATTATCAGCGATTTTGCCGGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCCCTTTTGCCAGTCGAAGAAATCAGAAGCCGAGTTTAATGATCGCGTCGCCGCCCATGGTGACAACATAGCTTGGCGGAGCTACATACATTGGATCGCCGCCAATTGCAGCAGCGGCTTCTTGCATCGTTTTTGCTCCATTCTTCGCAAATTTCATCCCGGCTGAGAATATGCATCACATCACCGCCa                                                                                                        <  2:374661/242‑1 (MQ=255)
                                                                                                                                                       aTGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGTTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACACCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCttt                                                                                           <  1:75941/251‑1 (MQ=255)
                                                                                                                                                                 aTTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTg                                                                                 >  2:510796/1‑251 (MQ=255)
                                                                                                                                                                       ccAGCGAGAAGTCGGCTTTTAGCGACAACTCCTTCACTTTTGACAGACGAATAACACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGCCAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCACGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAAcact                                                                           >  1:442495/1‑248 (MQ=255)
                                                                                                                                                                         aGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGATTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCAAATCACCACCAGCACCAGCCCttttttgcgatggaacacttc                                                                         >  2:346905/1‑236 (MQ=255)
                                                                                                                                                                         aGCGAAAAGTCGGCTTTTAGCAACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGACCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTc                                                                         <  2:155016/251‑1 (MQ=255)
                                                                                                                                                                          gCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCg                                                                        <  1:341597/251‑1 (MQ=255)
                                                                                                                                                                            gagaAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGACAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGc                                                                      >  2:551992/1‑251 (MQ=255)
                                                                                                                                                                              gaAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCtt                                                                    >  1:31326/1‑251 (MQ=255)
                                                                                                                                                                               aaGTCGGCTTTTAGCGACAAATCCTTCCCTTTTGCCAGACGAATAAAACAGAACCCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGCCGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACACCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCttt                                                                   <  1:611219/251‑1 (MQ=255)
                                                                                                                                                                                    gtcTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAAAAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGGCGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCGCCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGc                                                              <  1:45087/249‑1 (MQ=255)
                                                                                                                                                                                                         cACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACt                                         >  1:186839/1‑251 (MQ=255)
                                                                                                                                                                                                              ttGCCAGACAAAAAAAAAAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCAGCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCTCCGCCAGCACCAGCCCTTTTTTCAGAGTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGccc                                     <  1:619414/250‑1 (MQ=255)
                                                                                                                                                                                                               tacaaaaCGAAAAAAACAGAACCCGAGGTTAATGAGCGCTTCCCCGCCGATGTGGACAACATATCTGGGCAGAGCGACATACAGTGGATCGACGCCAGGTGCGGCAGCGTCTTCATGCATCGGTTTTTCGGCGTTCATCGCAAGGGACAGCCCTGCAGAGAAAAGGCCGCATATCACCGCCAGCACCAGCCCTGTGTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCg                                    <  2:294358/244‑1 (MQ=255)
                                                                                                                                                                                                               tGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTc                                                                                       >  1:569998/1‑199 (MQ=255)
                                                                                                                                                                                                               tGCCAGACGAATAAAACAGAAACCGAGGTTAATGAGCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTc                                                                                       <  2:569998/199‑1 (MQ=255)
                                                                                                                                                                                                                      cGAATAAAACAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGt                            >  1:283646/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaaaCAGAAAACGAGGTTAATGATAGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCAACGCAAAGTGCGGCAGCGGCTCCATGAACAGGTCTTGCGGCGTTCACCGCAAAGGCCACCTCGGCACAGAAAATGACGAACATCACCGCCAGCAACAGCCCTTTTTTCAGATTGCAGTCTTCGACATTAATGCCCATCTTTCGCACGTTCAACTGCGCGGCTCGCGTTACa                        >  2:538907/1‑250 (MQ=255)
                                                                                                                                                                                                                              aCAGAAACCGAGGTTAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATcc                    >  2:183224/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaaaCCGAGGTTAATGATCGCGTCGCCGCCCATGATGAAAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACACCCCGGCAGAAAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTa                 <  1:280024/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aaCCGAGGTGAATGAGCACGTAGCCGCCCGGGAGGACACCAGAGCTTGGCGAAGCGACATACAGTGGATGGATGCCAGTTGCGGCAGCGGCTTTATGCATGTGTTTTGCGTCGTTCACCGCAAAGAACATCCCGGCAGAGAAAATCCCGCCCACAACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGTCGCTAGTTACAATCCCTACg               <  2:613190/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gAGGATAATGATAGAGTCGCCGCCCATGATGACAACATAGTTTGGCAGAGCTACATACAGAGGATCGACGCCAAGTGCGCCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCAGCGTCAGCACCAGCCCTTTTTTCAGATTGCACTCTTCGGCTTTAATGCCCATCATGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCTTACGCCa            >  2:250029/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gggTTAATGATCGCGTCGCCGCCCATGGTGCCAACAGAGCTGGGCAGAGCGACCTACAGTGGTGCGACGCCAATTGCGGCAGCGGCTTCATGCATCGTTTTTGCGGCTTTCACCGCAAAGGACATCCCGGCAGAGAAAATTCCGCACATCTCCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAAt          <  2:439971/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ttAATGATCGCGTCTCCGCCCATGATGACAACATAGCTTGGCAAAGCGACATACAGTGGATCTACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGCGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGTCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTCCAATCCCTACGCCAATccc       >  2:654822/1‑249 (MQ=255)
                                                                                                                                                                                                                                           ttAATGATCGCGACGCCGGCCATGATGACAACATAGCTTGGCAGAGCGACATCCAATGGGTCGAAGCCAAGTACGGCAGAGGCTTCCTGCATCGGTTTTGCGGAGATCATCGAAAAGGACATCCCGGCAGAGAAAATGCCGCACATCAACGCCAGCACCAGACCTTTTTTCAAATTAAAACCTTCGGCTTTAATACCCATCTTGCGCTTTGTCAACTAGCCGGCGAGGGTTACAATGCCTACGCCAATc         >  2:580892/1‑249 (MQ=255)
                                                                                                                                                                                                                                            tAATGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAgc      >  1:209889/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tGATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCAcc                                                                                                   >  1:106329/1‑155 (MQ=255)
                                                                                                                                                                                                                                               tGATCGCGTCGCCGCCCAGGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCACCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACACCACCGCCAGCAcc                                                                                                   <  2:106329/155‑1 (MQ=255)
                                                                                                                                                                                                                                               tGATCGCGTCGCCGACCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGCCGCCAAGTGCGGCAGCGGCTTCATGAATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCATCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGAGCGAGCTACAATCCCTACGACAATCAGCGcc   >  1:214472/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGccc  >  1:571630/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gATCGCGTCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcc   >  1:175339/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gATCGCGTCGCCGCCCATGATGAAAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGGGCGGCAGCGGCTTCAGGCATCGGTTTTGCGGCTTTCATCGCAAAGGACATCCCGGCAGAGAAAAGGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGcca  <  1:479696/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
GACGATACCGCCGCACAATACATAGAAACTCATATGCAGCATCCAACTGATGTAGTCATACTGCGCCGGAATGCGGGCGTGGCCCCAGGCATAGAAAAAGAATTGCAGATACCACATCAACCCGCCCAGTGTCGAGAGTAACACATTGTGAATGATCAGCGATTTTGCCAGCGAGAAGTCGGCTTTTAGCGACAAATCCTTCACTTTTGCCAGACGAATAAAACAGAAACCGAGGTTAATGATCGCGCCGCCGCCCATGATGACAACATAGCTTGGCAGAGCGACATACAGTGGATCGACGCCAAGTGCGGCAGCGGCTTCATGCATCGGTTTTGCGGCGTTCATCGCAAAGGACATCCCGGCAGAGAAAATGCCGCACATCACCGCCAGCACCAGCCCTTTTTTCAGATTGAACTCTTCGGCTTTAATGCCCATCTTGCGCTCTTTCAACTGCCCGGCGCGAGTTACAATCCCTACGCCAATCAGCGCCA  >  NC_000913/4099608‑4100098

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