Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,736,609	G→A	P190L (CCA→CTA)	bax ←	putative glucosaminidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,736,609	0	G	A	96.7%	193.9 / ‑1.7	60	P190L (CCA→CTA)	bax	putative glucosaminidase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (27/31);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCT  >  NC_000913/3736512‑3736701                                                                                                  |                                                                                             gAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAgc                                                                                            <  2:1348585/100‑1 (MQ=255)    cGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCa                                                                                         >  1:203801/1‑100 (MQ=255)     gAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCAt                                                                                        >  1:211372/1‑100 (MQ=255)       agagTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTg                                                                                      >  1:802541/1‑100 (MQ=255)          gTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTa                                                                                   <  1:455876/100‑1 (MQ=255)          gTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTa                                                                                   >  1:73679/1‑100 (MQ=255)          gTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTa                                                                                   <  2:175818/100‑1 (MQ=255)          gTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTa                                                                                   <  2:1099600/100‑1 (MQ=255)           tAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTAc                                                                                  <  1:21889/100‑1 (MQ=255)               ccGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACg                                                                              <  1:52258/100‑1 (MQ=255)                  ggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTcc                                                                           <  1:1280608/100‑1 (MQ=255)                  ggTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTcc                                                                           >  1:1292345/1‑100 (MQ=255)                     gCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCttt                                                                        >  2:612687/1‑100 (MQ=255)                      cGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTc                                                                       >  2:1335789/1‑100 (MQ=255)                        tATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCAt                                                                     <  1:1365012/100‑1 (MQ=255)                        tATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCAt                                                                     <  2:370383/100‑1 (MQ=255)                         aTTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATg                                                                    <  1:302251/100‑1 (MQ=255)                         aTTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATg                                                                    <  1:767102/100‑1 (MQ=255)                          ttCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGc                                                                   >  1:1433321/1‑100 (MQ=255)                              gatctgTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCAttt                                                               <  2:1050497/95‑1 (MQ=255)                                 tAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCAt                                                            <  1:1208768/100‑1 (MQ=255)                                  aGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATg                                                           <  2:1375876/100‑1 (MQ=255)                                   gTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCAtttagatc                                                          >  1:1050504/1‑95 (MQ=255)                                   gTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGc                                                          <  2:580793/100‑1 (MQ=255)                                     gACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCg                                                        <  2:356651/100‑1 (MQ=255)                                          aGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGACCAg                                                   <  2:857152/100‑1 (MQ=255)                                              gcTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttg                                               >  1:398652/1‑100 (MQ=255)                                                 cACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttg                                            <  2:1292336/100‑1 (MQ=255)                                                  aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgt                                           >  2:1047439/1‑100 (MQ=255)                                                  aCCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgt                                           <  2:73679/100‑1 (MQ=255)                                                    cGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGgttgttgttg                                         <  2:303187/100‑1 (MQ=255)                                                       ttCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGttgttgttgttg                                      <  2:606272/100‑1 (MQ=255)                                                        tCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGCTTGTTGTTGTTgc                                     >  1:914295/1‑100 (MQ=255)                                                         cTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcg                                    >  1:1482740/1‑100 (MQ=255)                                                          tttGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgc                                   >  1:974044/1‑100 (MQ=255)                                                          tttGACAGAACTAAACTGTGAGTACCCTGTCACTTTACCTAGCGCATTGGTACATCGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgc                                   <  2:398952/100‑1 (MQ=255)                                                           ttGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTgcgcg                                  >  1:312347/1‑100 (MQ=255)                                                             gACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGcc                                >  2:400617/1‑100 (MQ=255)                                                              aCAGAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTTTTGTTGCGCGCCa                               <  2:974036/100‑1 (MQ=255)                                                                 gAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCt                            >  2:294866/1‑100 (MQ=255)                                                                 gAACTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCt                            >  1:158551/1‑100 (MQ=255)                                                                    cTAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCg                         >  2:178299/1‑100 (MQ=255)                                                                     tAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCTCGCCAGCTTCGa                        <  1:1030766/100‑1 (MQ=255)                                                                     tAAACTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGa                        <  1:1257019/100‑1 (MQ=255)                                                                      aaaCTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                       <  2:977012/100‑1 (MQ=255)                                                                      aaaCTGTGAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGAc                       >  2:786690/1‑100 (MQ=255)                                                                          tgtgAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTc                   <  1:376019/100‑1 (MQ=255)                                                                           gtgAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcc                  <  2:1309751/100‑1 (MQ=255)                                                                             gAGTACCCTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTcccc                <  1:531950/100‑1 (MQ=255)                                                                                  cccTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCg           >  2:1396540/1‑100 (MQ=255)                                                                                   ccTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGc          >  1:223208/1‑100 (MQ=255)                                                                                    cTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCt         >  2:50072/1‑100 (MQ=255)                                                                                    cTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCt         >  2:210855/1‑100 (MQ=255)                                                                                    cTTTCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCAAGCTTCGACGTTCCCCAACCGCt         >  2:935848/1‑100 (MQ=255)                                                                                     tttCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCtt        <  1:1381477/100‑1 (MQ=255)                                                                                       tCACTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTc      <  1:1232774/100‑1 (MQ=255)                                                                                         aCTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctg    >  2:1333344/1‑100 (MQ=255)                                                                                          cTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgc   <  1:1333353/100‑1 (MQ=255)                                                                                          cTTTACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctcc   >  2:684581/1‑98 (MQ=255)                                                                                           tttACCTAGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTctgct  <  2:1085710/100‑1 (MQ=255)                                                                                                  |                                                                                             GAACGAAGAGTAAGCCGGGTGCGTATTCAGGTTAGTGACATAGGCGCTCACCGATTCTTTGACAGAACTAAACTGTGAGTACCCTTTCACTTTACCTGGCGCATTGGTACAACGTCCTTTCATGCATTTCATGCCGAACAGGTTGTTGTTGTTGCGCGCCAGCTTCGACGTTCCCCAACCGCTTTCTGCT  >  NC_000913/3736512‑3736701

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: ‑