Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 903,337 C→T intergenic (+44/‑146) SAUSA300_RS04430 → / → SAUSA300_RS04435 Fe‑S cluster assembly protein SufB/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793903,3370CT93.3% 33.9 / ‑4.3 15intergenic (+44/‑146)SAUSA300_RS04430/SAUSA300_RS04435Fe‑S cluster assembly protein SufB/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/0);  total (14/1)
Fisher's exact test for biased strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CGGTCGAAATGAACCGTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTA  >  NC_007793/903241‑903433
                                                                                                |                                                                                                
cgGTCGAAATGAACCGTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTAt                                                                                               >  1:4975881/1‑100 (MQ=255)
          gAACCGTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCAtt                                                                                     >  1:4082568/1‑100 (MQ=255)
          agaCCGTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCAtt                                                                                     >  1:3174432/3‑100 (MQ=255)
            tgggTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCAtttt                                                                                   >  1:3556371/4‑100 (MQ=255)
                  tAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAcc                                                                             >  1:1218264/1‑99 (MQ=255)
                    tgcAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccat                                                                           >  1:1443175/3‑97 (MQ=255)
                    aTCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccat                                                                           >  1:1724857/1‑97 (MQ=255)
                     acAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccatt                                                                          >  1:399621/2‑96 (MQ=255)
                      cAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccattt                                                                         >  1:3761859/1‑95 (MQ=255)
                      cAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccattt                                                                         >  1:38147/1‑95 (MQ=255)
                      cAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccattt                                                                         >  1:5878064/1‑95 (MQ=255)
                       caGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccattta                                                                        >  1:1486312/2‑94 (MQ=255)
                        aGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccatttac                                                                       >  1:1585063/1‑93 (MQ=255)
                        aGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCTATTCTTAACTGTCGAAAACGGTTATTGATACCATTTTGATAccatttac                                                                       >  1:5115941/1‑93 (MQ=255)
                                                                                             tGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTa  <  1:1968248/100‑1 (MQ=255)
                                                                                                |                                                                                                
CGGTCGAAATGAACCGTTTAATCAAGTTCGAAATGGAAGGTAGTATTGGATAGCTTTAAACCGCGTTGTTAAGCCATTCTTGACTTCCGAAAATGGCTATTGATACCATTTTGATACTGAATATAACAAAAAGCCACATTACTGTGGCTTTTTTTGTTTTATAACTAAATCGGATTGATAGATAAGCTTTGTA  >  NC_007793/903241‑903433

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