Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,152,543 2 bp→GT intergenic (‑36/+112) SAUSA300_RS10965 ← / ← SAUSA300_RS10970 sucrose‑6‑phosphate hydrolase/LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,152,5430AG100.0% 27.5 / NA 13intergenic (‑36/+113)SAUSA300_RS10965/SAUSA300_RS10970sucrose‑6‑phosphate hydrolase/LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (1/12);  total (1/12)
*NC_0077932,152,5440AT92.3% 27.4 / ‑3.6 13intergenic (‑37/+112)SAUSA300_RS10965/SAUSA300_RS10970sucrose‑6‑phosphate hydrolase/LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (1/11);  total (1/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TGATAACGTTCTTCTCTAGTCCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTAATAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCATTTATAGTTATATGGTATCGGTTCCACA  >  NC_007793/2152476‑2152633
                                                                   ||                                                                                         
tGATAACGTTCTTCTCTAGTCCATTCGGTCATAATATTTACCCATTCCTCCTTTTATAATTTATGCTGTTAACCATTTATATTCTATTATAGCAAAtttt                                                            >  1:4330320/1‑100 (MQ=255)
        ttcttcTCTAGTCCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGt                                                    <  1:1745579/100‑1 (MQ=255)
        ccgatctcTAGTCCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGt                                                    <  1:1791812/96‑1 (MQ=255)
            gatctAGTCCTTTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGtggct                                                <  1:4615760/98‑5 (MQ=255)
             ctctAGTCCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTTCCTACGTTCttt                                               <  1:5002177/100‑1 (MQ=255)
             atctAGTTCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCttt                                               <  1:3024727/99‑1 (MQ=255)
                    ccgttCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACtttt                                        <  1:5700833/97‑1 (MQ=255)
                         cGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCGAATTTTACCTACGTTCTTTAACTTTTAACCt                                   <  1:2822753/100‑1 (MQ=255)
                          ggTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTa                                  <  1:5162207/100‑1 (MQ=255)
                                aataTTTACCCAGTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCAtt                            <  1:3806126/100‑1 (MQ=255)
                                            gTCCTCCTTTTATAATTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCATTTATAGTTATATg                <  1:3503849/100‑1 (MQ=255)
                                                          aTTTATGCTGTTAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCATTTATAGTTATATGGTATCGGTTCcaca  <  1:5208979/100‑1 (MQ=255)
                                                          aTTTATGCTGATAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCATTTATAGTTATATGGTATCGGTTCcagg  <  1:2288300/100‑3 (MQ=255)
                                                                   ||                                                                                         
TGATAACGTTCTTCTCTAGTCCATTCGGTCATAATATTTACCCAGTCCTCCTTTTATAATTTATGCTAATAACAATTTATATTCTATTATAGCAAATTTTACCTACGTTCTTTAACTTTTAACCTATCCATTTATAGTTATATGGTATCGGTTCCACA  >  NC_007793/2152476‑2152633

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: ‑