| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 4,053,569 | A→C | 9.5% | intergenic (+151/+78) | 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_000913 | 4,053,569 | 0 | A | C | 9.5% | 135.8 / 5.3 | 41 | intergenic (+151/+78) | hemN/yshB | coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein | 
| Reads supporting (aligned to +/- strand): ref base A (23/14); new base C (2/2); total (25/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.37e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01 | |||||||||||
GCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCA  >  NC_000913/4053431‑4053716                                                                                                                                           |                                                                                                                                                   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                                                                                                                                              <  2:716130/85‑1 (MQ=255)                                                             tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCatgat                                                                                                                                              >  1:716130/1‑85 (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                                                       <  3:746573/107‑1 (MQ=255)                                                                                                                              aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                                                       >  4:746573/1‑107 (MQ=255)                                                                                                                              aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg             >  1:166960/1‑149 (MQ=255)                                                                                                                              aTCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg             >  1:1406359/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            >  4:1118636/1‑149 (MQ=255)                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt            >  3:1587239/1‑149 (MQ=255)                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt            >  2:1183216/1‑149 (MQ=255)                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt            >  2:1079171/1‑149 (MQ=255)                                                                                                                               tCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgt            <  1:1079171/149‑1 (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)                                                                                                                                          gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCa  >  1:1781752/1‑149 (MQ=255)                                                                                                                                           aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc                                                <  3:401537/102‑1 (MQ=255)                                                                                                                                           aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc                                                >  4:401537/1‑102 (MQ=255)                                                                                                                                           |                                                                                                                                                   GCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCA  >  NC_000913/4053431‑4053716 | 
| 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: ‑ |