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,2710GA89.6% 129.9 / ‑4.1 48G134E (GGG→GAG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base A (16/27);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCC  >  NC_000913/469038‑469506
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
ccTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCCGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGGCCAATGGCGGCGGTCGCGTCGTGTTTGCCGCCGGAGAAGGGGGGCTGACGCtggtgg                                                                                                                                                                                                                              >  2:151472/1‑251 (MQ=255)
ccTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGTCCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCtggtgg                                                                                                                                                                                                                              <  1:91717/251‑1 (MQ=255)
                 ccaTGGTGCTGATGTCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGc                                                                                                                                                                                                             <  1:244492/251‑1 (MQ=255)
                        ggtGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGGGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCtg                                                                                                                                                                                                     <  1:246270/249‑1 (MQ=255)
                              tGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGGAGAGGTGCTGACGCTGGTGGTTTCCCTGGTGATGGGCGGCGCTGTGTTgg                                                                                                                                                                                               >  2:382628/1‑250 (MQ=255)
                                       ggTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGGGCGACCAATGACGTCGATCGCGTCGGGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGTTTATGATGt                                                                                                                                                                                      >  1:18177/1‑251 (MQ=255)
                                        gTTTATCTCCTGCGTTACGTCGGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTTCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGGTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTc                                                                                                                                                                                     <  1:222921/251‑1 (MQ=255)
                                                   gCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTa                                                                                                                                                                          <  1:306045/251‑1 (MQ=255)
                                                     gTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGATTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGc                                                                                                                                                                        <  1:438221/251‑1 (MQ=255)
                                                       tACGTCTGGCGGGAAATGCTGTTTGGTGCGTCTTATCAACTGGATGATGAACTGCGTGAAGATTATTACAGGCAGCTAAGCCAGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCACCTGCGGCCAATGGCGTCGATCGCGTCCGGTTTGCCACCGGGAGAGGGGGGCTGACGCTGGGGGGTTTACAGGGGATGGGCTGCGCTGGGTTGATTatggattttacgcaatatagctg                                                                                                                                                                      >  2:297625/1‑231 (MQ=255)
                                                             tGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGGTGAACTTCGTGACGATTATTAACGACACGTAAAACGGCAGCAGCCTGAGTTTTACCTGCGTCATAGCACCCGTGACCTCATGGCCCGTGGGACCAATGACGTCGATAGCGGCGGTTTTGCCGACGGGGAAAGGGAGCTGACGCCGGGGGATTCACTGGGGGTGGGCCGAGGCGGGGTTGTTTTGATGTCTACGCAAAtctgacggccgtt                                                                                                                                                                >  1:38137/1‑239 (MQ=255)
                                                              ggctggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGt                                                                                                                                                                 <  1:286985/245‑1 (MQ=255)
                                                                cGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCGGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGAc                                                                                                                                                             <  1:279696/251‑1 (MQ=255)
                                                                  ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGACGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                           >  1:231800/1‑251 (MQ=255)
                                                                     gtgtggagttGGTGGTCGTATCAACTGGCTGGTGAACTGAGTGAAGATTATTTCCGTCAGCTAAGCCGGCAGCAGTATGAGTTTTACCGTCGCGATTGCTGCGGTGACCTCATGGATCGTGCGTACAATTAAGTCGATTGCGTCGTGTTTGTCGCCGGAGAAGAGGTGCTGACTCTGGGGGATTCACTGGTGATGTGCTGCGCGTTGTTGACTATGATGTCTACGCAAATTAGCTGGCAGTTGCCCTTAt                                                                                                                                                        <  2:122266/242‑1 (MQ=255)
                                                                           gTTTGGTGCGTCTTAACAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAACAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTcc                                                                                                                                                   <  1:333605/250‑1 (MQ=255)
                                                                               ggTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTGTACCTGCGTCATCGCACCGGGGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGt                                                                                                                                               <  1:417282/250‑1 (MQ=255)
                                                                                           cAACCGGCTGTTGAACTGCGTGAAGATTATTACCGACAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATAGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGGCGATCGAGTCGTGTGTGCCGCCGGAGAAGAGGTGCTGACGCTGGAGGATTCACTGGTGATGGGCGGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCGGGCAGTTGACCGTATTTTCCCTGTGGCCGAGGCCAGt                                                                                                                                  <  2:191657/251‑1 (MQ=255)
                                                                                                       gAACGGCGTGAATATTATTACCGTCAGTTAAGCCGGCAGCATCCTTAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAt                                                                                                                      <  2:322258/251‑1 (MQ=255)
                                                                                                         aCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCTTCGTGTTTGCCGCCGGAGAAGGGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATgg                                                                                                                    >  2:345462/1‑250 (MQ=255)
                                                                                                              gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTgtgt                                                                                                                                                                                                  <  1:152577/168‑1 (MQ=255)
                                                                                                              gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGACGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGGtgtgg                                                                                                                                                                                                  >  2:152577/1‑167 (MQ=255)
                                                                                                                gAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAgcg                                                                                                             >  1:13382/1‑251 (MQ=255)
                                                                                                                gAAGATTATTACCGGCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTTACGCTGGTGGATTCACTGGGGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTGTCCCTGTTGCCGATGCCAGTGATGGCGATCATGAGTAAgcg                                                                                                             <  2:327304/251‑1 (MQ=255)
                                                                                                                           ccGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATg                                                                                                  >  1:272404/1‑251 (MQ=255)
                                                                                                                               cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCGCATGGCTCGTGCGACCAATGACGTCGAGCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCttt                                                                                              <  2:18177/251‑1 (MQ=255)
                                                                                                                                  cTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAAAAGCGTCGATCGAGTCGTGTTTGCCGCCGGAGGAGGGGGGCTGACGCTGGTGGATTCACTGGCGATGGGCTGCGCTGTGTTGATTATGATGTCTAAGCAAATTAGCTGGCAGTTGACATTATTTTCACTGTTGACGATACCAGTGATGGCGCTCATGATTTAGCGCACCGGCGAAGCTTTGc                                                                                            >  1:281856/1‑250 (MQ=255)
                                                                                                                                     aGCCGGCCGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGAAGTCGATCGCGTAGTGTTTGCCGCCGGAGAAGAGGCGCTGACGCTTGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCGGATGCCAGTGATGGGGATCATGATTAAGCGCAACGGGGATGCTTTTCATg                                                                                         >  2:402977/1‑250 (MQ=255)
                                                                                                                                        cGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGCGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACg                                                                                     >  1:51828/1‑251 (MQ=255)
                                                                                                                                             gcaTCCTGAGTTTTACCTGCGTCAGCGCACCGGTGACTTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGGCTACGCAAATTAGCTGGCAGGTGGCCGGATTGTCCCTGTGGCCGATGGCAGTGATGGCGATCATGATTAAGCGCAAGGGCGATGCGTTGCAGGAAGGCTTTa                                                                                <  2:61146/251‑1 (MQ=255)
                                                                                                                                                   tgcgTTTTACCCGCGTCATCGCACCGGTGACCTCATGGCTCGTCCGACCAATGACGTCGATCGCGTCGCGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACACAACTTAGCTGGCAGTTGAACTTATTTTCCCTGATGCCGCTGCCCGTGATGGCGATCATGCTTAAGCGCAACGGCGCTTATTTTCATGAACGCTTTAAGcagg                                                                          >  2:452733/4‑248 (MQ=255)
                                                                                                                                                         ttACCTGCGACATCGCAACGGTGACCTCATGGCTCGTGCGCCCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGGTGATTATGATGTCTACGCAAATTAGCTGGCAGGTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGAGCAACGGCGGTGCTTTGCATGAACGCCTTAAGCTGGCACAgg                                                                    >  2:83596/1‑251 (MQ=255)
                                                                                                                                                              tGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCa                                                                                                                       <  1:274147/195‑1 (MQ=255)
                                                                                                                                                              tGCGTCATCGCACAGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGAGTTTGCCGCCGGAGCAGAGGTGCTGACGCTGGGGGATTCACTGGTGCTGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGCTGACCTTATTTTCCCTGTTGCCGATGCCCGTGATGGCGATCa                                                                                                                       >  2:274147/1‑195 (MQ=255)
                                                                                                                                                               gCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGgcggcg                                                               <  1:151472/250‑1 (MQ=255)
                                                                                                                                                                 gTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCAATGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGttt                                                           >  1:132899/1‑251 (MQ=255)
                                                                                                                                                                      cGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCa                                                       <  2:358862/250‑1 (MQ=255)
                                                                                                                                                                       gCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGCAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACCGGCGGCGTTTTCCAGt                                                     >  2:338119/1‑251 (MQ=255)
                                                                                                                                                                        aaCCGGTAACCTCATGGCTCGTGCGACCAATGACGTCGGTCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCGGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGGTTTTccgtc                                                     >  2:312170/2‑248 (MQ=255)
                                                                                                                                                                                   tCATGGCTCGTGCGACCAATGACGTCGATCGCGGCGTGTTTGCCGACGGAGAAGAGGTGCTGACGCTGGTGGCTTCACTGGTGATGGGCTGCGATGTGTTGATTATGAAGTCTACGCACATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGAGCAAAGGCGATGCTTTGCATGAACGCTTTAAGCTGGCAAAGGCGGCGTTTTCCAGTCTTAATGACCGc                                         >  1:217724/1‑251 (MQ=255)
                                                                                                                                                                                    cATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCg                                          <  2:231800/249‑1 (MQ=255)
                                                                                                                                                                                       ggCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAcc                                      <  1:436573/250‑1 (MQ=255)
                                                                                                                                                                                                                cgcgTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCa             <  1:338119/250‑1 (MQ=255)
                                                                                                                                                                                                                cgcgGCGTGTTTGCCGGCGGAGAAGAGGTGCTGACGCTGTAGGATTCACTGGAGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCGGGCAGTTGACCTCATTTTCCCTGTTGCAGATGCCAGTGATTGCGATGAAGATTAAGCGCAACGGTGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGTGTTGTCCAGTCTTAATGACCGCACACAGGAAAGCCTCATCAGttgccgcaa            <  2:24395/251‑9 (MQ=255)
                                                                                                                                                                                                                 gcggcgtGTTTGCCGTCGGAGAAGAGGTGATGACGCTGGTTGATTCACTGGTGATGGGCTGCGCTGTGTTGAGTATGATGTCTACGCAAATTAGCTTGCAGTTGACCGTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAGGATTAAGCGCAACGGCGAGGTTTTGCATGAACGCGTTAAGCTGGCACAGGCGGCGGTTTCTAGTTTGATTGACCGCTCTTAGGAAAGCTTCACCATTATGCGCATg           <  2:309067/247‑1 (MQ=255)
                                                                                                                                                                                                                  cgtcgtGTTTTCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGTGCTGCGCTGGGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCGGGCACATGCGTCGTTTTCCAGTCTGAATGACCGCACTCGGGAAAGCCGCACCAGTATcgtcatga          <  2:442252/251‑8 (MQ=255)
                                                                                                                                                                                                                       tgtTTGCAGCCGGCGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTATAAGCAAATTAGCTGGCAGTTGACCTTATTTTACCTGGTGCCGATGCCAGTGATGGCGCTCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACAGCACCCAGGAAAGCCTCCCCAGTATCCGCATGATCaa      >  2:285882/1‑250 (MQ=255)
                                                                                                                                                                                                                           tGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGcc  <  1:12431/250‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
CCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCC  >  NC_000913/469038‑469506

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