Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,271 G→A G134E (GGG→GAG)  mdlA → ABC transporter family protein MdlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,2710GA96.7% 169.9 / NA 62G134E (GGG→GAG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (20/38);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGAT  >  NC_000913/469027‑469521
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
gccGGGCAGACTCTGATGGTGGTCGCCACCATGGCGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTAGTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACGGCGTGAAGATTATTACCGTCAGCTAAGCCGTCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGCGAAGAGGTg                                                                                                                                                                                                                                                        <  2:273298/248‑1 (MQ=255)
acCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTCCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGc                                                                                                                                                                                                                                                       <  1:297388/250‑1 (MQ=255)
        gtgccTGTTGTGGATCGCCACCATGGTGCTGTTTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGTCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGcggg                                                                                                                                                                                                                                              <  2:218101/248‑4 (MQ=255)
        gATCCTGATGTGGATCGCCACCATGGTGCTGATTTCCGTTGTGGTTTATCTCCTGCGTTACGTCGGGCGGGTATTGCTGTTTGGTGCGGCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCtgg                                                                                                                                                                                                                                              <  2:429028/251‑1 (MQ=255)
          tCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTTGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCtggt                                                                                                                                                                                                                                             <  2:500179/250‑1 (MQ=255)
           ccTGATGTGGATCGCCACCAGGGTGCTTGTTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAATAGGGGCGGACGCtggtgg                                                                                                                                                                                                                                           <  2:72507/251‑1 (MQ=255)
                             ccaTGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGc                                                                                                                                                                                                                          <  1:147342/251‑1 (MQ=255)
                               aTGGTGCTGATTGCGGGTGGGGTTTATCTCCTTCGTTTCGTCTGGCGGGTATTGCTGTTGGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGAGTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTTCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGTTGGTGGATTCACTGGTGATGGGCTg                                                                                                                                                                                                                        <  2:549155/251‑1 (MQ=255)
                               aTGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTg                                                                                                                                                                                                                        <  1:105841/251‑1 (MQ=255)
                                 ggTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGACCAATGGCGTCGATCGCGGCGTGTTTGCCGCCGGAGAAGAGGTGGTGGCGCTGGTGGATTCACTGGTGTTggggtggg                                                                                                                                                                                                                      >  1:328780/1‑246 (MQ=255)
                                             cGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGGTTCACTGGTGGTGGGGTGCGCTGGGTTgtttt                                                                                                                                                                                                          >  1:298695/1‑247 (MQ=255)
                                                tgtgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatg                                                                                                                                                                                                        <  2:12581/250‑1 (MQ=255)
                                                        aTCTCCTGCGTTACGTCTTGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGGTGTTTACg                                                                                                                                                                                               >  1:362658/1‑251 (MQ=255)
                                                            ccTGCGTTACGTCTGGCGGGTATTGCTGGTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGTCCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTTATGGGCTGCGCTGTGTTGATTATGATGTCTACGc                                                                                                                                                                                              <  1:491896/248‑1 (MQ=255)
                                                               gCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGTGGATTCCCTGGGGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt                                                                                                                                                                                         >  1:500734/1‑250 (MQ=255)
                                                                        cTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGGGGATTCACTGGGGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAACTGGccgg                                                                                                                                                                               >  1:222191/1‑248 (MQ=255)
                                                                           gCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGAGTCGTGTTTGCCGCCGGAGAAGGGGGGATGACGGTGGTGGATTCACTGGTGATGGGCTGCGCTGAGTTGATTATGATGTCTACGGAAATTAGCTGGCAGt                                                                                                                                                                               >  1:265305/1‑248 (MQ=255)
                                                                            cGGGTATTGCTTTTTGGTGCGTCTTTTCAACGGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGAc                                                                                                                                                                           <  2:500734/251‑1 (MQ=255)
                                                                              ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTTTTTGCCGCCGGGGGAGAGGTGCTGACGCTGGGGGAATCACTGGTGATGGGCTGCGCTGTGTTGATTAAGATGTCTACGCACATTAGATGGCAggttaact                                                                                                                                                                         >  1:184964/1‑244 (MQ=255)
                                                                               gTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCtt                                                                                                                                                                        <  1:350094/251‑1 (MQ=255)
                                                                                                    tATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGGCAGCTAAGCCGGCAGCATCCCGAGTTTTTCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATg                                                                                                                                                     <  1:454818/249‑1 (MQ=255)
                                                                                                       caaccggttGTTGAACTGCGTGAAGATTATTACCGGCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGGTGACCTTATTTTCCCTGTTGCCGATGCCAGt                                                                                                                                                <  2:231921/243‑1 (MQ=255)
                                                                                                                 ttGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGCTGCCAGTGCTGGCGATc                                                                                                                                      >  1:264047/1‑251 (MQ=255)
                                                                                                                   gAACTGCGTGAAGATTATTGCCGGCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCa                                                                                                                                     <  1:45619/250‑1 (MQ=255)
                                                                                                                      cTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGGAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGCTGCCAGTGGTGGCGATCATGAt                                                                                                                                 >  1:14872/1‑251 (MQ=255)
                                                                                                                            gAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAgcg                                                                                                                           <  1:14595/251‑1 (MQ=255)
                                                                                                                            gAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGAGTCTGGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGTGGATTCACTGGTGATGTGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTCGCTAGCCGTTGCCCTTATTTTCCCTGTTGCCGATTCCAGTTATGTCGCTCATGATTaatgg                                                                                                                           >  2:453621/1‑248 (MQ=255)
                                                                                                                                  tattaCCGTCAGCTAAGCCGGCATCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGCTCGACCAAGGACGTCGATCGCGTCGTGTTTGCCGCCGGATAAGAGGTGCTGTCGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATTATTAAGCGCAACg                                                                                                                      <  2:265305/250‑1 (MQ=255)
                                                                                                                                  tattaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACg                                                                                                                      <  1:304393/250‑1 (MQ=255)
                                                                                                                                  tattaCCGTCAGCTAAACCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGTTTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACgg                                                                                                                     >  1:411779/1‑251 (MQ=255)
                                                                                                                                   attaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGAGCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGc                                                                                                                    <  1:124229/251‑1 (MQ=255)
                                                                                                                                      aCCGTCAGCTAAGCCGGCAGCATACTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCTATCGGGTAGGGTTTGCCGCCGGGGAAGGGGTGCTGACGATGGTGGATTCACGGGTGATGGGCGGGGCTGGGTTGGTTATGCTGTCTACGCAAATTAGCTGGCCGGTGGCCTTCTTTTCCCTGTTGCCGGTGCCAGTGGTGGCGATCATGATTAAGCGCCACGGCg                                                                                                                   >  2:71313/1‑249 (MQ=255)
                                                                                                                                       ccGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATg                                                                                                                <  1:147095/251‑1 (MQ=255)
                                                                                                                                        cGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGGGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGc                                                                                                               >  1:152451/1‑251 (MQ=255)
                                                                                                                                          ccacgtAAGCCGCGAGCGTCCTTAGATTTACCATCGTCAGCGCACCGGTGACCTCACGGCTCGTGCGGTCCATGACGGCGAGCGCGGCGTGGTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGGCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTTCCGATGACAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCta                                                                                                             <  2:429749/246‑2 (MQ=255)
                                                                                                                                           cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCttt                                                                                                            <  1:28272/251‑1 (MQ=255)
                                                                                                                                              cTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGCAGAGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGc                                                                                                          >  1:102947/1‑250 (MQ=255)
                                                                                                                                                     tgcatcagccTGCGTTTTCCCTGCTTCATTGCACCGGCTACCTCATGGCTCGTGCGACCAATGACTTCGGTCGCGTCGTGTTTGCTGCTGGAGAGGAGGTGCTGACTCTGGTGGATTCACGTGTGATGGGCTGCGCTGTGTTTATTATGATGTCTACGCAAATTAGCTTGCAGTTGACCGTATTTTCCCTGTTGCCGAGGCCAGTGATTGCGAGCATGATTAAGCGCAACGGCGATGCTTTGCATGAACg                                                                                                   <  2:487750/242‑1 (MQ=255)
                                                                                                                                                     gacagcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGCTCGCGACCTGGTTGCCGCCGGTGAAGAGGGGCTGACGCTGGTGGATTCCCTGGTGATGGGCTGCGCTGTGTTGATTATGATGGCTACGCACATTAGCTGGGAGTTGTCCTTATTTTTCCTTTTTCCGGTGCACGGGGTGGCGATCTTGCTTTAGCTCAACGGGGAAGCTTTgatttgacgc                                                                                                  >  2:380373/3‑242 (MQ=255)
                                                                                                                                                                gAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTg                                                                                         <  1:71313/249‑1 (MQ=255)
                                                                                                                                                                gAGTTTTACCTGCGACATCGCCCTGGTGACCGCATGGCTCTGGCGACCAATGACGGTGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGCCGCTGGTGGATTCACTGGTGATGGGCTTCGCTGTGTTGATTATGATGTCTACGCAAATGAGCTGGCAGTTGACCTGATGTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTg                                                                                         <  2:509235/249‑1 (MQ=255)
                                                                                                                                                                gAGTGTTACCTGCTTCATCGCACCGGTGATCTCATTGTGCGTTGGTACAATTACGTCGATCGCGGTGGTTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGTTTTTCTGGTTGTGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTGCCCTGGTGCTGATGCCAGTGATGGCGATCATGATTAAGCGTAATGGCGATGCGTTGCATGAACGCtgcaagcgtg                                                                                        <  2:545980/250‑10 (MQ=255)
                                                                                                                                                                    tttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGcaca                                                                                    <  1:142605/250‑1 (MQ=255)
                                                                                                                                                                       aCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAggc                                                                                 <  1:442340/250‑1 (MQ=255)
                                                                                                                                                                        ccTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAggcgg                                                                               >  1:66855/1‑251 (MQ=255)
                                                                                                                                                                           gCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGt                                                                            <  2:222191/251‑1 (MQ=255)
                                                                                                                                                                                    cACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTc                                                                   <  1:350987/251‑1 (MQ=255)
                                                                                                                                                                                        ggTGTCCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGTGATGCTTTGCATGAACGCTTTaa                                                                                             <  2:309460/221‑1 (MQ=255)
                                                                                                                                                                                        ggTGACCTCATGGCTCTGGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGATTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTaa                                                               <  2:54517/251‑1 (MQ=255)
                                                                                                                                                                                        ggTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTaa                                                                                             >  1:309460/1‑221 (MQ=255)
                                                                                                                                                                                            aCCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTGTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                           <  1:453621/251‑1 (MQ=255)
                                                                                                                                                                                                             aCCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGATTTCAAGTCTTAATGACCGCACCCAAGCAAGccc                                          >  2:466537/1‑250 (MQ=255)
                                                                                                                                                                                                              ccAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTc                                         <  1:150959/251‑1 (MQ=255)
                                                                                                                                                                                                                aattacgttgttCGCGTCGTGTGGGCCGACGGAGAAGTGTCGCTGGCGAGGGTGGATGCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTCCCCTGTTGCCGATGCCAGTGATTGCGATTATGATTAAGCTCATCGGTGATGCTTTGCATGAACGCTTTAAGCTGGCACAGTCGGCGTTTATCAGTCTTAATGACTGCTCCCAGGAAAGCCTCAc                                       <  2:469305/240‑1 (MQ=255)
                                                                                                                                                                                                                   gACGTCGATCGCGTCGGGTTTTCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATTGGTTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTTTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCTCCCAGGAAAGCCTCACCa                                     <  2:14872/250‑1 (MQ=255)
                                                                                                                                                                                                                    aCGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                                   <  1:368455/251‑1 (MQ=255)
                                                                                                                                                                                                                         gATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCg                              >  1:496045/1‑251 (MQ=255)
                                                                                                                                                                                                                         gATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCg                              >  1:353667/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttGCCGCCGGGGAAGAGGTGCTAACGCTGGTGGATTCACTGGTGATGGGATGCGCTGTGTGGATTATGATGTCTAAGCACATTCGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATAACAGTGATGGCGATCATGATTAAGCGCAACAGCGATGCTTAGCATGAACGCTTTACGCTGGCACAGGCGGCGTTTTCCAGTATTAATGCCCGCACACAAGAAAGCCTCACCAGTCTCCGCATGATCACAGc                  >  2:497244/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tGCCGCCGGATAAGTGGTGTTGACGCTGTTTTATTCGCTGGTGTTGGGGTGCGCGGTGTTGATTATGTTGGCTACGCAAATTATCTGGCAGTTGACCTTATTTTCCCTGTTTACTATGCCAGTGATGGTGATCCTGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTACGCTGGCACAGGCGGCGTTTTCCAGTCGTAATGACCGCTCTCAGGAAAGCCTCACCAGTATCCGCATGATTAAAGcc                 <  2:184964/250‑1 (MQ=255)
                                                                                                                                                                                                                                              cgGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTa                                                                >  2:244989/1‑197 (MQ=255)
                                                                                                                                                                                                                                              cgGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTa                                                                <  1:244989/197‑1 (MQ=255)
                                                                                                                                                                                                                                                      gggTGCTGTCGCTGGTGGATTCACTGGTGAGGGGCTGCGCTGTGTTGTGTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGAt  <  2:298695/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGAT  >  NC_000913/469027‑469521

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