Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,913,228	C→T	D176N (GAT→AAT)	SAUSA300_RS09460 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,913,228	0	C	T	86.7%	17.7 / ‑3.2	15	D176N (GAT→AAT)	SAUSA300_RS09460	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (0/13);  total (1/14)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

GATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTT‑CCATTGGTGGT‑AAATCTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAGTTCTTGTTAATACTAAATGACT  >  NC_007793/1913153‑1913324                                                                              |                                                                                                  gaaccATAACAATATTTTGGGGTCTAGATTTATCAAAATCAAATGTATCAAATACTTTTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATAcc                                                                          <  1:2757993/97‑1 (MQ=255)   tCCATCACAATAGTTTGCGGTCTAGATTTATCAAAATCCAATGTAAAAAATACTTGTT‑CCATTGGTGGA‑AAATCTTAAAATTG‑ACCGCCACTAATACCAt                                                                        <  1:202928/100‑1 (MQ=255)    ccATCACAATATTTTGCGGTCTAGATTTATCAAAATCAAATGTATAAAATACTTGTT‑CCATTGGTGGT‑AAATTTTAAAATTG‑ACCGCCACTAATACCAtt                                                                       <  1:5851252/100‑1 (MQ=255)    ccATCACAATAGTTTGCGGTCTAGATTTATCAAAATCCAATGTATAAAATTCTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCAtt                                                                       <  1:2564393/100‑1 (MQ=255)     cagcacaATTGTTTGCGGTCGAGATTTATATAACACAAATATATC‑AATTGTTGTTGCCGTTGGTGGT‑AATCTTTTAAATTG‑TCCGCCACTAATACCAtct                                                                      <  1:2093211/97‑3 (MQ=255)      ttCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTT‑CCATTGGTGGT‑AAATCTTTAAATTG‑ACCGCCACTAATACCATtat                                                                     >  1:990220/2‑100 (MQ=255)             tAGTTTGCGTTTTAGTTTTATCTAAATCCAATGTAAAAAATTTTTTTT‑CCTTTGGGGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATg                                                              <  1:3689528/100‑1 (MQ=255)             tAGTTTGCGGTCTAGATTTATCTAAATCAAATGTATAAAATCCTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATg                                                              <  1:1744067/100‑1 (MQ=255)             tAGTTTGCGGTCTAGATTTATCAAAATCCAATGTATCAAATACTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATg                                                              <  1:3152049/100‑1 (MQ=255)               gTTGGCGGTTTAGATTTATTTAAATCCAATGTATAAAATATTTTTC‑CCATTGGTGTT‑AACTTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATgcg                                                            <  1:3549121/100‑3 (MQ=255)                 ttGCGGTCTATATTTATTAAAAGCAAATGTCTAAAATATTTTTT‑CCATTGGGGGA‑AAATCTTAAAATTG‑CCCGCCACTAATACCATTATAAACATGAcca                                                          <  1:1950083/100‑2 (MQ=255)                 ttGCGGTCAAGTTTTATCAAAATCCAATTTATCAAATACTTGTT‑CCATTGGGGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCt                                                          <  1:939250/100‑1 (MQ=255)                  ggCGGTCTAGATTTATCAAAATCGAATGTATAAAATATTTGTT‑CCATTGGGGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCtt                                                         <  1:5185849/99‑1 (MQ=255)                    cGGTTTATATTTTTATAAATCCAATGTTTAAAATTCTTGTT‑CCATTGGTGGT‑AAGTTTTTAATTTG‑ACCGCCACTAATACCATTATAAACATGACCtttt                                                       <  1:4078155/100‑1 (MQ=255)                    cGGTCTCGTTTTATCTAATTAAAAGGTATAAAATAGTTTTT‑CTTTTGGAGGA‑AAATTTTAAAATGA‑ACCGCCACTAATACCATTATAAACATGACCtttt                                                       <  1:2142340/100‑1 (MQ=255)                     ggTCTAGATTTATCTAAATCAAATGTATAAAATACTTGTT‑CCTTTGGTGGA‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTa                                                      <  1:673226/100‑1 (MQ=255)                       tCTAGATTTATATAAATCCAATTTATAAAATTCTTGTT‑CCATTGGTGTT‑AATTTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAAc                                                    <  1:5585894/100‑1 (MQ=255)                            tttttTCTAACCCCAATGTATCAAATATTTCTA‑CATTGGGTGGT‑AACTTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAAtta                                               <  1:1463646/99‑2 (MQ=255)                            ttttATTAAAATCAATTGTATAAAATTCTTGTT‑CATTTGGTGGA‑AAATTTTAAAATGG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:2505333/99‑1 (MQ=255)                            ttttATCTAAACCCAATGTATAAAATACTTTTT‑CCTTTGGTGGT‑AAATTTTAAAATGG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:5308397/99‑1 (MQ=255)                            ttttATCAAAATCAAATGTATAAAATTCTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:1047994/99‑1 (MQ=255)                            ttttATCAAAATCAAATGTATAAAAT‑TTTTTT‑TCATTGGGGGTAAAATCTTTAAATGG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:736924/99‑1 (MQ=255)                            atttactaaaaCCCAATGTATCAAAT‑CTTGTTCCCATTGGTGGT‑AATTCTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:3822053/92‑1 (MQ=255)                            aTTTATCTAAATCCAATTTATCAAATACTTTTT‑CAATTGGTGGT‑AAATTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:2477130/100‑1 (MQ=255)                            aTTTATCTAAATCCAATGTATAAAATACTTTTT‑CCTTTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:5471102/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATCAAATCTTTTTT‑CCTTTGGGGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:5383219/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATCAAATCTTTGTT‑CCATGGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:541393/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATCAAATCCTTGTT‑CCATTGGTGGA‑AAATCTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:3400740/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATAAAATTTTTGTT‑CCTTTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTCTAACAATTg                                               <  1:2649768/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATAAAATATTTTTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:2538878/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATAAAATACTTTTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:3094882/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTATAAAATACTTGTC‑CCATGGGTGCT‑AGATTTTTAAATTG‑ACCGCCACTTACACCATTATAAACATGACCTTTTAACAATTg                                               <  1:4660544/100‑1 (MQ=255)                            aTTTATCTAAATCAAATGTACAAAATACTTGTT‑CTATGGGGGGT‑AAATTTTAAAATGA‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:1589308/100‑1 (MQ=255)                            aTTTATCAAAATCCAATGTATAAAATCTTTTTT‑CCATTGGTGGA‑AAATTTTAAAATGG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:2600440/100‑1 (MQ=255)                            aTTTATCAAAATCAATTTTATAAAATTTTTTTT‑CCATGGGTGGA‑AAATTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTg                                               <  1:2223817/100‑1 (MQ=255)                            aTTTATCAAAATCAAATGTATAAAATACTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAAtcg                                               <  1:1640401/100‑3 (MQ=255)                             tttATCTAAATCAAATGTATAAAATCCTTGTT‑CCATTGGTGGT‑AAATTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTga                                              <  1:3816658/100‑1 (MQ=255)                              ttATCTAAATCCAATGTATCAAATACTTGTT‑CCTTTGGTGGT‑AAATTTTTAATTTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTgag                                             <  1:1598664/100‑1 (MQ=255)                                    aaaTCCAATGTATAAAATACTTGTT‑CCATTGGGGGT‑AAATTTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCag                                       <  1:4802303/100‑2 (MQ=255)                                                 aaaaTACTTGTT‑CTTTTGGTGTT‑AAAT‑TTTAAATTGAACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGAcc                          <  1:2817398/99‑1 (MQ=255)                                                 aaaaTACTTGTT‑CCATTGGTGGA‑AATTTTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGAcc                          <  1:2082967/99‑1 (MQ=255)                                                   aatctttGTT‑CCATTGGTGAA‑AATTCTTTAAATGG‑CCCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAg                        <  1:1128402/95‑1 (MQ=255)                                                   aaTTCTTGTT‑CCTTGGGGGGA‑AAATTTTTAATTTG‑ACCGCCACTAATGCCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAg                        <  1:2806941/100‑1 (MQ=255)                                                      aCTTATA‑CCATTGGTGGA‑AAATTTTTAAGTTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAGtta                     <  1:1402440/100‑2 (MQ=255)                                                             tcaaTTGGGGGT‑AATTCTTAAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAGTTCTTGTTa               <  1:1021786/97‑1 (MQ=255)                                                                           aaTCTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAGTTCTTGTTGATACTAAATGACt  <  1:4619602/100‑1 (MQ=255)                                                                              |                                                                                                  GATCCATCACAATAGTTTGCGGTCTAGATTTATCTAAATCCAATGTATCAAATACTTGTT‑CCATTGGTGGT‑AAATCTTTAAATTG‑ACCGCCACTAATACCATTATAAACATGACCTTTTAACAATTGAGAATCAATAATATAAAGACCAGTTCTTGTTAATACTAAATGACT  >  NC_007793/1913153‑1913324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑