Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,180,739 (G)6→7 coding (529/789 nt) ycfZ ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,180,7331.G97.3% 128.4 / ‑3.2 37coding (535/789 nt)ycfZputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (17/19);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CCATATTGGCAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCC‑‑GGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGATCTGCAT  >  NC_000913/1180604‑1180861
                                                                                                                                  |                                                                                                                                   
cgatatTGGCTATAGCCATCGCCCCTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCC‑‑ggg                                                                                                                                 <  2:301359/131‑1 (MQ=255)
    atTGGCAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTAGTGCTGGC‑‑ggggggg                                                                                                                            >  1:184083/1‑132 (MQ=255)
         cAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACGTAATGTGTTGGCGCAGCC‑‑GGGCGGGAGCAt                                                                                                                       >  1:301359/1‑133 (MQ=255)
         cAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCAt                                                                                                                       >  2:530499/1‑133 (MQ=255)
         cAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCAt                                                                                                                       >  1:464962/1‑133 (MQ=255)
          aaTAGTCATCGCCACTGCAACGGTTAACGCATAAATAAAACTAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCAtt                                                                                                                      >  2:157919/1‑133 (MQ=255)
                  tCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCt                                                                                                              >  2:5661/1‑133 (MQ=255)
                  tCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCt                                                                                                              >  1:364763/1‑133 (MQ=255)
                  tCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCt                                                                                                              >  1:467984/1‑133 (MQ=255)
                   cGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCtg                                                                                                             <  1:5661/133‑1 (MQ=255)
                       aCTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCtgttgt                                                                                                         >  1:57510/1‑133 (MQ=255)
                         tGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGGAGCATTTTCATGCTGTTGTc                                                                                                        >  1:437294/1‑133 (MQ=255)
                         tGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTcc                                                                                                       <  2:612637/133‑1 (MQ=255)
                         tGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTcc                                                                                                       >  2:561837/1‑133 (MQ=255)
                         tGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTcc                                                                                                       >  1:384109/1‑133 (MQ=255)
                               ttttAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATg                                                                                                 <  2:485953/133‑1 (MQ=255)
                                    aCGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGa                                                                                             <  2:437294/133‑1 (MQ=255)
                                     cGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGata                                                                                           <  1:511616/133‑1 (MQ=255)
                                      gCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGataa                                                                                          <  1:336371/133‑1 (MQ=255)
                                                          gtagtaAAAGCTCCTCGAAAATAATATACACGTCGCAGCCCTGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGa                                                                      <  2:184083/133‑1 (MQ=255)
                                                                      ccTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGTTGCGGATAGTACGTTAGTGGCGGCGGCGATTGTTTTAatat                                                          >  2:217683/1‑133 (MQ=255)
                                                                        tGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCt                                                        >  1:447341/1‑133 (MQ=255)
                                                                        tGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCt                                                        <  2:57510/133‑1 (MQ=255)
                                                                          gAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGTATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCtgt                                                      <  1:561837/133‑1 (MQ=255)
                                                                              aTAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCtt                                                  <  2:137212/133‑1 (MQ=255)
                                                                                   aTAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGa                                             >  2:571539/1‑133 (MQ=255)
                                                                                    tAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGCCGTCGATTGTTTTAATATCTGCTCTTGTAGaa                                            >  2:275563/1‑133 (MQ=255)
                                                                                    tAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGGGGCGGGGTGGATTGTTTTAATATCTGTGCTTGTAGaa                                            >  2:601881/1‑133 (MQ=255)
                                                                                        cGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGc                                        <  2:364763/133‑1 (MQ=255)
                                                                                             cAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTGTTAATATATGTACTTGTAGAACAGCCAATg                                   >  1:242466/1‑133 (MQ=255)
                                                                                                     gCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTg                           <  1:571539/133‑1 (MQ=255)
                                                                                                     gCGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTg                           >  1:435904/1‑133 (MQ=255)
                                                                                                      cGATAAAGCCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGa                          <  1:157919/133‑1 (MQ=255)
                                                                                                      cGATAAAACCTAATGTGTTGGTTCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGa                          <  1:275563/133‑1 (MQ=255)
                                                                                                      cGATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGa                          <  2:42614/133‑1 (MQ=255)
                                                                                                       gATAAAACCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGAt                         <  2:464962/133‑1 (MQ=255)
                                                                                                          aaaaCCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGGTTGTAGATAAGCCAATGCATTGTTGATTcc                      >  1:195969/1‑133 (MQ=255)
                                                                                                            aaCCTAATGTGTTGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCt                    <  2:467984/133‑1 (MQ=255)
                                                                                                                        tGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGAt        <  2:435904/133‑1 (MQ=255)
                                                                                                                        tGGTGCTGCCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGAt        <  1:530499/133‑1 (MQ=255)
                                                                                                                             ctgcCGGGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGATCTGCa   <  1:601881/133‑1 (MQ=255)
                                                                                                                              tgcC‑‑GGGGGGTAGCATTTTCATGCGGTTTTCCCGGATGCGGGTAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGATCTGCAt  <  2:195969/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
CCATATTGGCAATAGTCATCGCCACTGCAAATTTTAACGCATAAATAAAACGAGTCGTGTAGTAAAAGATCCTGGAAAATAAGATAAACGTCGCAGCCCAGGCGATAAAACCTAATGTGTTGGTGCTGCC‑‑GGGGGGAGCATTTTCATGCTGTTGTCCTGGATGCGGATAATACGTTAGTGGCGGCGTCGATTGTTTTAATATCTGTGCTTGTAGAACAGCCAATGCATTATTGATTCCCTGATACCAGAGATCTGCAT  >  NC_000913/1180604‑1180861

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