Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,579,248	C→T	T2T (ACG→ACA)	SAUSA300_RS07725 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,579,248	0	C	T	97.1%	86.6 / NA	35	T2T (ACG→ACA)	SAUSA300_RS07725	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (1/33);  total (2/33)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCAACAACTGTAAACCTTTGATTGCTTTTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCATGAT  >  NC_007793/1579153‑1579339                                                                                                |                                                                                            gCCTCAACAACTGTAAACCTTTGATTGCTATTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTct                                                                                         <  1:4831664/100‑2 (MQ=255) gCCTCAACAACTGTAAACCTTTGATTGCTATTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCa                                                                                         >  1:1542752/1‑100 (MQ=255)             gTAAACCTTTGATTGCTATTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGCGTCACTTCCCCAAAAc                                                                             >  1:1090787/1‑100 (MQ=255)                     ttGATTGCTATTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAACACCTCCTTgt                                                                     <  1:4338468/100‑2 (MQ=255)                                                tgtATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTaa                                          <  1:2917816/100‑2 (MQ=255)                                                             ccTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTTCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGt                             <  1:6432353/100‑1 (MQ=255)                                                              cTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTc                            <  1:471360/100‑1 (MQ=255)                                                              cTGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTc                            <  1:2176567/100‑1 (MQ=255)                                                              cTGTTGAATCTTTGAATGTTGTGACTAAGTACTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTATTTACACGTATCCTTTTCCTGCGTc                            <  1:4621717/100‑1 (MQ=37)                                                               tGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTcg                           <  1:5562565/100‑2 (MQ=255)                                                               tGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTcg                           <  1:5374596/100‑2 (MQ=255)                                                               tGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTcg                           <  1:389966/100‑2 (MQ=255)                                                               tGTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTca                           <  1:1127241/100‑2 (MQ=255)                                                               tGTTGAATCTTTGAATGTTGTGACTAAGTACTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGtat                           <  1:1697150/100‑3 (MQ=37)                                                                gTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTTGATCTAAAATGTCTTTACGCGTATCCTTTTCCTGCGTcgt                          <  1:5739696/100‑3 (MQ=37)                                                                gTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTcgg                          <  1:3211338/100‑3 (MQ=255)                                                                gTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTcgt                          <  1:2807061/100‑3 (MQ=255)                                                                gTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTcaa                          <  1:2661129/100‑3 (MQ=255)                                                                gTTGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCGTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTctg                          <  1:895227/100‑1 (MQ=255)                                                                gTTGAATCTTTGAAGGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATGTAAAATGTCTTTACACGTATCCTTTTCCTGCGTcag                          <  1:4811840/100‑3 (MQ=37)                                                                 ttGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCCAAGATGTCTTTACACGTATCCTTTTCCTGCGTctgc                         <  1:3743963/100‑1 (MQ=255)                                                                  tGAATCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTACCCTTTTCCTGCGTctgcc                        <  1:4411241/100‑2 (MQ=255)                                                                    aaTCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAGCCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGt                      <  1:6199413/100‑1 (MQ=255)                                                                     cttccgatctgttgtGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGtt                     <  1:4224362/91‑1 (MQ=25)                                                                     aTCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGtt                     <  1:4608897/100‑1 (MQ=255)                                                                     aTCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGtt                     <  1:5807359/100‑1 (MQ=255)                                                                     aTCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGtc                     <  1:1883232/100‑2 (MQ=255)                                                                      tCTTTGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGttg                    <  1:1213898/100‑2 (MQ=255)                                                                         ttGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCAt                 <  1:6252780/100‑1 (MQ=255)                                                                          tGAATGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCAta                <  1:6313222/100‑2 (MQ=255)                                                                            aaTGTTGTGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCta              <  1:1747848/100‑2 (MQ=255)                                                                             acgttgtGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTg             <  1:1373300/98‑1 (MQ=255)                                                                              tgttgtGACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTgg            <  1:256324/100‑2 (MQ=255)                                                                                 tgtgACTAAGTATTGTGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCtt         <  1:6488380/100‑1 (MQ=255)                                                                                        aaGTATTGTGTCATTTCCCCAAAACCTCCTCGACTCGATCTAAAATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCGTggt  <  1:6414682/100‑3 (MQ=255)                                                                                                |                                                                                            GCCTCAACAACTGTAAACCTTTGATTGCTTTTAGCTTTAGTTATGTGTGTATGCTTGCGTCCTGTTGAATCTTTGAATGTTGTGACTAAGTATTGCGTCACTTCCCCAAAACCTCCTTGACTCGATCTAAGATGTCTTTACACGTATCCTTTTCCTGCGTCTGCTGTTCCATCTTGTCTTTCATGAT  >  NC_007793/1579153‑1579339

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