Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,053,566:1 +A 10.0% intergenic (+148/+81) hemN → / ← yshB coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,053,5661.A10.0% 160.1 / 14.0 40intergenic (+148/+81)hemN/yshBcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (22/14);  new base A (2/2);  total (24/16)
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.78e-01

TGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTG  >  NC_000913/4053427‑4053713
                                                                                                                                            |                                                                                                                                                   
tGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑t                                                                                                                                                    >  1:1021523/1‑141 (MQ=255)
tGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑t                                                                                                                                                    <  2:1021523/141‑1 (MQ=255)
    gcTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                        <  1:932353/149‑1 (MQ=255)
     cTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatg                                                                                                                                             <  4:987190/143‑1 (MQ=255)
     cTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatg                                                                                                                                             >  3:987190/1‑143 (MQ=255)
               aaaGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgatgtc                                                                                                                                        <  2:1490871/138‑3 (MQ=255)
               aaaGTAGAAAACGAAAGCGCCATCAATGTCCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                        >  1:1490871/1‑138 (MQ=255)
                           gAAAGCGCCATCAATGTGCTGAGCACAAGATTGACTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCATGCTCA‑GATGCCTGATgcg                                                                                                                                 <  4:1183486/133‑1 (MQ=37)
                           gAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAGATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCATGCTCA‑GATGCCTGATgcg                                                                                                                                 >  3:1183486/1‑133 (MQ=38)
                                                   aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGc                                                                                                                                         >  4:1403010/1‑101 (MQ=255)
                                                   aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGc                                                                                                                                         <  3:1403010/101‑1 (MQ=255)
                                                    cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCATGCTCA‑GATGCCTGAt                                                                                                                                    <  2:90712/105‑1 (MQ=37)
                                                    cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTGGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCATGCTCA‑GATGCCTGAt                                                                                                                                    >  1:90712/1‑105 (MQ=25)
                                                       gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgatg                                                                                                                                          >  1:1067624/1‑96 (MQ=255)
                                                       gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgatg                                                                                                                                          <  2:1067624/96‑1 (MQ=255)
                                                                 tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgat                                                                                                                                           <  2:716130/85‑1 (MQ=255)
                                                                 tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgat                                                                                                                                           >  1:716130/1‑85 (MQ=255)
                                                                                            ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                       >  3:1325439/1‑62 (MQ=255)
                                                                                            ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                       <  4:1325439/62‑1 (MQ=255)
                                                                                                   ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                        >  4:632184/1‑54 (MQ=255)
                                                                                                   ttttGTAGGCAGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                        <  3:632184/54‑1 (MQ=255)
                                                                                                                                  aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                                                    <  3:746573/107‑1 (MQ=255)
                                                                                                                                  aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                                                    >  4:746573/1‑107 (MQ=255)
                                                                                                                                  aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg          >  1:1406359/1‑149 (MQ=255)
                                                                                                                                  aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg          >  1:166960/1‑149 (MQ=255)
                                                                                                                                   tcgggCATC‑TGAGCATGATGCTTGATGCGACGCTGGTCGCGTCTTATCAGGCCTAGGACTTTCTCCGCATCCGGCATAAGCAACGTGCTGCAGCTAAACAGCCCAGTCATTGTGGCACAGAGAGCGGCAGCAACGGCAGTTTGCGgtgt         >  4:1176546/4‑149 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                >  4:1450511/1‑94 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                <  3:1450511/94‑1 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt         >  2:1183216/1‑149 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt         >  4:1118636/1‑149 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt         <  1:1079171/149‑1 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt         >  3:1587239/1‑149 (MQ=255)
                                                                                                                                   tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt         >  2:1079171/1‑149 (MQ=255)
                                                                                                                                     aGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                                                                    >  3:1730335/1‑88 (MQ=255)
                                                                                                                                     aGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                                                                    <  4:1730335/88‑1 (MQ=255)
                                                                                                                                     aGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgtgt       >  4:195932/1‑149 (MQ=255)
                                                                                                                                      ggCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGACAGCAACGGCAGTTTGCGGtgtgtg      >  4:1814685/1‑149 (MQ=255)
                                                                                                                                       gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                                                      <  3:719158/100‑1 (MQ=255)
                                                                                                                                       gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                                                      >  4:719158/1‑100 (MQ=255)
                                                                                                                                         aTC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCa                                  >  2:352206/1‑118 (MQ=255)
                                                                                                                                         aTC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCa                                  <  1:352206/118‑1 (MQ=255)
                                                                                                                                          tC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTg  >  4:1036462/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTG  >  NC_000913/4053427‑4053713

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