Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,669,920	A→G	E110G (GAA→GGA)	treF →	cytoplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,669,920	0	A	G	95.9%	153.5 / ‑1.6	49	E110G (GAA→GGA)	treF	cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (20/27);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

TTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGC  >  NC_000913/3669824‑3670015                                                                                                 |                                                                                                ttCCCGACTGCGCACCTAAAATGGTCGCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGgaa                                                                                               <  1:139755/99‑1 (MQ=255)   cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                            >  2:1552154/1‑100 (MQ=255)   cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                            >  2:1279208/1‑100 (MQ=255)   cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                            >  1:406942/1‑100 (MQ=255)      gACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACt                                                                                         >  1:1136832/1‑100 (MQ=255)            gcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGc                                                                                   <  1:1362895/100‑1 (MQ=255)              aCCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTg                                                                                 >  1:998937/1‑100 (MQ=255)               ccTAAAATGAGCCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGCTTTTGACTTGCGCAAGTTTGTTGGAAGCCACTTCTGTCTGc                                                                                <  2:1130626/100‑1 (MQ=255)                    aaTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAg                                                                           <  1:1009832/100‑1 (MQ=255)                     aTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAgg                                                                          <  1:1334054/100‑1 (MQ=255)                       ggACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCCGGCTGCCGGAGGTc                                                                        <  2:861195/100‑1 (MQ=255)                        gACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCt                                                                       >  1:181614/1‑100 (MQ=255)                        gACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCt                                                                       >  1:811080/1‑100 (MQ=255)                           ccGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACt                                                                    <  2:908776/100‑1 (MQ=255)                                ggATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGGCTACTCCAg                                                                >  2:190238/1‑99 (MQ=255)                                 gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                              <  1:8517/100‑1 (MQ=255)                                     tCTTAATCCGCTCCCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                          <  1:1348618/100‑1 (MQ=255)                                     tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                          <  1:631480/100‑1 (MQ=255)                                     tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                          >  2:827698/1‑100 (MQ=255)                                     tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                          <  1:876876/100‑1 (MQ=255)                                      cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                         >  2:139966/1‑100 (MQ=255)                                      cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                         <  2:578959/100‑1 (MQ=255)                                      cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                         >  2:864952/1‑100 (MQ=255)                                       ttAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatg                                                        <  1:950273/100‑1 (MQ=255)                                          ctCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCACGTTTGTTGGACACCACTTCTTGCTGCCGGAGGTCTACTCCAGCGaccatg                                                        >  2:1090557/2‑92 (MQ=255)                                          aTCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgtat                                                     >  2:1294892/1‑100 (MQ=255)                                            ccGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATc                                                    >  1:1173011/1‑99 (MQ=255)                                             cGCTACCGTAAATTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCgg                                                  <  2:544563/100‑1 (MQ=255)                                             cGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCgg                                                  >  1:216328/1‑100 (MQ=255)                                                tACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGAcc                                               <  1:281595/100‑1 (MQ=255)                                                tACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGAcc                                               <  2:1033790/100‑1 (MQ=255)                                                  ccGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCg                                             >  2:1558119/1‑100 (MQ=255)                                                  ccGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCg                                             <  2:175744/100‑1 (MQ=255)                                                           gcgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTcc                                    >  2:1421241/1‑100 (MQ=255)                                                            cgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                   <  1:897428/100‑1 (MQ=255)                                                             gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAATATTCCCt                                  <  1:827703/100‑1 (MQ=255)                                                             gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                  <  2:652567/100‑1 (MQ=255)                                                               cgtcatcgtcatTTTGACTTGCGCAAGTTTGTTGGAACCCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGa                                <  2:938416/100‑1 (MQ=255)                                                                gTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGaa                               <  2:1485601/100‑1 (MQ=255)                                                                gTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGaa                               <  1:1250183/100‑1 (MQ=255)                                                                   aTCGTGATTTTGACTTGCGCAAGTTTGTTGGATACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAag                             <  1:1552933/99‑1 (MQ=255)                                                                      gTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCa                         <  2:354998/100‑1 (MQ=255)                                                                      gTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCa                         <  2:835144/100‑1 (MQ=255)                                                                          ttttGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATc                     >  1:1270065/1‑100 (MQ=255)                                                                                 ttGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGc              <  2:414096/100‑1 (MQ=255)                                                                                         gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCg      >  1:1090719/1‑100 (MQ=255)                                                                                          tttGTTGGAAACCATTTCTGGCTGCCGGAGGTCTTCGCCAGCGAGTATGTATCGGACCCGCAAAATGCCCTGAAAGAGCATATCGACCAGCTGTGGCCgg     <  2:763634/100‑1 (MQ=255)                                                                                             gttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGc  >  2:1243990/1‑100 (MQ=255)                                                                                             gttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGc  >  2:1155151/1‑100 (MQ=255)                                                                                             gttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCAGTGGCCGGTGc  <  2:1213133/100‑1 (MQ=255)                                                                                                 |                                                                                                TTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGC  >  NC_000913/3669824‑3670015

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