Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,368,101 T→C 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,1010TC94.1% 76.6 / NA 34A21A (GCT→GCCarnAfused UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (15/17);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTG  >  NC_000913/2367965‑2368236
                                                                                                                                        |                                                                                                                                       
tACACCCATGTCCGCGCCTGCCCCCGCTATTTTGTCCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGCAACCCGTCGTTTGTGCCTACCACGATCTGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGt                                                                                                                                   <  2:691704/143‑1 (MQ=255)
 acacCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAGACCGTCGTTTTTGCCTACCACGCTCTGGGATGCCTCGGTAGTGAAGGCCTGCGGGCCGCCGGtt                                                                                                                                  >  1:905618/1‑143 (MQ=255)
  cacCGGTGGGCGCGCCGGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTa                                                                                                                                 <  2:270067/143‑1 (MQ=255)
   acCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCTGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTAc                                                                                                                                <  1:706625/143‑1 (MQ=255)
     cGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGtctctc                                                                                                                              >  1:168828/1‑138 (MQ=255)
       ctGTCCGCGTCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGaaa                                                                                                                            <  2:940028/142‑1 (MQ=255)
              cgcccgccCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGCTGCCTCGGTATTGCAGCCCTGCTGGCCGCCGGTTACGCAATCAgcgc                                                                                                                     >  1:958758/1‑143 (MQ=255)
                   gcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTACTGGCCGCCGGTTACGAAATTAGCGCCAttt                                                                                                                >  2:932988/1‑143 (MQ=255)
                                 gTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGTCCGCCGGTTTCGAAATTAGTGCCATTTTTACCAATACTGAt                                                                                                  >  1:465149/1‑143 (MQ=255)
                                                cGTCCATCCAGCAAGGAAAATGAATAATGAAGACCGGCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                                                                   <  2:467554/143‑1 (MQ=255)
                                                      cccaGCAAGGAAATTGAATAATGAAAACCGTCGTTTTTGCCTCCCACGATATGGGATGCCTCGGTATTGAAGTCTTGTTGGTCGTCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCttt                                                                             <  1:673257/142‑1 (MQ=255)
                                                        caGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCACTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATGCCGGTGAAAAAGGCTtata                                                                           >  2:921280/1‑143 (MQ=255)
                                                         aGCAAGGAAGATGAAGAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTAt                                                                          <  2:685677/143‑1 (MQ=255)
                                                               aaaaaTAAATAATGAAAACCGTCCTTTTTCCCTACCACGATATGGTATGCCTCGGTATTGAAGCCCTTCTGGCTTCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCg                                                                    <  1:361520/142‑1 (MQ=255)
                                                                 aaaTGAATTACGACAAGCATCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGt                                                                  <  2:391190/143‑1 (MQ=255)
                                                                     gAATAATGAAAACCGTCTTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGTTGTTTGTTATGAAATTAGAGCCATTTTTACACAAACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                                                              <  1:85951/143‑1 (MQ=255)
                                                                             tcgaCCGGCGTTGGTGCGTACGACGATGTGGGATTCGTCGGGATTGGAGCGCTTCTGGGCGCTGGTTTTGACATTGGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc                                                      <  2:299768/140‑1 (MQ=255)
                                                                             gtaaCCGTCGTTTTTGCCTACTCCGATATGGCATGCCTTGGTCTTGAAGCCCTGCTGGCCGCCGGTTACGACAGTAGCGCCATGTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc                                                      <  2:661218/141‑1 (MQ=255)
                                                                             gcatccGTCGTTTTTGCCTGCCACGATCTGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAACTTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc                                                      <  2:898939/139‑1 (MQ=255)
                                                                                       tttttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATCCTGCTAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGacgga                                            >  1:197900/1‑139 (MQ=255)
                                                                                           tGCCGACGACGATCGTGGATGCCTCTGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCGTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCa                                        <  2:149595/143‑1 (MQ=255)
                                                                                              cTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATATTTTTGGTGAAAGAGCCATTAATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTc                                     <  1:783608/143‑1 (MQ=255)
                                                                                                  gtgcagctGGGACGCCTCGGTATCGACCCGCTGCTGGCCGCCGGCTACGAAATTAGCGCCATTTTTACCCATACTGATACTCCCGGTGCAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGt                                 <  1:826137/136‑1 (MQ=255)
                                                                                                  cACGATATGGGATGCCTCGGAATTGAAGCCCTGCAGGCCGCCGGTTACGAAAATAACGACATATTTACGCATAGTGATAATCCAGGTGAAAAAGCCAGTGATGGGACGGTGGCTCGAGTGGCAGCAGACAAAGACAGTGCGGt                                 >  2:51061/1‑143 (MQ=255)
                                                                                                   aCGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGCCTGGCGGCGGCACGAGGCCTTCCGCtt                                >  1:359448/1‑143 (MQ=255)
                                                                                                                  tCGGTATTGAAGCCCTGCTGGCCGCCGATTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAAc                 >  1:392568/1‑143 (MQ=255)
                                                                                                                   cgctATTGAAGGCCTGCTGGCCGCCGGTTACGAGATGTGCGTCATTTTTACCCATACTGGTAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCTGTTTATGCGCCGGATAACg                <  2:140488/140‑1 (MQ=255)
                                                                                                                    ggTATTGAAGCCCTGCTGGCCGCCGGGTACGAAATTAGCGCCATTGTTACCCATACTGATAATCCCGGTGAAAGAGCCTTGTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGACATTCCGGTTTATGCGGCGGagagggt               >  2:980604/1‑137 (MQ=255)
                                                                                                                        ttGAAGCCCTGCTGGCCGCCGGTTACGACATTAGCGCCATTTTTTCCCATTCTGATAATCCCGGTGATAATGTCTTTTATTGTTCTGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAat           <  1:598383/143‑1 (MQ=255)
                                                                                                                        ttGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGGTTATGCGCCGGATAACGTTAat           >  1:151985/1‑143 (MQ=255)
                                                                                                                           aaGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAGGAGGCATTCCGGTTTATGCGCCGGATAGCGTTAatcat        >  2:699769/1‑143 (MQ=255)
                                                                                                                            aGCCCTGCAGGCCGACGGTTACGAAATTAGCGCCATTTATACCCATACTGATAATCCCGGTAAAAAAGCTTTTTATGGTTTGGTGGCTAGTCTGGCGGCGAAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAatcatc       <  1:53841/143‑1 (MQ=255)
                                                                                                                                ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCt   >  1:334768/1‑143 (MQ=255)
                                                                                                                                 tgctgGCCGCCGGTTACGAAATTAGCGCCCATCATACCCATACTGATAATCCCGGTGAAAAAGCCTGTTATGGGTCGGTGGCTCGTCTGGCGGCGGAAAGAGGACTACCGGTTTACGCGCCGCCGAACGTGACTCATcccgtg  >  2:819126/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTG  >  NC_000913/2367965‑2368236

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: