Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 238,472 C→G 7.1% A89P (GCC→CCC)  ykfM ← lethality reduction protein, putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913238,4720CG 7.1% 180.1 / 3.5 54A89P (GCC→CCC) ykfMlethality reduction protein, putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base G (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAAGAAAGGACATGGAATCCGGAGTTAAAAACTGTTTGTTTTAGAATATTTAATTATTTAACTGCTGAATCTTCCTTGCGAGAAATATATCCATTATTACATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATA  >  NC_000913/238202‑238750
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
tGCAAGAAAGGACATGGAATCCGGAGTTAAAAACTGTTTGTTTTAGAATATTTAATTATTTAACTGCTGAATCTTCCTTGCGAGAAATATATCCATTATTACATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACt                                                                                                                                                                                                                                                                              >  1:1148864/1‑281 (MQ=255)
                                              aaTATTTAATTATTTAACTGCTGAATCTTCCTTGCGAGAAATATATCCATTATTACATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTaaaa                                                                                                                                                                                                                                >  2:713607/1‑281 (MQ=255)
                                                       ttattTAACTGCTGAATCTTCCTTTGGAGAAAAATATCCATTATTACATTTCATGCCATTTTTATAATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATAAGTATTGACATATCCTGTGACAAAAGGGGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTATTATCGTAAGTAAAATTACATGAat                                                                                                                                                                                                                      <  2:1148863/281‑1 (MQ=255)
                                                                  ctgAATCTTCCTTGCGAGAAATATATCCATTATTACATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAg                                                                                                                                                                                                            >  1:191747/1‑281 (MQ=255)
                                                                                    aaaTTTTTCCATTCTTACATTTCATGCCATTTTTAATATAGATTTCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGAAAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCg                                                                                                                                                                                          <  2:191747/281‑1 (MQ=255)
                                                                                          atCCATTATTGCATTTCATGCCATTTTTAATAGAGATTCCTCCTGGGACAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGcac                                                                                                                                                                                    <  1:621954/281‑1 (MQ=255)
                                                                                                   taCATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAAc                                                                                                                                                                           >  1:1137613/1‑281 (MQ=255)
                                                                                                            tGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCgtaagtaa                                                                                                                                                                                                                                  >  2:810762/1‑217 (MQ=255)
                                                                                                            tGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCgtaagtaa                                                                                                                                                                                                                                  <  1:810763/217‑1 (MQ=255)
                                                                                                                             aTTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCg                                                                                                                                                                                                                                                                                       <  2:142525/147‑1 (MQ=255)
                                                                                                                             aTTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCg                                                                                                                                                                                                                                                                                       >  1:142525/1‑147 (MQ=255)
                                                                                                                               tGCTCCTGGTAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTTTTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAATCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTgg                                                                                                                                               <  1:384882/281‑1 (MQ=255)
                                                                                                                                 cTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAATCACAGAATCAACATTCTACAATCATAAAATATTAACGTGGAg                                                                                                                                             >  2:24754/1‑281 (MQ=255)
                                                                                                                                 cTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAg                                                                                                                                             <  1:713607/281‑1 (MQ=255)
                                                                                                                                 cTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAg                                                                                                                                             >  1:469871/1‑281 (MQ=255)
                                                                                                                                           aaTTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTttatta                                                                                                                                   <  2:1137612/281‑1 (MQ=255)
                                                                                                                                             ttCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTattatt                                                                                                                                 <  2:469871/281‑1 (MQ=255)
                                                                                                                                                                                     atatCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCAt                                                                                         >  2:985116/1‑281 (MQ=255)
                                                                                                                                                                                        tCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACa                                                                                                                                                                          <  2:464948/197‑1 (MQ=255)
                                                                                                                                                                                        tCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACa                                                                                                                                                                          >  1:464948/1‑197 (MQ=255)
                                                                                                                                                                                           tgtgACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGGGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGcac                                                                                                                                                                                    >  2:262942/1‑184 (MQ=255)
                                                                                                                                                                                           tgtgACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGGGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGcac                                                                                                                                                                                    <  1:262942/184‑1 (MQ=255)
                                                                                                                                                                                             tgACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTtta                                                                                                                                      <  2:177231/228‑1 (MQ=255)
                                                                                                                                                                                             tgACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTtta                                                                                                                                      >  1:177231/1‑228 (MQ=255)
                                                                                                                                                                                                       aGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAaata                                                                       >  1:953681/1‑281 (MQ=255)
                                                                                                                                                                                                                    gCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATaaa                                                          >  1:338453/1‑281 (MQ=255)
                                                                                                                                                                                                                      tATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGa                                                                                                                                                                                                                                                    <  2:692879/93‑1 (MQ=255)
                                                                                                                                                                                                                      tATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGa                                                                                                                                                                                                                                                    >  1:692879/1‑93 (MQ=255)
                                                                                                                                                                                                                             aTAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGGGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTa                                                                                                         <  1:713616/225‑1 (MQ=255)
                                                                                                                                                                                                                             aTAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGGGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTACGGTa                                                                                                         >  2:713616/1‑225 (MQ=255)
                                                                                                                                                                                                                                   aTTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATAttt                                                                                                                                                     >  1:313678/1‑175 (MQ=255)
                                                                                                                                                                                                                                   aTTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATAttt                                                                                                                                                     <  2:313678/175‑1 (MQ=255)
                                                                                                                                                                                                                                     tAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTcc                                                                                                                                                                    >  1:86633/1‑158 (MQ=255)
                                                                                                                                                                                                                                     tAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTcc                                                                                                                                                                    <  2:86633/158‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAACTTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGcac                                                                                                                                                                                    <  1:911043/125‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGcac                                                                                                                                                                                    >  2:911042/1‑125 (MQ=255)
                                                                                                                                                                                                                                                        aGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATag                                                                                                                         >  1:359452/1‑182 (MQ=255)
                                                                                                                                                                                                                                                        aGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATag                                                                                                                         <  2:359452/182‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTaa                      >  1:793689/1‑281 (MQ=255)
                                                                                                                                                                                                                                                               ttGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCa                 <  1:24754/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTc                                                                                                              <  1:648943/183‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTc                                                                                                              >  2:648943/1‑183 (MQ=255)
                                                                                                                                                                                                                                                                   tATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTa           <  2:338453/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTGTCTGTGGCGAAATCGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATATAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTATCCAATCATAAAATATCTACGTGCTGCATTTTATTATTTAATATAGAGGTATAACACTGGTAAAACATAAAGAAGCATGCAATGAAGGGAAAAATAATAAGGCAATAAAAAACATAGAAAGAAACTCTTTTAATATAACATTTAAACGTAGGATTaa          >  2:162978/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                       tccgtgGCGAAATTGACTAACTTCGTTTTTTTGATTAAGAATTATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTACGGTAAAAAACAAAGAAGAATTGAATGCAGGGAAAAATAATATGGACATAAAAAACATAGAAAGAAACTCTTTTAATGAAAAATGTAAACGCAAGGTTAATca       >  2:649354/4‑280 (MQ=255)
                                                                                                                                                                                                                                                                       tCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATc        >  2:760845/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                       tCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATc        >  1:462864/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                         tgtgGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAg                                                                                                                                                                                  <  2:274661/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tgtgGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAg                                                                                                                                                                                  >  1:274661/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                         tgtgGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTc     >  1:465571/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                          gtgGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATc                                                                                                                                                                                                                                          <  1:833138/51‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gtgGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATc                                                                                                                                                                                                                                          >  2:833137/1‑51 (MQ=255)
                                                                                                                                                                                                                                                                            ggCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCata  >  1:186142/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                            ggCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCata  >  1:1079326/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                      
TGCAAGAAAGGACATGGAATCCGGAGTTAAAAACTGTTTGTTTTAGAATATTTAATTATTTAACTGCTGAATCTTCCTTGCGAGAAATATATCCATTATTACATTTCATGCCATTTTTAATATAGATTGCTCCTGGGAAAATTCCCGTATTATTACCGCAGGGGTGATAATTAGTATTGACATATCCTGTGACAAAAGGAGCTATTAAAGGTGCTATTACGATAGCTATTAGTAAAAATATAAGAGTTAGCTGTATTGTTATGTCTGTGGCGAAATTGACTACCTTCGTTTTTTTGATTAAGAATGATTTTATTATCGTAAGTAAAATTACATGAATATTTAAAAAGGAAAACGACATGAAACCGAAGCACAGAATCAACATTCTCCAATCATAAAATATTTCCGTGGAGCATTTTATTATTGAATATAGAGGTTTAACTCCGGTAAAAAACAAAGAAGCATTGAATGCAGGGAAAAATAATATGGCCATAAAAAACATCGAAAGAAACTCTTTTAATTTAACATGTAAACGCATGGTTAATCCTCATA  >  NC_000913/238202‑238750

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