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

GTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTG  >  NC_000913/469138‑469627
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gTTGAACTGCGTTAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGTCAGTTGACCTTATTTTCCCTGTTTCCGATTCCAGTGATGACGAt                                                                                                                                                                                                                                                  <  1:131318/251‑1 (MQ=255)
  tGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCTTGTTGCCGATGCCAGTGGTGACGAT‑ct                                                                                                                                                                                                                                               >  2:240909/1‑250 (MQ=255)
   gaccTGCTTGAAGTTTATTACCGTCAGCTTAGCCGGCAGCATCCTTAGTTTTACCTGCGTCATCGCGCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTAGTTTGTGCCGCCGGATAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGTTCTGCGCTGTGTTGATGATGATTTCTACGCAAATTTTATCGCAGTTGCCCTTATTTTCCCTGTTCCCGATGCCAGTGATGGCGAT‑CAt                                                                                                                                                                                                                                              <  2:458354/248‑1 (MQ=255)
        gCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGAt                                                                                                                                                                                                                                           <  2:104125/249‑1 (MQ=255)
         cGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTaa                                                                                                                                                                                                                                        <  2:235971/251‑1 (MQ=255)
          gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTGTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAg                                                                                                                                                                                                                                       <  2:128730/251‑1 (MQ=255)
                attnttACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAAc                                                                                                                                                                                                                                 <  2:412263/251‑1 (MQ=255)
                attattACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGAGCGACCAATGACGTCCATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTCCGCAAATTAGCTGGCCGTTGACCTTTTTTTCCCTGTTGACGATGCCAGTGAGGACGAT‑CATGATTAAGCGCAAc                                                                                                                                                                                                                                 >  2:110904/1‑251 (MQ=255)
                   attaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGc                                                                                                                                                                                                                              >  1:395387/1‑251 (MQ=255)
                         gTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCt                                                                                                                                                                                                                        <  2:252782/251‑1 (MQ=255)
                           cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GGCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGGGATGGGCTGCGCTGTGTTGATTATGATGTCTCCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCACCGGCGATGCttt                                                                                                                                                                                                                      >  2:144294/1‑251 (MQ=255)
                                aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCATGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCa                                                                                                                                                                                                                   <  1:168585/249‑1 (MQ=255)
                                                gAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTg                                                                                                                                                                                                   <  2:387750/249‑1 (MQ=255)
                                                                aTCGCACCGGTGACCTCATGGCTCGTGCGCCCAATGCAGTCGATCGC‑GTCGGGTTTGACGCCGGAGAAGGGGTGCTGTCGCTGGTGGGTTCACTGGTGATGGGCTACGCTGTGGTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTTCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTTAGCTgaaacatgccgcggttttc                                                                                                                                                                                 >  2:64239/1‑233 (MQ=255)
                                                                                catagCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCa                                                                                                                                                                 <  1:360324/247‑1 (MQ=255)
                                                                                   ggCTCGTGCGTCCTATGACGTCGATCGCAGT‑GTTGTTGCCGCCGGAGAAGGGGTTCTTACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGGTGATTATGATTTCTACGCAAATTGGCTGACAGTTGACATTATTTTCCCTTTTGCCAATGGCAGTGATGACGAT‑CATGATTAAGCGTAACAGCGATGTTTGGCATGAACGCTTTAAGCTTACACAGGCTGTGTTTTCTAGTCTTAATGACCGCAccc                                                                                                                                                              <  2:435958/251‑1 (MQ=255)
                                                                                   ggCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGATAATGGGTGCTTACGCTGGTGTATTCACTGGTGATCGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTATCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCAcc                                                                                                                                                               <  1:224306/250‑1 (MQ=255)
                                                                                    gCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCa                                                                                                                                                             <  1:125591/251‑1 (MQ=255)
                                                                                    gCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCa                                                                                                                                                             <  1:92809/251‑1 (MQ=255)
                                                                                       ggTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                                                                                                                                                          <  2:308060/250‑1 (MQ=255)
                                                                                            gACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAATGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                                                                                     >  1:310550/1‑251 (MQ=255)
                                                                                                  tGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAcc                                                                                                                                                >  1:217739/1‑250 (MQ=255)
                                                                                                  ggACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAATCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACGCAGGAAAGCCTCACCa                                                                                                                                               <  2:269501/250‑1 (MQ=255)
                                                                                                                        tGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCt                                                                                                                          >  2:207417/1‑251 (MQ=255)
                                                                                                                            gccgGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTgg                                                                                                                      >  2:321814/1‑251 (MQ=255)
                                                                                                                                     ggggTGCTGACGCTGGTGGATTCCCTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGa                                                                                                                >  2:170059/1‑248 (MQ=255)
                                                                                                                                      gggTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGAt                                                                                                            <  1:166936/251‑1 (MQ=255)
                                                                                                                                       ggTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATc                                                                                                           <  2:396325/251‑1 (MQ=255)
                                                                                                                                          gCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGcc                                                                                                        >  1:407414/1‑251 (MQ=255)
                                                                                                                                          gCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGc                                                                                                         <  1:431296/250‑1 (MQ=255)
                                                                                                                                             gACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGGGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCGGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCTCCAGTATCCGCATGATCAAAGCCGTTGGTCTGGAAGATCGCCAg                                                                                                      <  2:30933/250‑1 (MQ=255)
                                                                                                                                                   ggtggATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCg                                                                                                <  2:374281/250‑1 (MQ=255)
                                                                                                                                                      ggATTCACTGGTGAAGGGCTGCGCTGTGTTGATTATGATGACTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTAt                                                                                            >  2:134159/1‑251 (MQ=255)
                                                                                                                                                           cACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGcc                                                                                       >  2:199212/1‑251 (MQ=255)
                                                                                                                                                                 tGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCtt                                                                                                                                                                                                          <  1:295863/130‑1 (MQ=255)
                                                                                                                                                                 tGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCtt                                                                                                                                                                                                          >  2:295863/1‑130 (MQ=255)
                                                                                                                                                                  gATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGAt                                                                                 <  1:62823/250‑1 (MQ=255)
                                                                                                                                                                      ggCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGa                                                                            <  1:110904/251‑1 (MQ=255)
                                                                                                                                                                        cTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGa                                                                            >  1:48445/1‑249 (MQ=255)
                                                                                                                                                                                        atgatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCaaaaaaa                                                           >  2:197927/1‑250 (MQ=255)
                                                                                                                                                                                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGaa                                                                           <  2:220166/226‑1 (MQ=255)
                                                                                                                                                                                                tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGaa                                                                           >  1:220166/1‑226 (MQ=255)
                                                                                                                                                                                                 aCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGCCGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCa                                                                                                       >  1:52813/1‑197 (MQ=255)
                                                                                                                                                                                                 aCGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCa                                                                                                       <  2:52813/197‑1 (MQ=255)
                                                                                                                                                                                                        ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGt                                          >  2:353207/1‑251 (MQ=255)
                                                                                                                                                                                                         caGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGt                                          >  1:162049/2‑250 (MQ=255)
                                                                                                                                                                                                             tGGCAGTTGACCTTATTTTGACTGTTGCCGCTGCCAGCGATGAC‑ATACATGATTAAGCGCACCGGCGATGCTTTGCATGAGCGCTTTAATCTGACACAGGCGGCGTTTTCCCGTCTTAGTGACCGTACCCAGCGCAGCCTCACCAGTCTCCGCATGTTTTAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATTCCGAAGATCCCAGCAAATATAAAATGGGGATGGCGCGTAttt                                      >  2:209655/1‑249 (MQ=255)
                                                                                                                                                                                                                  gTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTc                                <  2:61014/251‑1 (MQ=255)
                                                                                                                                                                                                                          tATTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTcg                          <  1:333143/249‑1 (MQ=255)
                                                                                                                                                                                                                           aTTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                                                                                              >  1:60342/1‑180 (MQ=255)
                                                                                                                                                                                                                           aTTTTCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                                                                                              <  2:60342/180‑1 (MQ=255)
                                                                                                                                                                                                                            ttttCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGc                                  >  1:183796/1‑239 (MQ=255)
                                                                                                                                                                                                                            ttttCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGGTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGc                                  <  2:183796/239‑1 (MQ=255)
                                                                                                                                                                                                                             tttCCCTGTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGTAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATGCCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTcgaccc                      <  1:57435/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gTTGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCAt                <  1:240909/249‑1 (MQ=255)
                                                                                                                                                                                                                                     ttGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCt              >  1:382612/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ttGCCGATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACAGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAACACATGCGGGGGGAGCGTATTGATGCTCGTTTCGACCCGACCATCt              >  2:183436/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gATGCCAGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcg        <  1:96310/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aGTGATGACGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGACACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCATTTCGACCCGACCATCTATATCGCGATTg  <  1:386001/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGAT‑CATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTG  >  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: ‑