Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 774,306 C→T T101M (ACG→ATG)  SAUSA300_RS03755 → aminodeoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793774,3060CT94.1% 84.9 / ‑2.2 36T101M (ACG→ATG) SAUSA300_RS03755aminodeoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (28/4);  minor base A (2/0);  total (32/4)
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.43e-01

AACGATAGTGTCACCTTTACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAACGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCAGAGGACGG  >  NC_007793/774211‑774372
                                                                                               |                                                                  
aaCGATAGTGTCACCTTTACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACa                                                                >  1:79163/1‑100 (MQ=255)
           cACCTTTACCTGTGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCTAACTTTTATAAAAAATAAAATGACAGATCGGAAgag                                                     >  1:5291697/1‑98 (MQ=255)
           cACCTTTACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATAACAGAACGAAAGCa                                                     >  1:4836576/1‑100 (MQ=255)
           cACCTTTACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAACGACAGagatcggaag                                                     >  1:6156451/1‑91 (MQ=255)
               tttACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTa                                                 >  1:5590420/1‑100 (MQ=255)
                ttACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAGATGACAGATCGgaaaagcacac                                                >  1:4000825/1‑92 (MQ=255)
                ttACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAg                                                >  1:2330039/1‑100 (MQ=255)
                ttACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAg                                                >  1:3224693/1‑100 (MQ=255)
                 tACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGAAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGc                                               >  1:6145881/1‑100 (MQ=255)
                 gaCCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAACCCAACTTTTATAAAAAATAAAGTGACAGATCGGAAgagcacacg                                               >  1:2085358/2‑92 (MQ=255)
                  aCCTGAGAAAAGTAGTTTTACAGTAAAAGGTGAAGTGTTGCCCAAAGAAGATAATCCAACTAGAAGAAAAAAAAAAGTGACAGATCGTAAgagcacacgt                                              >  1:1969137/1‑91 (MQ=255)
                  aCCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAATAGACAGATCGGAAgagcacacgt                                              >  1:4875212/1‑91 (MQ=255)
                  aCCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAAAGACAGATCGGAAgagcacacgt                                              >  1:3681579/1‑91 (MQ=255)
                      gagaAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGcacaca                                          >  1:5138271/1‑100 (MQ=255)
                                gtTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATTAAAc                                <  1:265519/100‑1 (MQ=255)
                                gtTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAAc                                <  1:2101613/100‑1 (MQ=255)
                                gtTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCAGGAAAc                                <  1:711485/100‑1 (MQ=255)
                                gtTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCACGAAAc                                <  1:4330326/100‑1 (MQ=255)
                                    tACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGAt                            >  1:4600242/1‑100 (MQ=255)
                                     aCAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGAtt                           >  1:5783479/1‑100 (MQ=255)
                                     aCAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCAAGAAACAGAtc                           >  1:2121513/1‑99 (MQ=255)
                                       aGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCAGGAAACAGAtcgg                         >  1:1915960/1‑97 (MQ=255)
                                                      cttGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCa          >  1:4593230/2‑100 (MQ=255)
                                                       ttGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCag         >  1:3682714/1‑100 (MQ=255)
                                                        tGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAACGGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCaga        >  1:5426264/1‑100 (MQ=255)
                                                        tGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCaga        >  1:3016840/1‑100 (MQ=255)
                                                        tGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCaga        >  1:5887488/1‑100 (MQ=255)
                                                        tGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCaga        >  1:5379639/1‑100 (MQ=255)
                                                        tGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCaga        >  1:5270191/1‑100 (MQ=255)
                                                         gCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCagag       >  1:1330524/1‑100 (MQ=255)
                                                         gCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCagag       >  1:5282653/1‑100 (MQ=255)
                                                         gCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCagag       >  1:594603/1‑100 (MQ=255)
                                                         gCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCagag       >  1:178829/1‑100 (MQ=255)
                                                         gCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCagag       >  1:1612065/1‑100 (MQ=255)
                                                          cccAAAGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCAGAgg      >  1:6375569/1‑100 (MQ=255)
                                                              aaGTAGTTAATCCAACTTTTATAAAAAATAAAATGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCAGAGGACgg  >  1:4128653/1‑100 (MQ=255)
                                                                                               |                                                                  
AACGATAGTGTCACCTTTACCTGAGAAAAGTAGTTTTACAGCAAAATTTCAAGTGTTGCCCAAAGTAGTTAATCCAACTTTTATAAAAAATAAAACGACAGAACGAAAGCATTTAGCACACAATCATGAAACAGATTTAATATTGCTAACTTCAGAGGACGG  >  NC_007793/774211‑774372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑