Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 979,291 G→A A366A (GCG→GCA)  SAUSA300_RS04805 → peptide ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793979,2910GA97.1% 89.9 / ‑5.1 35A366A (GCG→GCA) SAUSA300_RS04805peptide ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/9);  total (26/9)
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.98e-01

ACACTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCGCACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAATATATC  >  NC_007793/979193‑979390
                                                                                                  |                                                                                                   
acacTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAtcgg                                                                                                    >  1:1470020/1‑99 (MQ=255)
acacTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAtcgg                                                                                                    >  1:1608224/1‑99 (MQ=255)
acacTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAtcgg                                                                                                    >  1:4919018/1‑99 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAtcggaa                                                                                                  >  1:5059046/1‑97 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCacac                                                                                                  >  1:3920171/1‑100 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCacac                                                                                                  >  1:3582637/1‑100 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCacac                                                                                                  >  1:3430535/1‑100 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCacac                                                                                                  >  1:1017566/1‑100 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCGCAc                                                                                                  >  1:5783250/1‑100 (MQ=255)
  acTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGACTATGCGAGTACCATAAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCacac                                                                                                  >  1:3981932/1‑100 (MQ=255)
   atatcAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACt                                                                                                 >  1:5349799/5‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTTAAGATCGAAGAGCacacgtc                                                                                               >  1:3047997/1‑97 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:3374648/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:922951/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:463696/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:358076/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:1314465/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:1609771/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:2068094/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:3010514/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:2087353/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:2508211/1‑100 (MQ=255)
     aaCAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCCTTAAAATATAATCCTAAAGAAGCAAGAGCACACTgg                                                                                               >  1:5390543/1‑100 (MQ=255)
        aGCCAAAGGAATTGCGAAAGCGCCTAATGACAAAGATTATGCAAGTACCATGAATTCGCCTTTAAAAGATAATCCTAAAGAAGCAAGAGCAAACTGGGAc                                                                                            <  1:2501203/100‑1 (MQ=255)
         gCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACa                                                                                           <  1:4170044/100‑1 (MQ=255)
                  aTTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTaaaa                                                                                  <  1:3402634/100‑1 (MQ=255)
                   ttGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTaaaaa                                                                                 <  1:5726701/100‑1 (MQ=255)
                         aaGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGtt                                                                           <  1:2497098/100‑1 (MQ=255)
                           cgatcTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTaa                                                                         <  1:2224810/96‑2 (MQ=255)
                                   tGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAAt                                                                 >  1:1866825/1‑100 (MQ=255)
                                    ggCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATg                                                                >  1:3124870/1‑100 (MQ=255)
                                             tATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACAt                                                       >  1:2395121/1‑100 (MQ=255)
                                                 cGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCTAGAGCACACTGGGACAAAGCTAATAAAGAGTTAGGTAAAATTGAAGTAACAtttta                                                   >  1:3251001/1‑99 (MQ=255)
                                                    gTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAAt                                                >  1:193337/1‑100 (MQ=255)
                                                      aCCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGTATACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGa                                              >  1:5414695/1‑100 (MQ=255)
                                                      aCCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGa                                              >  1:5658536/1‑100 (MQ=255)
                                                                      tAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACAcc                              <  1:4789740/100‑1 (MQ=255)
                                                                      tAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACAcc                              <  1:5518902/100‑1 (MQ=255)
                                                                      tAAAATATAATCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGACACAcc                              <  1:2015499/100‑1 (MQ=255)
                                                                              aaTCCTAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGAtcgga                      >  1:766222/1‑96 (MQ=255)
                                                                                   tAAAGAAGCAAGAGCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGAtcggaagagc                 >  1:3198258/1‑91 (MQ=255)
                                                                                             cgagCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAat       >  1:3289986/2‑100 (MQ=255)
                                                                                             agagCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAat       >  1:3356239/1‑100 (MQ=255)
                                                                                              gagCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAata      >  1:4337474/1‑100 (MQ=255)
                                                                                               agCACACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAatat     >  1:2909644/1‑100 (MQ=255)
                                                                                                  gcACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAATATATc  >  1:555815/1‑100 (MQ=255)
                                                                                                  gcACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAATATATc  >  1:3255417/1‑100 (MQ=255)
                                                                                                  gcACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAATATATc  >  1:3802050/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACACTAACAGCCAAAGGAATTGCGAAAGCGCCTAATGGCAAAGATTATGCGAGTACCATGAATTCGCCTTTAAAATATAATCCTAAAGAAGCAAGAGCGCACTGGGACAAAGCTAAAAAAGAGTTAGGTAAAAATGAAGTGACATTTTCAATGAACACAGAAGATACACCAGATGCAAAAATATCTGCTGAATATATC  >  NC_007793/979193‑979390

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