Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,965,820	G→A	A20A (GCG→GCA)	trxA →	thioredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,965,820	0	G	A	95.6%	145.4 / ‑1.6	45	A20A (GCG→GCA)	trxA	thioredoxin 1
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (22/21);  total (23/22)
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.04e-01

ACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGG  >  NC_000913/3965726‑3965915                                                                                               |                                                                                                aCCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACgg                                                                                            <  1:206094/100‑1 (MQ=255)   cAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACgggg                                                                                          >  2:980697/1‑100 (MQ=255)        cGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATc                                                                                     >  2:394779/1‑100 (MQ=255)         gCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATcc                                                                                    >  2:410342/1‑100 (MQ=255)         gCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATcc                                                                                    <  1:157219/100‑1 (MQ=255)          ccAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCt                                                                                   >  1:365452/1‑100 (MQ=255)               cTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCtcgtcg                                                                              >  1:1151358/1‑100 (MQ=255)                    tCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGAttt                                                                          >  1:1063575/1‑99 (MQ=255)                     ccTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCt                                                                        >  1:688671/1‑100 (MQ=255)                      cTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTg                                                                       <  1:656090/100‑1 (MQ=255)                        gtgGAGTTATATATGAGCGATAAAATTATTCCCCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTggg                                                                     <  1:655603/100‑1 (MQ=255)                        gtgGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTggg                                                                     <  1:176087/100‑1 (MQ=255)                          ggAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCa                                                                   <  1:785577/100‑1 (MQ=255)                          ggAGTTATATATGAGCGATAAAATGATGCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGc                                                                    <  1:305760/99‑1 (MQ=255)                                atatatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAgtgg                                                             <  2:1063565/100‑1 (MQ=255)                                   tatGAGCGATAAAATTATTCACGTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                          <  2:881382/100‑1 (MQ=255)                                   tatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                          <  2:614983/100‑1 (MQ=255)                                   tatGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGc                                                          <  2:11117/100‑1 (MQ=255)                                    atGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCg                                                         >  2:815827/1‑100 (MQ=255)                                     tGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCgg                                                        <  1:821439/100‑1 (MQ=255)                                      gAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGt                                                       >  2:666697/1‑100 (MQ=255)                                       aGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTc                                                      >  1:47886/1‑100 (MQ=255)                                            tAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGc                                                 <  2:956069/100‑1 (MQ=255)                                            tAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGc                                                 <  1:1181599/100‑1 (MQ=255)                                            tAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGc                                                 >  2:65758/1‑100 (MQ=255)                                               aaTTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCaaa                                              <  2:1190158/100‑1 (MQ=255)                                                     tCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATc                                        >  2:404960/1‑100 (MQ=255)                                                     tCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATc                                        >  1:736761/1‑100 (MQ=255)                                                        ccTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcc                                     >  2:1147248/1‑100 (MQ=255)                                                         ctgactgaCGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGccc                                    <  2:480471/100‑1 (MQ=255)                                                         ctgactgaCGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGccc                                    >  1:352544/1‑100 (MQ=255)                                                          tgactgacGACAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGcccc                                   <  1:815833/100‑1 (MQ=255)                                                                 cgacAGTTTTGACACGGATGTACTCAAAGCATACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGGTTCTg                            >  2:1250199/1‑100 (MQ=255)                                                                  gacAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTgg                           <  2:736756/100‑1 (MQ=255)                                                                   acAGTTTTGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGa                          <  1:719904/100‑1 (MQ=255)                                                                          tGACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTTGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATc                   >  2:253711/1‑100 (MQ=255)                                                                           gACACGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCg                  >  2:535851/1‑100 (MQ=255)                                                                            acacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGc                 >  1:470960/1‑100 (MQ=255)                                                                             cacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCt                >  1:761652/1‑100 (MQ=255)                                                                             cacGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCt                <  2:716304/100‑1 (MQ=255)                                                                              acGGATGTACTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTg               <  1:1236240/100‑1 (MQ=255)                                                                                      aCTCAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAatat       >  1:5881/1‑100 (MQ=255)                                                                                         cAAAGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg    >  1:17997/1‑100 (MQ=255)                                                                                         cAAAGCAGACGGGGCGATCCTCGTCGATTTCTAGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAg    >  1:1238257/1‑100 (MQ=255)                                                                                           aaGCAGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAggg  <  2:222501/100‑1 (MQ=255)                                                                                               |                                                                                                ACCAACACGCCAGGCTTATTCCTGTGGAGTTATATATGAGCGATAAAATTATTCACCTGACTGACGACAGTTTTGACACGGATGTACTCAAAGCGGACGGGGCGATCCTCGTCGATTTCTGGGCAGAGTGGTGCGGTCCGTGCAAAATGATCGCCCCGATTCTGGATGAAATCGCTGACGAATATCAGGG  >  NC_000913/3965726‑3965915

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