Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 90,315 (TA)8→6 coding (870‑873/1752 nt) SAUSA300_RS00415 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00779390,3020A.96.5% 100.5 / ‑2.1 29coding (857/1752 nt)SAUSA300_RS00415hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (28/0);  total (29/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_00779390,3030T.96.5% 100.6 / ‑2.4 29coding (858/1752 nt)SAUSA300_RS00415hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (28/0);  total (29/0)
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.40e-01
*NC_00779390,3040A.96.5% 101.0 / ‑2.7 29coding (859/1752 nt)SAUSA300_RS00415hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (28/0);  total (29/0)
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.40e-01
*NC_00779390,3050T.96.5% 100.4 / ‑2.1 29coding (860/1752 nt)SAUSA300_RS00415hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (28/0);  total (29/0)
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.90e-01

TATCGAATTCTTCAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAATATATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAATAG  >  NC_007793/90203‑90391
                                                                                                   ||||                                                                                      
tgccGAATTCCGCAACTACAGAAAATATATCTACAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCaa                                                                                           >  1:5038410/4‑100 (MQ=255)
           tcAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatat                                                                                >  1:4197187/1‑100 (MQ=255)
            cAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatata                                                                               >  1:1142878/1‑100 (MQ=255)
            cAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatata                                                                               >  1:4126665/1‑100 (MQ=255)
            cAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatata                                                                               >  1:3656550/1‑100 (MQ=255)
             aaCTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatataa                                                                              >  1:1678146/1‑99 (MQ=255)
              aCTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAatatatatatataaa                                                                             >  1:2595904/1‑98 (MQ=255)
                                ataGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATAc                                                       >  1:1388836/1‑100 (MQ=255)
                                ataGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATAc                                                       >  1:4903898/1‑100 (MQ=255)
                                ataGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATAc                                                       >  1:5448461/1‑100 (MQ=255)
                                ataGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATAc                                                       >  1:2913005/1‑100 (MQ=255)
                                   actatatCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACAtt                                                    >  1:5562223/3‑100 (MQ=255)
                                                   tAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCg                                    >  1:1029952/1‑100 (MQ=255)
                                                    aGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑‑‑ATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACt                                 >  1:1142419/1‑100 (MQ=255)
                                                         tAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAATATATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGAc                                  >  1:1612975/1‑100 (MQ=255)
                                                          taGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATc                             >  1:3047421/2‑100 (MQ=255)
                                                          aaGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATc                             >  1:5301943/1‑100 (MQ=255)
                                                          aaGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACAGATc                             >  1:4476333/1‑100 (MQ=255)
                                                          aaGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATACAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATc                             >  1:2323682/1‑100 (MQ=255)
                                                          aaGTAATATGATACTCTTGGTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATGAGATGCAACTGATc                             >  1:3085067/1‑100 (MQ=255)
                                                                 atGATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGt                      >  1:1246895/1‑100 (MQ=255)
                                                                   gATACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTg                    >  1:5139912/1‑100 (MQ=255)
                                                                   caTACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTg                    >  1:3907694/2‑100 (MQ=255)
                                                                    gtACTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGa                   >  1:4629271/2‑100 (MQ=255)
                                                                      aCTCTTGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGAtt                 >  1:4572586/1‑100 (MQ=255)
                                                                           tGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATc            >  1:3859078/1‑100 (MQ=255)
                                                                           tGTTAGATGATATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATc            >  1:4079309/1‑100 (MQ=255)
                                                                                 atgatATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGaa      >  1:2179855/1‑100 (MQ=255)
                                                                                   tatatGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGaaaa    >  1:1162884/2‑99 (MQ=255)
                                                                                   gatATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAAt    >  1:434779/1‑100 (MQ=255)
                                                                                   gatATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAAt    >  1:3945753/1‑100 (MQ=255)
                                                                                   gatATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAAt    >  1:2746513/1‑100 (MQ=255)
                                                                                   gatATGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAAt    >  1:2091175/1‑100 (MQ=255)
                                                                                     tatGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAATAg  >  1:2472647/1‑100 (MQ=255)
                                                                                     tatGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAATAg  >  1:839426/1‑100 (MQ=255)
                                                                                     tatGCTCGACGTCA‑‑‑‑ATATATATATATAAATTCTTTTAAAAATAATACACTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAATAg  >  1:745833/1‑100 (MQ=255)
                                                                                                   ||||                                                                                      
TATCGAATTCTTCAACTACAGAAAATATATCTATAGATATATCAAGATTTTTAGAACTAAGTAATATGATACTCTTGTTAGATGATATGCTCGACGTCAATATATATATATATATAAATTCTTTTAAAAATAATACATTAGAGACATTTGATGCGACTGATCCTGAAGTTGATTACATCACTGAAATAG  >  NC_007793/90203‑90391

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