Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,452 G→T F742L (TTC→TTA fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,4520GT97.1% 217.8 / ‑2.0 69F742L (TTC→TTAfecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (42/25);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTC  >  NC_000913/4514216‑4514683
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
aaaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTcac                                                                                                                                                                                                                           <  2:655180/251‑1 (MQ=255)
 aaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTcaca                                                                                                                                                                                                                          <  1:647653/251‑1 (MQ=255)
   cgggCCGTGGTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAACACCAGCAGCAGGCCTGCAAAAAGAAAACGGAGAAATGCCAACATAAACACAGTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACAcc                                                                                                                                                                                                                        <  2:241438/248‑1 (MQ=255)
      gCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGa                                                                                                                                                                                                                      >  2:702305/1‑250 (MQ=255)
      gCCGGGTTCGGCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTCCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGaa                                                                                                                                                                                                                     <  2:393144/251‑1 (MQ=255)
                   tGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTc                                                                                                                                                                                                           >  1:642650/1‑248 (MQ=255)
                   tGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTc                                                                                                                                                                                                           >  1:288323/1‑248 (MQ=255)
                           ggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCAt                                                                                                                                                                                                <  2:101974/251‑1 (MQ=255)
                                        tGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCTCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCaa                                                                                                                                                                                     >  2:67586/1‑249 (MQ=255)
                                          gCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTc                                                                                                                                                                                 >  1:47510/1‑251 (MQ=255)
                                          gCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTc                                                                                                                                                                                 >  1:489884/1‑251 (MQ=255)
                                             gaacACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCGTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCAt                                                                                                                                                                               <  1:236449/247‑1 (MQ=255)
                                                  ccAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGTGCCAAAGTCATAcaccc                                                                                                                                                                         >  2:391002/1‑247 (MQ=255)
                                                  ccAGCAGCAGGACTGCAAAAAGAAAACGGATTAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATGCAGCGTGCGCGGGTGGCCTGCATAGATGCCGTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACgtca                                                                                                                                                                         <  2:482455/251‑4 (MQ=255)
                                                            gCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTATTTAACACAGAAAGCCAGGTTCAGATCAACCCACTGCGGGACAAAGTAATACGCAACGCGTGCGcc                                                                                                                                                               >  1:127744/1‑251 (MQ=255)
                                                                   aaaaGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGc                                                                                                                                                        >  1:529003/1‑251 (MQ=255)
                                                                         aaaCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGa                                                                                                                                                    >  1:69606/1‑249 (MQ=255)
                                                                                 aaaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAAt                                                                                                                                          <  1:213401/251‑1 (MQ=255)
                                                                                     gCCAACATAATCAAATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCCGAACTTCCACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATCAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGCATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCCGAGCATGAAGCCGTGAATGCgg                                                                                                                                      >  2:296759/1‑251 (MQ=255)
                                                                                       cAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGcc                                                                                                                                    >  1:371039/1‑251 (MQ=255)
                                                                                         aCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGAt                                                                                                                                                                                                       >  2:485040/1‑182 (MQ=255)
                                                                                         aCATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGAt                                                                                                                                                                                                       <  1:485040/182‑1 (MQ=255)
                                                                                              aTCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTg                                                                                                                             >  1:326749/1‑251 (MQ=255)
                                                                                                cacaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGCCCCATGCATATACAGCGTGCGCGGCTGGCCTGCATAGCTGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGAGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGCGCATGAAGCCGGGAATGCGGCCGGTACTGcc                                                                                                                           >  1:29591/1‑251 (MQ=255)
                                                                                                cacaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGGCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACACCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGcc                                                                                                                           >  2:228439/1‑251 (MQ=255)
                                                                                                cacaTTCCAACTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCATTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCAGGCGCCCAAGAGCATGAAGCCGGGAATGCGGCCGGTACTGcc                                                                                                                           >  2:537809/1‑251 (MQ=255)
                                                                                                   aTTCCAGCTAAAAGCCCGTCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCGTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTc                                                                                                                        <  2:252289/251‑1 (MQ=255)
                                                                                                      ccAGCGAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGGCATAAGAGCGGATGAAATAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGgc                                                                                                                     <  2:410713/251‑1 (MQ=255)
                                                                                                          cTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGCTGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGctc                                                                                                                 >  2:295456/1‑251 (MQ=255)
                                                                                                           tAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGctct                                                                                                                >  1:180111/1‑251 (MQ=255)
                                                                                                             aaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCTTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGctctat                                                                                                              <  2:449433/251‑3 (MQ=255)
                                                                                                              aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGGGCGCCCCAGAGCATGAAGCCGGGAATGCGTCCGGGACTGCCGTCGGCGCTCTCtt                                                                                                             >  2:164843/1‑251 (MQ=255)
                                                                                                                gCCCGGCAATCCGGGCGTTCACACATCAGACCTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGGCTAAGATGTTTTTCACACCGAATGCCAGATTCAGATCTGCCCTCTGCGGGCCAAAGTCATACGCCCCGCGTGCGCCCCCGCGCATTAAGCCGGGGATGCGGACGGCACTTCCGTCGGCGCTCTCTTTc                                                                                                           >  2:616538/1‑251 (MQ=255)
                                                                                                                 cccGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTcc                                                                                                          >  2:177750/1‑250 (MQ=255)
                                                                                                                    ggccaGCCGGGCGTTAACACATCACAACTTCAACGATCCCGGCATATACAGCGTGCTCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGGAGTCCTGGTCTAAGATGTTTTTCACACCGAATTCCATGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCTGTACTGCCGTCTGCGCTCTCTTTCACCg                                                                                                       <  2:678870/247‑1 (MQ=255)
                                                                                                                       aaGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGAACTGCCGTCGGCGCTCTCTTTCACCGTAt                                                                                                    >  1:398703/1‑251 (MQ=255)
                                                                                                                       aaGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCAGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCCCCGttt                                                                                                    >  2:300820/1‑249 (MQ=255)
                                                                                                                           cGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGc                                                                                                >  2:353162/1‑251 (MQ=255)
                                                                                                                            gggCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCg                                                                                               >  1:248246/1‑251 (MQ=255)
                                                                                                                               cGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGtt                                                                                             >  1:224235/1‑250 (MQ=255)
                                                                                                                                       caTCAGAACTTCAACGACCCCTGCATATCCAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAAAGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGAGTTCTCCGCaa                                                                                     >  2:33801/1‑250 (MQ=255)
                                                                                                                                       caTCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCaaa                                                                                    <  1:7480/251‑1 (MQ=255)
                                                                                                                                                      cccccttGCATATACCGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGaa                                                                     <  2:322215/245‑1 (MQ=255)
                                                                                                                                                             gCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTg                                                              >  1:351325/1‑251 (MQ=255)
                                                                                                                                                               atatACCGCGTGCGCGGTTGGACTGCATAGATGCCTTGGTTTTTGTCGTCATAAGAGCGGATGAAGTAGTCTTTGTCTAAGATGTTTTTCACTCCGAATGCCAGGTTCTGATCTGCCATCTGCGGGCCAAAGTCATGCGCCCCGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGCCCGGAAATCGCTGtt                                                            <  2:501301/251‑1 (MQ=255)
                                                                                                                                                                  taCAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTc                                                           >  1:387511/1‑249 (MQ=255)
                                                                                                                                                                    cAGCGTGCGAGGCTGGCCTGCCTAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGCAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCTTTATCCGCAACCTGGCTGGACTGGaactaccttttcggat                                                       >  2:296717/1‑237 (MQ=255)
                                                                                                                                                                        gTGCGCTGCTGGCCTGCAGTGATGCCTTTGTGTGTGGCGTCATAAGAGCGGATGAAGTAGGCCTGGTCTAAGATGTTTTTCACCCCGAATGCCCGGTTCAGATCGGCCATCTGCGGGCCAAAGTCATCCGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGAGTTAACCGCAAACTGGCTGGACTGGACATCGCTGTTCAGATTGaa                                                   <  2:182770/251‑1 (MQ=255)
                                                                                                                                                                         tGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCCCCGTATTGGCGTTATCCGCAACCTGGCTGGACTGGACATCGCTGTTCCGATTTAAc                                                  >  1:420773/1‑251 (MQ=255)
                                                                                                                                                                               gcgggCCTGCATAGATGCCTTTGTTGTTGGCGTCATAAGAGCGGATGAAGTAGGCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTGCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGGATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAg                                            <  2:169158/248‑1 (MQ=255)
                                                                                                                                                                                             aTGCCTTTGTTGTTGTCGTCATAAGAGCGGATGACGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCAGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTTGACGTCCAGTTTCCTGGATTGtt                              >  2:663897/1‑250 (MQ=255)
                                                                                                                                                                                               gCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTAATGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGCATGATAGTCCAGTTTCCTGGCTTGTAGt                            >  2:529322/1‑251 (MQ=255)
                                                                                                                                                                                                   ttgttgttgTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTcc                          <  1:353162/249‑1 (MQ=255)
                                                                                                                                                                                                    tgttgttgtCGTCATAAGGGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTCCACACCGAAGGCCAGGTTCAGATCTGCCAGCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGTCTTGTAGTCCAc                        <  2:272849/250‑1 (MQ=255)
                                                                                                                                                                                                      ttgttgaCGTCATAAGAACGGATGAAGTCGTCCTGGTCTAAGATGTTTTTCACACGGAATGCCAGGTTCAGTTCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCTCCCGAGCATGAAGCAGGGAATGCGGCCGGTACTGCGGTCGGCGCTCTCTTTACCCTTTTTGGCGATATCCGCAACCTGGCTGTACTGGCCATCGCTGTTCAGATTGACAGGCTAGTTTCCTAGCTTGTCGTCAACGcc                     >  2:392504/1‑251 (MQ=255)
                                                                                                                                                                                                        gttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGTGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTTGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCa                   <  1:651604/251‑1 (MQ=255)
                                                                                                                                                                                                        gttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCa                   <  1:67586/251‑1 (MQ=255)
                                                                                                                                                                                                        gttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAAGGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCa                   <  2:288323/251‑1 (MQ=255)
                                                                                                                                                                                                           tgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGc                 >  1:427551/1‑251 (MQ=255)
                                                                                                                                                                                                               gtcATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCCCACCGAATGACAGGTTCAGATCTGCCATCTGCGGGCCCAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTTGCGTTCTCCGCAACCTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCAAGCCCcgcgagc             >  2:34908/1‑248 (MQ=255)
                                                                                                                                                                                                               gtcATAAGAGCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCATGGCTTGTAGTCCACGCCCAGCGTGc             >  1:615636/1‑251 (MQ=255)
                                                                                                                                                                                                                    cagagCGGATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGAGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGGGCCttt         >  2:202043/2‑250 (MQ=255)
                                                                                                                                                                                                                     agcggggATGAAGTAGTACTGGTCTAAGAGGTTTGTCAAACCGCATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCAAGAAGCCGGGAATGCGGCCGGTACTGCCGGCGGCGCTACCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAATGTCCAGTTTCCTGGCGTGTAGCCCACGCCCAGCGTGCCTTTAt       <  2:47510/246‑1 (MQ=255)
                                                                                                                                                                                                                          ggATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGt                                           <  1:548728/210‑1 (MQ=255)
                                                                                                                                                                                                                          ggATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCCCGCCCAGCGTGCCTTTATGTTTc  >  2:61229/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGATCAGCTTGAACGTCCAGt                                           >  2:548728/1‑210 (MQ=255)
                                                                                                                                                                                                                          ggATGAAGTAGTCCTGGTCTAAGATGTTTTTCACACCGAATGCCAGGGTCAGATCTGCCATCTACGGGCCAAAGTCATACGCAACGCGTGCGCCCCAGAACCTGAAGCCGGCAATGCGGCCGGTACTGCCGTCGGCGCGCTCTTTCACCGTATTGGCGTTTTCCGCAACCTGGCTGGCCTGGAAGTCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGGAGTTCACCACCAGCGGTACTTTATGtttt  >  2:698172/1‑250 (MQ=255)
                                                                                                                                                                                                                             tGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCa                                                          >  1:31681/1‑192 (MQ=255)
                                                                                                                                                                                                                             tGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCGGTTCa                                                          <  2:31681/192‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                       
AAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTC  >  NC_000913/4514216‑4514683

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: