Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	529,067	C→T	P112L (CCC→CTC)	SAUSA300_RS02520 →	ribonuclease M5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	529,067	0	C	T	98.1%	143.8 / NA	52	P112L (CCC→CTC)	SAUSA300_RS02520	ribonuclease M5
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/6);  minor base G (1/0);  total (46/6)
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.99e-01

ATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGAATTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGTAGAGAAATTTTAAG  >  NC_007793/528970‑529163                                                                                                  |                                                                                                 atTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCActct                                                                                                >  1:5757832/1‑100 (MQ=255)  gtGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCtt                                                                                               >  1:1769186/2‑100 (MQ=255)      agtaactaaaGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTtgat                                                                                           >  1:3627158/8‑100 (MQ=255)             aaaGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCtt                                                                                    >  1:1203291/1‑100 (MQ=255)                gCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATg                                                                                 <  1:1834830/100‑1 (MQ=255)                   aaaaataaaaaAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAAt                                                                              >  1:6219496/1‑100 (MQ=255)                   aaaaataaaaaAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAAt                                                                              >  1:2425113/1‑100 (MQ=255)                         aaaaaaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTg                                                                        >  1:1320521/1‑100 (MQ=255)                         aaaaaaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTg                                                                        >  1:1448148/1‑100 (MQ=255)                            tcaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATa                                                                     >  1:3091640/3‑100 (MQ=255)                            taaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATa                                                                     >  1:3694354/2‑100 (MQ=255)                            aaaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATa                                                                     >  1:2277724/1‑100 (MQ=255)                             aaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATaa                                                                    >  1:876208/1‑100 (MQ=255)                             aaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATaa                                                                    >  1:6010173/1‑100 (MQ=255)                             aaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATaa                                                                    >  1:217276/1‑100 (MQ=255)                             aaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACGCTTTGATGAAGCTTATGAATCAATTGATaa                                                                    >  1:1984597/1‑100 (MQ=255)                              aGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATaaa                                                                   >  1:2940749/1‑100 (MQ=255)                              aGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATaaa                                                                   >  1:5549222/1‑100 (MQ=255)                               gggAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAAt                                                                  >  1:4256178/1‑100 (MQ=255)                                    gaTTGGTGTTGAACATGCCGAATTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCtgtg                                                             >  1:5206854/2‑100 (MQ=255)                                     tttGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGc                                                            >  1:5264468/2‑100 (MQ=255)                                      ttGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCt                                                           <  1:470808/100‑1 (MQ=255)                                       tGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTa                                                          >  1:6339563/1‑100 (MQ=255)                                       tGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTa                                                          >  1:3052207/1‑100 (MQ=255)                                        ggTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTaa                                                         >  1:2165631/1‑100 (MQ=255)                                         gtgtTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAAt                                                        <  1:641908/100‑1 (MQ=255)                                             tGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATagag                                                    >  1:2169589/1‑100 (MQ=255)                                             gctaCATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATagag                                                    >  1:3615773/4‑100 (MQ=255)                                              gAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGt                                                   >  1:5773212/1‑100 (MQ=255)                                               gaCATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGtt                                                  >  1:4253974/2‑100 (MQ=255)                                                 cttGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTgg                                                >  1:1139886/3‑100 (MQ=255)                                                 aaTGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTgg                                                >  1:5831675/2‑100 (MQ=255)                                                  aTGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTggg                                               >  1:4637105/1‑100 (MQ=255)                                                    gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTgggcc                                             >  1:3949569/1‑98 (MQ=255)                                                    gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt                                             >  1:3667482/1‑100 (MQ=255)                                                    gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt                                             >  1:3646595/1‑100 (MQ=255)                                                    gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt                                             >  1:6383569/1‑100 (MQ=255)                                                    gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTGATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt                                             >  1:3255965/1‑100 (MQ=255)                                                            aattGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAAttgtt                                     >  1:6618362/1‑100 (MQ=255)                                                            aattGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAAttgtg                                     >  1:2849210/1‑99 (MQ=255)                                                            aattGATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAAttgtc                                     >  1:354882/1‑99 (MQ=255)                                                                gATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTGATAGAGTTGGGGTTAATTGTTGGGa                                 >  1:1333816/1‑100 (MQ=255)                                                                gATATTAAAGAAGCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAAttgtcagga                                 <  1:2344804/100‑6 (MQ=255)                                                                          aagCGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTACTTGTCGGGAAAGAtcggaa                       >  1:310630/1‑95 (MQ=255)                                                                            ccGTTAATGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGGGTTGGGGTTAATTGTTGGGAAAGAtcggaaga                     >  1:6268899/2‑93 (MQ=255)                                                                                 aaTGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCGATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGt                >  1:1636841/1‑100 (MQ=255)                                                                                 aaTGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGt                <  1:6183531/100‑1 (MQ=255)                                                                                 aaTGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGt                >  1:1453030/1‑100 (MQ=255)                                                                                 aaTGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTGAATTGTTGGGAAAGATGCAAGGCGCCGt                <  1:554754/100‑1 (MQ=255)                                                                                   tGCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGTag              >  1:3033341/1‑100 (MQ=255)                                                                                    gCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGTaga             >  1:1335565/1‑100 (MQ=255)                                                                                    gCATGTTAGTTCACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATACAAGGCGCCGTaga             >  1:4896589/1‑100 (MQ=255)                                                                                               cACTCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTCTTAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGTAGAGAAATTTTAAg  >  1:3257397/1‑100 (MQ=255)                                                                                                  |                                                                                                 ATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGAATTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTAATTGTTGGGAAAGATGCAAGGCGCCGTAGAGAAATTTTAAG  >  NC_007793/528970‑529163

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