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

ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAG  >  NC_000913/469027‑469504
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
acCGGGCAGATCCTGATGTGGATCGCCACCATGGGGCTGATTGCCGTTGTGTGTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGc                                                                                                                                                                                                                                      <  2:296529/250‑1 (MQ=255)
  cGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGa                                                                                                                                                                                                                                   <  1:510583/251‑1 (MQ=255)
       aGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCtg                                                                                                                                                                                                                              <  1:428949/251‑1 (MQ=255)
           ccTGATGTGGATCGCCACCAAGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGGTTGGTGCGTCTTAGCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGGCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTTCCGCtggtg                                                                                                                                                                                                                           <  2:255568/250‑1 (MQ=255)
            cTGATGTGGATCGCCAACATGGTGCTGATTGCCGTTGGGGTTTATATCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGAGCGTCTTATCAACTGGGTGGTGAAATGCGTGAAGATTATGACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACAGGTTACCTCATGGCTCGTGGGAACAAAGGCGTCGAACGC‑GTCGTGTTTGGC‑CCGGCAGGAGGGGTGCGGACggaggtggg                                                                                                                                                                                                                         >  2:428708/1‑243 (MQ=255)
               aTGTGGGTCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTTACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTc                                                                                                                                                                                                                      <  2:419204/251‑1 (MQ=255)
                            ccaTGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGCTCGC‑GTCGGGTTTGCCGCCCG‑AGAAGGGGTGCTGACGCTGGTGGATTCACTGGGGGTGGGc                                                                                                                                                                                                         >  2:549688/1‑251 (MQ=255)
                                      gATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCGGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTgtgt                                                                                                                                                                                               <  2:107139/251‑1 (MQ=255)
                                         tGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTCTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGCACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGGCCAATGACGTCGATCGC‑GTCGTGTTTTCCGCCGG‑AGGAGAGGTGCTGACGCTGGTGGATTCCCTGGGGATTGGCTGCGCtgtggtgg                                                                                                                                                                                            >  2:523295/1‑247 (MQ=255)
                                          gCCGTTGTGGTTTGTCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCGGCGTCATCGCACCGGTGTCATCATGGCTCGCGGGGCCAATGTCGTCGGTCGC‑GTCGTGGTTGCCGCCGG‑GGGAGGGGTGCGGCCGCGGGTGGGTTCACCGGGGTTGGGCGGCGCGGTGtttat                                                                                                                                                                                           >  2:109565/1‑248 (MQ=255)
                                             gTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTat                                                                                                                                                                                        <  1:19487/251‑1 (MQ=255)
                                                gtgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgat                                                                                                                                                                                     >  1:63322/1‑251 (MQ=255)
                                                  ggTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGGTTATGATGt                                                                                                                                                                                   >  1:216179/1‑251 (MQ=255)
                                                  ggTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGGGGTGCTGACGCTGGGGGATTCCCTGGTGATGGGCTGCGCTGTGGTGTTTAtgttgt                                                                                                                                                                                   >  1:15349/1‑251 (MQ=255)
                                                   gTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGGGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGCCCCATGCCGTCGATCCC‑CGCGGGTTTTCCGCCGG‑GGAAGAGGGGGTTCCGCTGGTGGGTTCACTGGTGGTGGGCCgcgccgggtttgttatggtggc                                                                                                                                                                                  >  2:344860/1‑233 (MQ=255)
                                                     ttATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGCTTATTACCGTCAGCTAAACCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGGCGTCGCTCGC‑GTCGTGTTTGCCGCCGG‑AGGAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTTTGATGTCt                                                                                                                                                                                 >  1:259905/1‑250 (MQ=255)
                                                       aTCTCCTGCTGTACGGTTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTGGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAc                                                                                                                                                                               <  2:12056/250‑1 (MQ=255)
                                                             tGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAtt                                                                                                                                                                        >  1:585807/1‑251 (MQ=255)
                                                              gCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTa                                                                                                                                                                       <  1:195109/251‑1 (MQ=255)
                                                                   aCGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGGGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGGGTTGATTATGATGTCTACGCAAATTAGCTgg                                                                                                                                                                  >  1:299621/1‑251 (MQ=255)
                                                                   aCGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTgg                                                                                                                                                                  <  1:100475/251‑1 (MQ=255)
                                                                         ggCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGtg                                                                                                                                                             <  2:418441/250‑2 (MQ=255)
                                                                            gggTAGTGCTGTTTGGTGCGTCGTATCAACGGGCTGTTGAACTGCGTGAAGATGATTACCGTCAGCTAAGCCGGCAGCATGCGGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCGCGTGCGACCAAGGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCGGCGTTGATTATGATGTCTACGCAGAACAGCTGGCAGTTGAcc                                                                                                                                                         <  2:98757/251‑1 (MQ=255)
                                                                             ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGCCCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑GGAAGAGGGGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                        >  1:59285/1‑251 (MQ=255)
                                                                                aTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAATATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCGGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAGTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTAt                                                                                                                                                     <  2:223733/251‑1 (MQ=255)
                                                                                 ttGCTGTTTGGTGCGTCTTATCAACTTGCTTTTGAACTGCGTGAAGATTATTTCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTAt                                                                                                                                                     <  1:193754/250‑1 (MQ=255)
                                                                                         tGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTTGCAGTTGACCTTATTTTCCCTg                                                                                                                                             <  1:446051/250‑1 (MQ=255)
                                                                                           gTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCAGCGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGCCCAATTCCGTCGATCGCCGT‑GTGTTTGCCGCCGG‑AGTAGGGGGGCGGACGCGGGGGGATTCACTGGTGATGGGCTGCGCTGGGGTGATTTTGCTGTCTACGCCCATTAGCTGCCAGGTGACCTTATTTTCCCTGTTg                                                                                                                                          >  2:475615/1‑251 (MQ=255)
                                                                                            tGCGTCTTATCAACTGGCTGTTGAACTGCGTTAAGATTATTACCGTCAGCTAAGCCGGCAGCATCATGAGTTTTACCTGCGTCATCGCAACGGTGACATCATGGCTCGTGCGCTCAATGCCGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGGAGAGGTGCTGACGCTGGTGGATTTACTGGTGGTGGGCTGCGCTGTGTTGTTTATGATGTCTACGCAAATTAACTGGCAGATTAACTTGTTTTCCCTGTTGc                                                                                                                                         >  2:207382/1‑251 (MQ=255)
                                                                                             gCGTCTTATCACCTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGACCAATGACGTCGCTCGC‑GCCGTGGTTGCCGCCGG‑GGCAGAGGGGCTGACGCTGGTGGGTTCCCTTGGGATGGGCTGCGCTGGGTTGATTTTGATTTCTTAGACAATTAGCGGGCAGTTGACCTTATTTTTCCTGTTGcc                                                                                                                                        >  2:80100/1‑251 (MQ=255)
                                                                                             gCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTTCCTGTTGcc                                                                                                                                        >  1:440857/1‑251 (MQ=255)
                                                                                             gCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGcc                                                                                                                                        <  1:227810/251‑1 (MQ=255)
                                                                                                           ggCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATg                                                                                                                           <  1:486880/250‑1 (MQ=255)
                                                                                                              tGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCgg                                                                                                                       >  2:548974/1‑250 (MQ=255)
                                                                                                               gtgtaaCTGCGTGAAGATTATTACCGTCAGCTAAGACGGCAGCATCCTGAGTTTTACCTGCTTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATGAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCg                                                                                                                        <  2:50671/245‑1 (MQ=255)
                                                                                                                  gAACTGCGTGAATATTACTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCTCCTGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGTTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATGATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCa                                                                                                                    <  2:13507/250‑1 (MQ=255)
                                                                                                                  gAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑CGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAt                                                                                                                   >  2:508436/1‑251 (MQ=255)
                                                                                                                    aCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGa                                                                                                                 <  1:422038/251‑1 (MQ=255)
                                                                                                                         gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTa                                                                                                              >  1:446970/1‑249 (MQ=255)
                                                                                                                                    taCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCAGCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCg                                                                                                  <  1:350529/250‑1 (MQ=255)
                                                                                                                                           aGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTg                                                                                          >  1:178336/1‑251 (MQ=255)
                                                                                                                                               aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATg                                                                                      >  1:580928/1‑251 (MQ=255)
                                                                                                                                                  ccGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAAc                                                                                   <  1:483459/251‑1 (MQ=255)
                                                                                                                                                    ggCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGc                                                                                 >  1:416594/1‑251 (MQ=255)
                                                                                                                                                          aTCCTGAGTGTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAg                                                                           <  2:89611/251‑1 (MQ=255)
                                                                                                                                                            ccTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑GGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGCTTAAGCGCAACGGCGATGCTTTGCCTGAACGCTTTAAGcc                                                                         >  2:26632/1‑250 (MQ=255)
                                                                                                                                                               gAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTgg                                                                       <  2:580928/250‑1 (MQ=255)
                                                                                                                                                                  ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGcac                                                                    >  1:137413/1‑250 (MQ=255)
                                                                                                                                                                       ccTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAggcgg                                                              >  1:492333/1‑251 (MQ=255)
                                                                                                                                                                            gTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGttt                                                         >  1:134368/1‑251 (MQ=255)
                                                                                                                                                                             tCATCGCCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GGCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATTAACGCTTTAAGCTGGCACAGGCGGCGttt                                                         <  2:544236/250‑1 (MQ=255)
                                                                                                                                                                              caacGCACCGGTGCTCTCATGGTGCGTGCGACCAATGACTTCGAGCGC‑GTCGTGTTTGCCGCCGG‑ATAAGAGGTGCTGACGCTGGTGGATTCACTTGTGATGGGCTGCTCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGATCGCTTTAAGCTGGCACAGGCGGCGgttt                                                        <  2:400519/247‑1 (MQ=255)
                                                                                                                                                                               aTCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTcc                                                      <  2:216179/251‑1 (MQ=255)
                                                                                                                                                                                      cGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCGGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTa                                               <  2:145089/251‑1 (MQ=255)
                                                                                                                                                                                       ggTGACCTCATGGCTCGTGCGACCAAAGACGTCGCTCAC‑GTCGTGTTTGCCGCCGG‑AGACGGGGTGCTGAAGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCCGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCCGTCTTa                                               >  2:390775/1‑250 (MQ=255)
                                                                                                                                                                                        gtcaCCTCATGGCTCGTTCGACCAATGAAGTCGATAGA‑GTCGTGTTTGCCGCCGG‑AGAAGAGGCGCTGACCCTGGTGGATTAATTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCCGCTGACCTTATTTTCCCTGTTGCCGCTGCCAGAGATGGCGATCATGATTAAGCACAACGGCGGTGCTTTGCCTGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                             >  2:353835/4‑251 (MQ=255)
                                                                                                                                                                                          gACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGa                                           <  1:569084/251‑1 (MQ=255)
                                                                                                                                                                                          gACCTCATGGCTCGTGCGACCAAGGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATg                                            <  2:164840/250‑1 (MQ=255)
                                                                                                                                                                                                       gTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGCGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGa                               >  2:586721/1‑250 (MQ=255)
                                                                                                                                                                                                               aaTGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                      >  2:75878/1‑251 (MQ=255)
                                                                                                                                                                                                               aaTGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                      >  1:500529/1‑251 (MQ=255)
                                                                                                                                                                                                                aTGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAcc                     <  2:440857/251‑1 (MQ=255)
                                                                                                                                                                                                                   aCGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                  >  1:141721/1‑251 (MQ=255)
                                                                                                                                                                                                                       cGATCGC‑GTCTTGTTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATcc              <  1:481668/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCa     >  2:178083/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gtgtTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCa     >  2:174359/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gtTTGCCGCCGG‑AGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAg  <  1:523295/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
ACCGGGCAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGC‑GTCGTGTTTGCCGCCGG‑AGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAG  >  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: ‑