Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,368,101 T→C 100% A21A (GCT→GCC arnA → fused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,1010TC97.1% 66.0 / ‑6.0 26A21A (GCT→GCCarnAfused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (16/8);  minor base A (0/1);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGG  >  NC_000913/2367963‑2368239
                                                                                                                                          |                                                                                                                                          
tCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGACTAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGccgccg                                                                                                                                        >  1:535574/1‑143 (MQ=255)
              ggtgcccccccCCGCTCTTTTGTTCAGCAAGTTATCCGTCCCTCCAGCATGGAAAATGAATAATTAAAACCGTCGTTTTTGCCTATCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAgc                                                                                                                          <  1:480773/140‑1 (MQ=255)
                               ttttGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATAc                                                                                                         <  2:331573/143‑1 (MQ=255)
                                 ttctttcgCAACTTATCCTGCCATCAACGAAGGAGACTGACTAATGACAACCGTTGTTTTTGCCTCCCACCATTAGCGTTGCCTCGTTATTGAAGCCCTGCTTGCTGCCGTTTATGAATTTAGCGCCATTTTTACCCATACTg                                                                                                       <  2:658712/136‑1 (MQ=255)
                                      cAGCAAGTTATCCGTCCATCCAGCAAGGAAACTGAATAATGAAAACCGTCGTTTTTGCCTACCACGCTCTGGGCTGCCTCGGTATCGAAGCCCTGCTGGCCGCCGGTTACGAAAATAGCGACATTTTTACAAATACTGaaaaa                                                                                                  >  1:191702/1‑142 (MQ=255)
                                        gCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGGTACGCAGTTAGCGCCATTTTTACCCATACTGATAATcc                                                                                                >  1:490336/1‑143 (MQ=255)
                                                   gTCCATCCAGCAAGGAAAATGAATGATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGCAGCCCTGCTGGCCGCCGCTTACGAACTTAGCGCCATTTTTACCCATACTGATAAACACGGTaaaaaaa                                                                                     >  1:13860/1‑142 (MQ=255)
                                                                ggAAAATGAAGAATGAAAACCGTCGTTTTTGCCTACCTCGATATGGGATGCCTCTGTATTGAAGCCTTGTTGGTCGTCGGTTACGAAATTAGCGCCAATTATAGCCATACTGATAAGCCCGGTGAAAAAGCCTTTTATGGTTc                                                                        <  1:69883/143‑1 (MQ=255)
                                                                      tGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTCCGAAATTAGCGCTATTTTTACACATAATGATAATCACGGTGAAAAAGCGTTTTATGGTTCGGTGGc                                                                  >  1:388183/1‑143 (MQ=255)
                                                                                 aaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTTCTCTTACTGATAATGCCGGTGAAAAAGCCTGTTATGGTTCGGTGGCTCGTCTggcgg                                                       <  1:557041/143‑1 (MQ=255)
                                                                                      tcgtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGaa                                                  >  2:451227/1‑143 (MQ=255)
                                                                                         gtttgGCCTACCACGATATCGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGATATTATCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTGGATGGTTCGGTGGCTCGTCTGGCGGCGGAAaga                                               <  1:130984/142‑1 (MQ=255)
                                                                                             tGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCATTGTGTGGTTCGGTGGCTCGTCTGGCGGCGGAAACAGGCa                                           >  1:191985/1‑143 (MQ=255)
                                                                                                 tACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTcc                                       >  2:392611/1‑143 (MQ=255)
                                                                                                 tACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGCCGGGCAGAGGCGTTcc                                       >  1:82145/1‑143 (MQ=255)
                                                                                                     aCGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGCAAGAGGCATTCCGGtt                                   >  1:380856/1‑143 (MQ=255)
                                                                                                        ttatGGGATGCCTCGGTATTGAAGCCCTGCTTGCCGCCGGTTTTGTAGTTAGCCCCATTTTTTTTCTTACTAATAAACCCGGTGAAAAAGCCTTTTATGGATCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTAt                                <  1:122789/142‑1 (MQ=255)
                                                                                                                tGCCTCGGTATTGAAGCCCGGCTATCCTCTGATTACGATAGTAGTGGCATTATCACCCACACTGAGAATCGCTGTTGAAAAGCCTTTGATGGGGCGGTTGGTTGTCTGTCGTCGGACAGGGGCATGGGGGTTTATGCGCCGGa                        >  2:435682/1‑143 (MQ=255)
                                                                                                                   cTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATac                     >  2:380961/1‑142 (MQ=255)
                                                                                                                   cTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATaa                     >  1:275886/1‑143 (MQ=255)
                                                                                                                    tCGGTATTGAGGCCCTGCTGGCCGCGGGTTACGAAATTAGCGCCATTTTTACCCAGGCTGAGAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGTCATTCCGGTTTATGCGCCGGATAAc                    <  2:527313/143‑1 (MQ=255)
                                                                                                                    tCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAAc                    >  1:64076/1‑143 (MQ=255)
                                                                                                                       gTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTCTGGTTCGGTGGCTCGTCTGGCGGCGGAAACAGGCATTCCGGTTTATGCGCCGGATAACGtt                 >  1:64037/1‑143 (MQ=255)
                                                                                                                         aTTGAATCCCTCCGCGCATCCGGTGACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCTGTGAAAAAGCCTTTTAGGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTaa               <  2:490336/143‑1 (MQ=255)
                                                                                                                             aaGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTCTTATGGTTCGGTGGCTCGTCTGCCTGCGAAAAGAGGCATTCCAGTTTATGCGGCGGCACACGATAatcat           >  2:44574/1‑143 (MQ=255)
                                                                                                                                cccTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGACGTAAAGACGCATTCCGGGTTATGCGCCGGATAACGTTAATCATCCg        >  1:148825/1‑143 (MQ=255)
                                                                                                                                 ccTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGTTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGc       <  1:523291/143‑1 (MQ=255)
                                                                                                                                  ctcccggccgCCGATTACGAACTTAACGCCATTGTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCCGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCt      <  2:611210/138‑1 (MQ=255)
                                                                                                                                    cctggccgCCGGTTACGAAATGCGCGCCATTTTGACCCATACTGATCATCCCGCCGAGAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCtgt    <  2:458854/142‑1 (MQ=255)
                                                                                                                                      tggccgCCGGTTACGAAATTAGCGCCATTTGTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTgtgg  >  1:232738/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGG  >  NC_000913/2367963‑2368239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: