Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 754,285 (A)6→7 intergenic (‑77/‑132) SAUSA300_RS03645 ← / → SAUSA300_RS03650 hypothetical protein/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793754,2791.A81.2% 33.4 / 3.4 16intergenic (‑71/‑138)SAUSA300_RS03645/SAUSA300_RS03650hypothetical protein/MFS transporter
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (10/3);  total (12/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 = 6.79e-01

GATTACTTTTTTGCAATTTCATATGATCAATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAG‑AAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTGGAAA  >  NC_007793/754187‑754367
                                                                                             |                                                                                        
gATTACTTTTTTGCAATTTCATATGATCAATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAG‑aaaaaaa                                                                                   >  1:4505270/1‑99 (MQ=255)
     cTTTTTTGCACTTTCATATGATCAATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTg                                                                               >  1:6569968/1‑100 (MQ=255)
     cTTTTTTGCAATTTCATATGATCAATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTg                                                                               >  1:5383076/1‑100 (MQ=255)
                         aTCTATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTa                                                           <  1:6240647/100‑1 (MQ=255)
                            taTCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATa                                                        >  1:5404862/2‑100 (MQ=255)
                             aTCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAG‑AAAAAATATTGTCAAACGTAGCAATGTTGTAATACa                                                      >  1:5642219/1‑100 (MQ=255)
                              tCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACa                                                      >  1:6611520/1‑100 (MQ=255)
                               ccccTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAG‑AAAAAATATTGTCAAATGTAGCAATGTTGTAATACAat                                                    >  1:4726054/1‑100 (MQ=255)
                                 ccTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAata                                                   >  1:2779040/1‑100 (MQ=255)
                                  cTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAatat                                                  <  1:5159045/100‑1 (MQ=255)
                                                         tataCAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAg                           >  1:6229421/1‑100 (MQ=255)
                                                          ataCAATTAAATGGAAAATAGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGc                          >  1:1795532/1‑100 (MQ=255)
                                                                             aGTGATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGt       >  1:824979/1‑100 (MQ=255)
                                                                              gTGATAATTACAAAG‑AAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTgg     <  1:831035/100‑1 (MQ=255)
                                                                                gATAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTGGa    >  1:4037639/1‑100 (MQ=255)
                                                                                 ctAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTGGaa   >  1:750864/2‑100 (MQ=255)
                                                                                  tAATTACAAAGAAAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTGGaaa  <  1:3887235/100‑1 (MQ=255)
                                                                                             |                                                                                        
GATTACTTTTTTGCAATTTCATATGATCAATCCCCTTTATTTTAATATGTCATTAATTATACAATTAAATGGAAAATAGTGATAATTACAAAG‑AAAAAATATTGTCAAATGTAGCAATGTTGTAATACAATATAGAAACTTTTTACGAATATTTAGCATGAATTGCAATCTGTCGTGGAAA  >  NC_007793/754187‑754367

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