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,2710GA95.8% 148.2 / NA 49G134E (GGG→GAG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/29);  minor base C (2/0);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTC  >  NC_000913/469034‑469513
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
aGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCtg                                                                                                                                                                                                                                         >  1:206330/1‑251 (MQ=255)
                 gCCACCATGGCGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGGACTGCGTGAAGCTTATTACCGTCAGCTAAGCCGGCCGCCTCCTGAGCTTTACCTGCGTCATCGCACAGGTGACCTCATCGCTCGCCCGGCCTATCACGTCGATCGCGCCGCGGTTTGCGACGGAGAAGCGGAGCTGCCGCTGGCCGGTTCACTgg                                                                                                                                                                                                                            >  1:157593/1‑247 (MQ=255)
                               ggtttccGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGTGCTGCGCtgtg                                                                                                                                                                                                           <  2:314426/245‑1 (MQ=255)
                               gATTGCCGGTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTAGCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCtgt                                                                                                                                                                                                            <  1:78313/249‑1 (MQ=255)
                                               tATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTGGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTTGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTa                                                                                                                                                                                           <  1:351202/250‑1 (MQ=255)
                                                aTCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTgcgc                                                                                                                                                                                                               >  1:288795/1‑229 (MQ=255)
                                                aTCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTgcgc                                                                                                                                                                                                               <  2:288795/229‑1 (MQ=255)
                                                        cGTTACGGCTGGCGGGTATTGCGTTTTGGTGCGTCTTATCAACTGGCTGTTAAACTGTGTTAATATTAGTTTCGACGGCGAAACCGGCAGCATCCTGAGTTTTTCCTGCGGCATCGCACCGGTGACCTCATGGCTCGGGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCTGAGAAGTTGTGATGACGCTGGTGTATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAg                                                                                                                                                                                 <  2:326015/251‑1 (MQ=255)
                                                                     gggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTACGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGCGGAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAAGCAAATTAGCTGGCAGTTGAcc                                                                                                                                                                   >  1:25074/1‑251 (MQ=255)
                                                                      ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGAACCATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGGCGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                                  >  1:138927/1‑251 (MQ=255)
                                                                      ggTATTGCTGTTTGGTGCGACTTATCAACGGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTAATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                                  <  1:401588/251‑1 (MQ=255)
                                                                       gTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCtt                                                                                                                                                                 >  1:12218/1‑251 (MQ=255)
                                                                         aTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTAt                                                                                                                                                               <  1:79393/251‑1 (MQ=255)
                                                                              tGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTcc                                                                                                                                                          >  1:51125/1‑251 (MQ=255)
                                                                                   ggTGCGTCTTATCAACTGGCTGTTGAACTGCGTTAAGTTTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGt                                                                                                                                                      <  2:287815/250‑1 (MQ=255)
                                                                                    gagcgccTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTAAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATAACGTCAATCGCGTCGCGATTGCCGCCGGCGAAGAGGCGCTGCCGCTGGGGGGTTCACTGGGGATGGGCTGGGGTGTGGTGATTATGATGTGAACGAAAACTAGCTGGCGAGTGACCTTTTTTTCCCTGTTg                                                                                                                                                    >  2:241519/7‑251 (MQ=255)
                                                                                               cAACTGGCTGTTGAACTGCGTGAAGATTATTTCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGt                                                                                                                                         <  1:286332/251‑1 (MQ=255)
                                                                                                    ggCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGAt                                                                                                                                      <  1:113171/249‑1 (MQ=255)
                                                                                                            aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATg                                                                                                                            >  1:163004/1‑251 (MQ=255)
                                                                                                            aaCTGAGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGGGATGGCGATCATg                                                                                                                            >  2:21383/1‑251 (MQ=255)
                                                                                                              atGCGTGAAGATTATTGACGTCAGCTAAGCCTGCAGCAGCCTGAGTTTTTCCAGCGTCAGCGCACCGGTGACCGCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGTCGATCATGAt                                                                                                                          <  2:124715/250‑1 (MQ=255)
                                                                                                                  gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGGTGTCTACGCAAATTAGCTGGCCGTTGACCTTATTTTCCCTGTTGCCGGTGCCAGTGATGGCGATCATGATTa                                                                                                                        >  2:140155/1‑249 (MQ=255)
                                                                                                                          tattaCCGTCAGCTAAGCCGGCAGCATCCTTAGTTTTCCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACg                                                                                                               <  2:108432/250‑1 (MQ=255)
                                                                                                                          tattaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGCAGAGGAGCTGACGCGGGTGGATTTCCTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTTGCTGGCAGTTGACCTTATTTTGCCTGTTGCCGATGCCAGTGATGGcaggtcggaagagcgtctggtg                                                                                                              >  2:295194/1‑230 (MQ=255)
                                                                                                                          tattaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGc                                                                                                                                   <  1:295194/230‑1 (MQ=255)
                                                                                                                          tattaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACg                                                                                                               <  1:327119/250‑1 (MQ=255)
                                                                                                                                cGTCAGCTCAGCCGGAAGCATCCTGAGTTTTACCTGCGTCATCGCTCCGGTGACCTCATGGCTCGCGCGACCAATGCCGTCGATCGCGTCGTGTTTGACGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCCAATTAGCTGGCCGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCCATCATGATTAAGCGCAACGGCGCTGc                                                                                                        >  2:202851/1‑251 (MQ=255)
                                                                                                                                   cAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCGCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCttt                                                                                                     <  2:217619/251‑1 (MQ=255)
                                                                                                                                        aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCAt                                                                                                 <  1:298341/250‑1 (MQ=255)
                                                                                                                                              gcagcaTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCACGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGCGAAGAGGTGCTGAAGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGACTACGCAAATTAACTGGCAGGTGACATTATTTTACCTGTTGCCGATGCAAGTGGTGGCGATCATGATTAAGTGCCACGGCGATGCTTTGCCTGAACGCt                                                                                          >  2:309449/1‑251 (MQ=255)
                                                                                                                                                   aTCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGAGGCTTTGCATGAACGCTTTAAg                                                                                     <  1:88597/251‑1 (MQ=255)
                                                                                                                                                          gTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGcac                                                                              <  1:202851/251‑1 (MQ=255)
                                                                                                                                                           ttttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGcac                                                                              >  1:113063/1‑250 (MQ=255)
                                                                                                                                                             ttACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAg                                                                            <  1:246703/250‑1 (MQ=255)
                                                                                                                                                                       cATCGCACCGGTTGCCTCAGGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGGCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGGGATGGGCTGCGCTTGGTTGATTATGATGGCTACGCAAATTTGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGGGATGGCGATCATGATTAAGCGCAACGGCGATGCTTGGCATGAACGCTTTAAGCTGGCACAGGCGGCGtttt                                                                  <  2:157593/250‑1 (MQ=255)
                                                                                                                                                                            cACCGGTGACCTCATGGATCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTGCCAGTc                                                            <  1:100495/251‑1 (MQ=255)
                                                                                                                                                                              ccGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCtt                                                          >  1:305274/1‑251 (MQ=255)
                                                                                                                                                                              ccGGTAACCTCAAGGCGCGTGGGACCAATGAAGTCGCACGCCTCGTGTTTGCAGCCGGAGAAGCGGTGCTGACGCTGGTGGGTTTACTAGTGATGGGCTGGGCTGTGTTGATTATGATGTCTACGCAAATTAGCTTGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCCGTGATGGCGATCATGATTAAGCGCACCGGCGCTGCTTTTCATGCACGCTTTAAGCTGGCACCGGCGGCGTTTTCCCGTCtt                                                          >  2:230256/1‑251 (MQ=255)
                                                                                                                                                                               cGGTTACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTTGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGGCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTa                                                         <  1:144606/251‑1 (MQ=255)
                                                                                                                                                                                              tCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                         >  1:148807/1‑251 (MQ=255)
                                                                                                                                                                                                  gCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAg                                     <  1:129827/251‑1 (MQ=255)
                                                                                                                                                                                                      ccAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTc                                 <  1:142195/251‑1 (MQ=255)
                                                                                                                                                                                                                          gtgtTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCaa             <  1:21383/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggtTTGCCGCCGGAGAATAGGTGCTGACTCTGGTGGGTTCACTGGTGATGGGCTGCGCTGTGTTGACTATGATTTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTTTTTGCGATGCCAGTGATGGCGATCATGATTAAGCGCATCGGCGATTCTTTGCATGAACGTTTTAAGCTGTCACAGGCGGCGTTTTCCAGTCTTAAGTGCCGCACCCAGTAAAGCCTGGCCAGTATCCGCGTGATCaaa            <  2:77408/250‑1 (MQ=255)
                                                                                                                                                                                                                             gtTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAg           <  1:202139/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ccgccgGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCttt      <  2:206330/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cgccgGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTg     <  2:81453/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccgGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGt   <  1:322770/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cgGAGAAGAGGTGATGAAGCTGGTGGATTCACTGGAGATGGGCTGCGCTGGGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAGGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGAACAGGCGGCGTTTTACAGGATTAATGACCGCACCCAGGAAAGCCTAACCAGTATCCGCATGATCACAGCCTTGGGtg  >  2:93264/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
AGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTC  >  NC_000913/469034‑469513

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