Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,384 G→A A172T (GCG→ACG)  mdlA → ABC transporter family protein MdlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,3840GA98.3% 188.2 / NA 59A172T (GCG→ACG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (24/33);  total (25/34)
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.62e-01

GTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTG  >  NC_000913/469138‑469627
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gTTGAACTGCGTTAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGTCAGTTGACCTTATTTTCCCTGTTTCCGATTCCAGTGATGACGAt                                                                                                                                                                                                                                                  <  1:131318/251‑1 (MQ=255)
  tGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCTTGTTGCCGATGCCAGTGGTGACGATct                                                                                                                                                                                                                                               >  2:240909/1‑250 (MQ=255)
   gaccTGCTTGAAGTTTATTACCGTCAGCTTAGCCGGCAGCATCCTTAGTTTTACCTGCGTCATCGCGCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTAGTTTGTGCCGCCGGATAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGTTCTGCGCTGTGTTGATGATGATTTCTACGCAAATTTTATCGCAGTTGCCCTTATTTTCCCTGTTCCCGATGCCAGTGATGGCGATCAt                                                                                                                                                                                                                                              <  2:458354/248‑1 (MQ=255)
        gCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGAt                                                                                                                                                                                                                                           <  2:104125/249‑1 (MQ=255)
         cGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTaa                                                                                                                                                                                                                                        <  2:235971/251‑1 (MQ=255)
          gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTGTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAg                                                                                                                                                                                                                                       <  2:128730/251‑1 (MQ=255)
                attnttACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAAc                                                                                                                                                                                                                                 <  2:412263/251‑1 (MQ=255)
                attattACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGAGCGACCAATGACGTCCATCGCGTCGTGTTTGCCGCCGGAGAAGGGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTCCGCAAATTAGCTGGCCGTTGACCTTTTTTTCCCTGTTGACGATGCCAGTGAGGACGATCATGATTAAGCGCAAc                                                                                                                                                                                                                                 >  2:110904/1‑251 (MQ=255)
                   attaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGc                                                                                                                                                                                                                              >  1:395387/1‑251 (MQ=255)
                         gTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCt                                                                                                                                                                                                                        <  2:252782/251‑1 (MQ=255)
                           cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGGCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGGGATGGGCTGCGCTGTGTTGATTATGATGTCTCCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCACCGGCGATGCttt                                                                                                                                                                                                                      >  2:144294/1‑251 (MQ=255)
                                aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCATGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCa                                                                                                                                                                                                                   <  1:168585/249‑1 (MQ=255)
                                                gAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTg                                                                                                                                                                                                   <  2:387750/249‑1 (MQ=255)
                                                                aTCGCACCGGTGACCTCATGGCTCGTGCGCCCAATGCAGTCGATCGCGTCGGGTTTGACGCCGGAGAAGGGGTGCTGTCGCTGGTGGGTTCACTGGTGATGGGCTACGCTGTGGTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTTCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTTAGCTgaaacatgccgcggttttc                                                                                                                                                                                 >  2:64239/1‑233 (MQ=255)
                                                                                catagCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                                                                                                 <  1:360324/247‑1 (MQ=255)
                                                                                   ggCTCGTGCGTCCTATGACGTCGATCGCAGTGTTGTTGCCGCCGGAGAAGGGGTTCTTACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGGTGATTATGATTTCTACGCAAATTGGCTGACAGTTGACATTATTTTCCCTTTTGCCAATGGCAGTGATGACGATCATGATTAAGCGTAACAGCGATGTTTGGCATGAACGCTTTAAGCTTACACAGGCTGTGTTTTCTAGTCTTAATGACCGCAccc                                                                                                                                                              <  2:435958/251‑1 (MQ=255)
                                                                                   ggCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGATAATGGGTGCTTACGCTGGTGTATTCACTGGTGATCGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTATCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCAcc                                                                                                                                                               <  1:224306/250‑1 (MQ=255)
                                                                                    gCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCa                                                                                                                                                             <  1:125591/251‑1 (MQ=255)
                                                                                    gCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCa                                                                                                                                                             <  1:92809/251‑1 (MQ=255)
                                                                                       ggTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                                                                                                                                                          <  2:308060/250‑1 (MQ=255)
                                                                                            gACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAATGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                                                                                     >  1:310550/1‑251 (MQ=255)
                                                                                                  tGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAcc                                                                                                                                                >  1:217739/1‑250 (MQ=255)
                                                                                                  ggACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAATCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACGCAGGAAAGCCTCACCa                                                                                                                                               <  2:269501/250‑1 (MQ=255)
                                                                                                                        tGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCt                                                                                                                          >  2:207417/1‑251 (MQ=255)
                                                                                                                            gccgGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTgg                                                                                                                      >  2:321814/1‑251 (MQ=255)
                                                                                                                                     ggggTGCTGACGCTGGTGGATTCCCTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGa                                                                                                                >  2:170059/1‑248 (MQ=255)
                                                                                                                                      gggTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGAt                                                                                                            <  1:166936/251‑1 (MQ=255)
                                                                                                                                       ggTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATc                                                                                                           <  2:396325/251‑1 (MQ=255)
                                                                                                                                          gCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGcc                                                                                                        >  1:407414/1‑251 (MQ=255)
                                                                                                                                          gCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGc                                                                                                         <  1:431296/250‑1 (MQ=255)
                                                                                                                                             gACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGGGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCGGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCTCCAGTATCCGCATGATCAAAGCCGTTGGTCTGGAAGATCGCCAg                                                                                                      <  2:30933/250‑1 (MQ=255)
                                                                                                                                                   ggtggATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCg                                                                                                <  2:374281/250‑1 (MQ=255)
                                                                                                                                                      ggATTCACTGGTGAAGGGCTGCGCTGTGTTGATTATGATGACTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTAt                                                                                            >  2:134159/1‑251 (MQ=255)
                                                                                                                                                           cACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGcc                                                                                       >  2:199212/1‑251 (MQ=255)
                                                                                                                                                                 tGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCtt                                                                                                                                                                                                          <  1:295863/130‑1 (MQ=255)
                                                                                                                                                                 tGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCtt                                                                                                                                                                                                          >  2:295863/1‑130 (MQ=255)
                                                                                                                                                                  gATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAt                                                                                 <  1:62823/250‑1 (MQ=255)
                                                                                                                                                                      ggCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGa                                                                            <  1:110904/251‑1 (MQ=255)
                                                                                                                                                                        cTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGa                                                                            >  1:48445/1‑249 (MQ=255)
                                                                                                                                                                                        atgatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCaaaaaaa                                                           >  2:197927/1‑250 (MQ=255)
                                                                                                                                                                                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGaa                                                                           <  2:220166/226‑1 (MQ=255)
                                                                                                                                                                                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGaa                                                                           >  1:220166/1‑226 (MQ=255)
                                                                                                                                                                                                 aCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGCCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCa                                                                                                       >  1:52813/1‑197 (MQ=255)
                                                                                                                                                                                                 aCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCa                                                                                                       <  2:52813/197‑1 (MQ=255)
                                                                                                                                                                                                        ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGt                                          >  2:353207/1‑251 (MQ=255)
                                                                                                                                                                                                         caGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGt                                          >  1:162049/2‑250 (MQ=255)
                                                                                                                                                                                                             tGGCAGTTGACCTTATTTTGACTGTTGCCGCTGCCAGCGATGACATACATGATTAAGCGCACCGGCGATGCTTTGCATGAGCGCTTTAATCTGACACAGGCGGCGTTTTCCCGTCTTAGTGACCGTACCCAGCGCAGCCTCACCAGTCTCCGCATGTTTTAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATTCCGAAGATCCCAGCAAATATAAAATGGGGATGGCGCGTAttt                                      >  2:209655/1‑249 (MQ=255)
                                                                                                                                                                                                                  gTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTc                                <  2:61014/251‑1 (MQ=255)
                                                                                                                                                                                                                          tATTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTcg                          <  1:333143/249‑1 (MQ=255)
                                                                                                                                                                                                                           aTTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                                                                                              >  1:60342/1‑180 (MQ=255)
                                                                                                                                                                                                                           aTTTTCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                                                                                              <  2:60342/180‑1 (MQ=255)
                                                                                                                                                                                                                            ttttCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGc                                  >  1:183796/1‑239 (MQ=255)
                                                                                                                                                                                                                            ttttCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGGTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGc                                  <  2:183796/239‑1 (MQ=255)
                                                                                                                                                                                                                             tttCCCTGTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGTAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATGCCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTcgaccc                      <  1:57435/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCAt                <  1:240909/249‑1 (MQ=255)
                                                                                                                                                                                                                                     ttGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCt              >  1:382612/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ttGCCGATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACAGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAACACATGCGGGGGGAGCGTATTGATGCTCGTTTCGACCCGACCATCt              >  2:183436/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gATGCCAGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcg        <  1:96310/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aGTGATGACGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCATTTCGACCCGACCATCTATATCGCGATTg  <  1:386001/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTG  >  NC_000913/469138‑469627

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