Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,066,175 T→C intergenic (‑411/+130) cobU ← / ← insH‑6 cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,066,1750TC96.3% 83.6 / ‑4.0 27intergenic (‑411/+130)cobU/insH‑6cobinamide‑P guanylyltransferase/cobinamide kinase/CP4‑44 prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (13/13);  total (13/14)
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.31e-01

ACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCC  >  NC_000913/2066039‑2066298
                                                                                                                                        |                                                                                                                           
aCGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCggg                                                                                                                          >  1:830103/1‑140 (MQ=255)
          tgtTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTATGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                <  1:64686/140‑1 (MQ=255)
                        ggCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                  <  2:328379/140‑1 (MQ=255)
                        ggCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                  <  2:124352/140‑1 (MQ=255)
                                       gTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                   >  2:511399/1‑140 (MQ=255)
                                          aGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCGGTCATTTCATTTCTTCATGt                                                                                >  2:618507/1‑140 (MQ=255)
                                                         cTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                 <  2:400546/140‑1 (MQ=255)
                                                          tGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtc                                                                >  2:203843/1‑140 (MQ=255)
                                                           gTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                               <  2:247471/140‑1 (MQ=255)
                                                                  aaCTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                        <  1:511399/140‑1 (MQ=255)
                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                              <  1:34359/140‑1 (MQ=255)
                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                              >  2:334886/1‑140 (MQ=255)
                                                                            acaacGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                              <  2:830103/140‑1 (MQ=255)
                                                                                 gTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                         <  2:610959/140‑1 (MQ=255)
                                                                                        aaTTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                  <  1:203843/140‑1 (MQ=255)
                                                                                            aTCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                              >  2:185780/1‑140 (MQ=255)
                                                                                             tCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                             >  1:815783/1‑140 (MQ=255)
                                                                                              ccGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGt                            <  2:815783/140‑1 (MQ=255)
                                                                                               cGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGtt                           >  2:372150/1‑140 (MQ=255)
                                                                                                  aaTGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTTAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttct                        >  1:575166/1‑140 (MQ=255)
                                                                                                    tggtggCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                      <  1:703151/140‑1 (MQ=255)
                                                                                                       tggCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                   >  1:107565/1‑140 (MQ=255)
                                                                                                       tggCGACTATGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                   >  1:465996/1‑140 (MQ=255)
                                                                                                                tGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg            >  1:623319/1‑138 (MQ=11)
                                                                                                                tGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg            <  2:623319/138‑1 (MQ=14)
                                                                                                                tGCACTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGc          >  2:140052/1‑140 (MQ=14)
                                                                                                                        ggAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc  <  1:109710/140‑1 (MQ=32)
                                                                                                                                        |                                                                                                                           
ACGGGCATTGTGTTCTGAAATCCAGGCTATTGATTCAACGTCAGCGACAGATCTGCGCTGTAAAAGAACTTGTAAAACAACGTAATAGAATTATCCGGAATGGTGGCGACTATGCACTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCC  >  NC_000913/2066039‑2066298

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