Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,053,573 Δ1 bp 8.7% intergenic (+155/+74) 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,5730T. 8.7% 157.6 / 8.9 46intergenic (+155/+74)hemN/yshBcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (26/16);  new base . (2/2);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGC  >  NC_000913/4053431‑4053721
                                                                                                                                               |                                                                                                                                                    
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                                                                                                                                                  <  2:1067624/96‑1 (MQ=255)
                                                   gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgatg                                                                                                                                                  >  1:1067624/1‑96 (MQ=255)
                                                             tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgat                                                                                                                                                   >  1:716130/1‑85 (MQ=255)
                                                             tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgat                                                                                                                                                   <  2:716130/85‑1 (MQ=255)
                                                                                        ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                               <  4:1325439/62‑1 (MQ=255)
                                                                                        ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                               >  3:1325439/1‑62 (MQ=255)
                                                                                               ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                                >  4:632184/1‑54 (MQ=255)
                                                                                               ttttGTAGGCAGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGcc                                                                                                                                                <  3:632184/54‑1 (MQ=255)
                                                                                                                              aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                                                            >  4:746573/1‑107 (MQ=255)
                                                                                                                              aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                                                            <  3:746573/107‑1 (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                                                                        <  3:1450511/94‑1 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                        >  4:1450511/1‑94 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt                 >  2:1079171/1‑149 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt                 <  1:1079171/149‑1 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt                 >  4:1118636/1‑149 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt                 >  2:1183216/1‑149 (MQ=255)
                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt                 >  3:1587239/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                                                              >  4:719158/1‑100 (MQ=255)
                                                                                                                                   gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                                                              <  3:719158/100‑1 (MQ=255)
                                                                                                                                     aTC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCa                                          <  1:352206/118‑1 (MQ=255)
                                                                                                                                     aTC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCa                                          >  2:352206/1‑118 (MQ=255)
                                                                                                                                      tC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTg          >  4:1036462/1‑149 (MQ=255)
                                                                                                                                          gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCa       >  1:1781752/1‑149 (MQ=255)
                                                                                                                                           aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc                                                     >  4:401537/1‑102 (MQ=255)
                                                                                                                                           aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc                                                     <  3:401537/102‑1 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   >  4:803811/1‑135 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   <  3:803811/135‑1 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAc     >  3:1208173/1‑149 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAc     <  4:1587239/149‑1 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACCGCAGTTTGCGGTGTGTGGCTGTCAAc     <  3:1176546/149‑1 (MQ=255)
                                                                                                                                            gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGACAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAc     <  3:1814685/149‑1 (MQ=255)
                                                                                                                                             cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACt    >  1:1732058/1‑149 (MQ=255)
                                                                                                                                             cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACt    <  2:166960/149‑1 (MQ=255)
                                                                                                                                             cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACt    >  1:1764347/1‑149 (MQ=255)
                                                                                                                                             cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACt    >  4:1851773/1‑149 (MQ=255)
                                                                                                                                               tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   <  2:1726385/132‑1 (MQ=255)
                                                                                                                                               tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   >  1:1726385/1‑132 (MQ=255)
                                                                                                                                               tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   >  1:1081208/1‑132 (MQ=255)
                                                                                                                                               tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgt                   <  2:1081208/132‑1 (MQ=255)
                                                                                                                                               tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTgc  <  3:195932/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAACTGC  >  NC_000913/4053431‑4053721

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