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,2710GA88.1% 172.8 / ‑1.5 67G134E (GGG→GAG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (8/0);  new base A (26/33);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 4.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAG  >  NC_000913/469027‑469504
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
acCGGGCAGATCCTGATGTGGATCGCCACCATGGGGCTGATTGCCGTTGTGTGTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGc                                                                                                                                                                                                                                      <  2:296529/250‑1 (MQ=255)
  cGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGa                                                                                                                                                                                                                                   <  1:510583/251‑1 (MQ=255)
       aGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCtg                                                                                                                                                                                                                              <  1:428949/251‑1 (MQ=255)
           ccTGATGTGGATCGCCACCAAGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGGTTGGTGCGTCTTAGCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGGCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTTCCGCtggtg                                                                                                                                                                                                                           <  2:255568/250‑1 (MQ=255)
            cTGATGTGGATCGCCAACATGGTGCTGATTGCCGTTGGGGTTTATATCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGAGCGTCTTATCAACTGGGTGGTGAAATGCGTGAAGATTATGACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACAGGTTACCTCATGGCTCGTGGGAACAAAGGCGTCGAACGCGTCGTGTTTGGCCCGGCAGGAGGGGTGCGGACggaggtggg                                                                                                                                                                                                                         >  2:428708/1‑243 (MQ=255)
               aTGTGGGTCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTTACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTc                                                                                                                                                                                                                      <  2:419204/251‑1 (MQ=255)
                            ccaTGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGCTCGCGTCGGGTTTGCCGCCCGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGGGGTGGGc                                                                                                                                                                                                         >  2:549688/1‑251 (MQ=255)
                                      gATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCGGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTgtgt                                                                                                                                                                                               <  2:107139/251‑1 (MQ=255)
                                         tGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTCTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGCACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGGCCAATGACGTCGATCGCGTCGTGTTTTCCGCCGGAGGAGAGGTGCTGACGCTGGTGGATTCCCTGGGGATTGGCTGCGCtgtggtgg                                                                                                                                                                                            >  2:523295/1‑247 (MQ=255)
                                          gCCGTTGTGGTTTGTCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCGGCGTCATCGCACCGGTGTCATCATGGCTCGCGGGGCCAATGTCGTCGGTCGCGTCGTGGTTGCCGCCGGGGGAGGGGTGCGGCCGCGGGTGGGTTCACCGGGGTTGGGCGGCGCGGTGtttat                                                                                                                                                                                           >  2:109565/1‑248 (MQ=255)
                                             gTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTat                                                                                                                                                                                        <  1:19487/251‑1 (MQ=255)
                                                gtgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgat                                                                                                                                                                                     >  1:63322/1‑251 (MQ=255)
                                                  ggTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGGTTATGATGt                                                                                                                                                                                   >  1:216179/1‑251 (MQ=255)
                                                  ggTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGGGGATTCCCTGGTGATGGGCTGCGCTGTGGTGTTTAtgttgt                                                                                                                                                                                   >  1:15349/1‑251 (MQ=255)
                                                   gTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGGGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGCCCCATGCCGTCGATCCCCGCGGGTTTTCCGCCGGGGAAGAGGGGGTTCCGCTGGTGGGTTCACTGGTGGTGGGCCgcgccgggtttgttatggtggc                                                                                                                                                                                  >  2:344860/1‑233 (MQ=255)
                                                     ttATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGCTTATTACCGTCAGCTAAACCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGGCGTCGCTCGCGTCGTGTTTGCCGCCGGAGGAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTTTGATGTCt                                                                                                                                                                                 >  1:259905/1‑250 (MQ=255)
                                                       aTCTCCTGCTGTACGGTTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTGGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAc                                                                                                                                                                               <  2:12056/250‑1 (MQ=255)
                                                             tGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt                                                                                                                                                                        >  1:585807/1‑251 (MQ=255)
                                                              gCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTa                                                                                                                                                                       <  1:195109/251‑1 (MQ=255)
                                                                   aCGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGGGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGGGTTGATTATGATGTCTACGCAAATTAGCTgg                                                                                                                                                                  >  1:299621/1‑251 (MQ=255)
                                                                   aCGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTgg                                                                                                                                                                  <  1:100475/251‑1 (MQ=255)
                                                                         ggCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGtg                                                                                                                                                             <  2:418441/250‑2 (MQ=255)
                                                                            gggTAGTGCTGTTTGGTGCGTCGTATCAACGGGCTGTTGAACTGCGTGAAGATGATTACCGTCAGCTAAGCCGGCAGCATGCGGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCGCGTGCGACCAAGGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCGGCGTTGATTATGATGTCTACGCAGAACAGCTGGCAGTTGAcc                                                                                                                                                         <  2:98757/251‑1 (MQ=255)
                                                                             ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGCCCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGGGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                        >  1:59285/1‑251 (MQ=255)
                                                                                aTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAATATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCGGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAGTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTAt                                                                                                                                                     <  2:223733/251‑1 (MQ=255)
                                                                                 ttGCTGTTTGGTGCGTCTTATCAACTTGCTTTTGAACTGCGTGAAGATTATTTCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTAt                                                                                                                                                     <  1:193754/250‑1 (MQ=255)
                                                                                         tGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTTGCAGTTGACCTTATTTTCCCTg                                                                                                                                             <  1:446051/250‑1 (MQ=255)
                                                                                           gTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCAGCGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGCCCAATTCCGTCGATCGCCGTGTGTTTGCCGCCGGAGTAGGGGGGCGGACGCGGGGGGATTCACTGGTGATGGGCTGCGCTGGGGTGATTTTGCTGTCTACGCCCATTAGCTGCCAGGTGACCTTATTTTCCCTGTTg                                                                                                                                          >  2:475615/1‑251 (MQ=255)
                                                                                            tGCGTCTTATCAACTGGCTGTTGAACTGCGTTAAGATTATTACCGTCAGCTAAGCCGGCAGCATCATGAGTTTTACCTGCGTCATCGCAACGGTGACATCATGGCTCGTGCGCTCAATGCCGTCGATCGCGTCGTGTTTGCCGCCGGAGGAGAGGTGCTGACGCTGGTGGATTTACTGGTGGTGGGCTGCGCTGTGTTGTTTATGATGTCTACGCAAATTAACTGGCAGATTAACTTGTTTTCCCTGTTGc                                                                                                                                         >  2:207382/1‑251 (MQ=255)
                                                                                             gCGTCTTATCACCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGACCAATGACGTCGCTCGCGCCGTGGTTGCCGCCGGGGCAGAGGGGCTGACGCTGGTGGGTTCCCTTGGGATGGGCTGCGCTGGGTTGATTTTGATTTCTTAGACAATTAGCGGGCAGTTGACCTTATTTTTCCTGTTGcc                                                                                                                                        >  2:80100/1‑251 (MQ=255)
                                                                                             gCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTTCCTGTTGcc                                                                                                                                        >  1:440857/1‑251 (MQ=255)
                                                                                             gCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGcc                                                                                                                                        <  1:227810/251‑1 (MQ=255)
                                                                                                           ggCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATg                                                                                                                           <  1:486880/250‑1 (MQ=255)
                                                                                                              tGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCgg                                                                                                                       >  2:548974/1‑250 (MQ=255)
                                                                                                               gtgtaaCTGCGTGAAGATTATTACCGTCAGCTAAGACGGCAGCATCCTGAGTTTTACCTGCTTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATGAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCg                                                                                                                        <  2:50671/245‑1 (MQ=255)
                                                                                                                  gAACTGCGTGAATATTACTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCTCCTGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGTTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATGATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCa                                                                                                                    <  2:13507/250‑1 (MQ=255)
                                                                                                                  gAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGCGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAt                                                                                                                   >  2:508436/1‑251 (MQ=255)
                                                                                                                    aCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGa                                                                                                                 <  1:422038/251‑1 (MQ=255)
                                                                                                                         gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTa                                                                                                              >  1:446970/1‑249 (MQ=255)
                                                                                                                                    taCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCAGCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCg                                                                                                  <  1:350529/250‑1 (MQ=255)
                                                                                                                                           aGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTg                                                                                          >  1:178336/1‑251 (MQ=255)
                                                                                                                                               aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATg                                                                                      >  1:580928/1‑251 (MQ=255)
                                                                                                                                                  ccGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAAc                                                                                   <  1:483459/251‑1 (MQ=255)
                                                                                                                                                    ggCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGc                                                                                 >  1:416594/1‑251 (MQ=255)
                                                                                                                                                          aTCCTGAGTGTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAg                                                                           <  2:89611/251‑1 (MQ=255)
                                                                                                                                                            ccTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGCTTAAGCGCAACGGCGATGCTTTGCCTGAACGCTTTAAGcc                                                                         >  2:26632/1‑250 (MQ=255)
                                                                                                                                                               gAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTgg                                                                       <  2:580928/250‑1 (MQ=255)
                                                                                                                                                                  ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGcac                                                                    >  1:137413/1‑250 (MQ=255)
                                                                                                                                                                       ccTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAggcgg                                                              >  1:492333/1‑251 (MQ=255)
                                                                                                                                                                            gTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGttt                                                         >  1:134368/1‑251 (MQ=255)
                                                                                                                                                                             tCATCGCCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGGCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATTAACGCTTTAAGCTGGCACAGGCGGCGttt                                                         <  2:544236/250‑1 (MQ=255)
                                                                                                                                                                              caacGCACCGGTGCTCTCATGGTGCGTGCGACCAATGACTTCGAGCGCGTCGTGTTTGCCGCCGGATAAGAGGTGCTGACGCTGGTGGATTCACTTGTGATGGGCTGCTCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGATCGCTTTAAGCTGGCACAGGCGGCGgttt                                                        <  2:400519/247‑1 (MQ=255)
                                                                                                                                                                               aTCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTcc                                                      <  2:216179/251‑1 (MQ=255)
                                                                                                                                                                                      cGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCGGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTa                                               <  2:145089/251‑1 (MQ=255)
                                                                                                                                                                                       ggTGACCTCATGGCTCGTGCGACCAAAGACGTCGCTCACGTCGTGTTTGCCGCCGGAGACGGGGTGCTGAAGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCCGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCCGTCTTa                                               >  2:390775/1‑250 (MQ=255)
                                                                                                                                                                                        gtcaCCTCATGGCTCGTTCGACCAATGAAGTCGATAGAGTCGTGTTTGCCGCCGGAGAAGAGGCGCTGACCCTGGTGGATTAATTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCCGCTGACCTTATTTTCCCTGTTGCCGCTGCCAGAGATGGCGATCATGATTAAGCACAACGGCGGTGCTTTGCCTGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                             >  2:353835/4‑251 (MQ=255)
                                                                                                                                                                                          gACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGa                                           <  1:569084/251‑1 (MQ=255)
                                                                                                                                                                                          gACCTCATGGCTCGTGCGACCAAGGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATg                                            <  2:164840/250‑1 (MQ=255)
                                                                                                                                                                                                       gTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGCGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                               >  2:586721/1‑250 (MQ=255)
                                                                                                                                                                                                               aaTGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                      >  2:75878/1‑251 (MQ=255)
                                                                                                                                                                                                               aaTGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                      >  1:500529/1‑251 (MQ=255)
                                                                                                                                                                                                                aTGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAcc                     <  2:440857/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                  >  1:141721/1‑251 (MQ=255)
                                                                                                                                                                                                                       cGATCGCGTCTTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATcc              <  1:481668/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCa     >  2:178083/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCa     >  2:174359/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gtTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAg  <  1:523295/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAG  >  NC_000913/469027‑469504

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