Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,100,348	G→T	C127* (TGC→TGA)	ansB ←	L‑asparaginase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,100,348	0	G	T	97.5%	126.3 / ‑2.6	40	C127* (TGC→TGA)	ansB	L‑asparaginase 2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (19/20);  total (20/20)
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.21e-01

ATCACTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCGCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCA  >  NC_000913/3100214‑3100481                                                                                                                                       |                                                                                                                                      atcaCTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCGCAttt                                                                                                                                  >  1:504844/1‑140 (MQ=255)     cTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCAttt                                                                                                                                  <  2:820028/136‑1 (MQ=255)     cTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCAttt                                                                                                                                  >  1:820028/1‑136 (MQ=255)      tACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGTGCCGACCATCACCACCGGTTTGTCTCATTTCACCg                                                                                                                             <  1:515535/140‑1 (MQ=255)       aCCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGt                                                                                                                            <  2:352566/140‑1 (MQ=255)         cAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCa                                                                                                                          <  1:640834/140‑1 (MQ=255)           gCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAgg                                                                                                                        >  2:112211/1‑140 (MQ=255)           gCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAgg                                                                                                                        >  2:941793/1‑140 (MQ=255)           gCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAgg                                                                                                                        >  2:359315/1‑140 (MQ=255)           gCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAgg                                                                                                                        >  1:827441/1‑140 (MQ=255)           gCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAgg                                                                                                                        <  1:493733/140‑1 (MQ=255)              cGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCg                                                                                                                     <  1:359315/140‑1 (MQ=255)              cGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCg                                                                                                                     <  2:307053/140‑1 (MQ=255)                        ggCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAg                                                                                                           <  2:630744/140‑1 (MQ=255)                          cGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCa                                                                                                         <  1:102764/140‑1 (MQ=255)                                     tCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCAt                                                                                              <  2:843552/140‑1 (MQ=255)                                     tCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCAt                                                                                              <  2:827441/140‑1 (MQ=255)                                     tCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCAt                                                                                              <  2:577324/140‑1 (MQ=255)                                       aGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCg                                                                                            <  1:569821/140‑1 (MQ=255)                                         cTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCgtg                                                                                          >  1:57871/1‑140 (MQ=255)                                                 aCTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTAcc                                                                                  <  1:112211/140‑1 (MQ=255)                                                 aCTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTAcc                                                                                  <  1:941793/140‑1 (MQ=255)                                                 aCTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTAcc                                                                                  <  1:621719/140‑1 (MQ=255)                                                           tataCAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCgtgt                                                                                         <  1:757768/123‑1 (MQ=255)                                                           tataCAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCgtgt                                                                                         >  2:757768/1‑123 (MQ=255)                                                             taCAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGAc                                                                      >  1:732447/1‑140 (MQ=255)                                                                  gTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGc                                                                 >  2:437840/1‑140 (MQ=255)                                                                     gAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGCTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGt                                                              >  1:571811/1‑140 (MQ=255)                                                                       aTGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGaa                                                                   <  1:614793/133‑1 (MQ=255)                                                                       aTGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGaa                                                                   >  2:614793/1‑133 (MQ=255)                                                                       aTGCACCGTGTGCGCTCATAGACGTGGACGGACGCATTGCGGCGACCATCACCACCGGTTTGTCTCATTGCACCGTCAGGTCGCGGAAGTAAGCAGTTTCTTCCAGCGTGTCGGAACCGTGGGTAATGACGAAGCCGTCg                                                            >  2:665280/1‑140 (MQ=255)                                                                        tGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGCGGGTAATGACGAAGCCGTCgg                                                           >  1:923614/1‑140 (MQ=255)                                                                          gACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTc                                                         >  1:674437/1‑140 (MQ=255)                                                                                cTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCg                                                   >  2:247403/1‑140 (MQ=255)                                                                                cTGCCCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGGAATGACGAAGCCGTCGGTCTTATCg                                                   >  2:125051/1‑140 (MQ=255)                                                                                                ggacggacGCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAAttt                                   >  2:467110/1‑140 (MQ=255)                                                                                                       cgCATTGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGc                            <  1:467110/140‑1 (MQ=255)                                                                                                           ttGCGCCGACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGAAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAgt                        >  2:454306/1‑140 (MQ=255)                                                                                                                  gACCATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGCGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGcc                 <  2:923614/140‑1 (MQ=255)                                                                                                                    ccATCACCACCGGTTTGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAg               >  1:123102/1‑140 (MQ=255)                                                                                                                                 tttGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCa  >  2:388722/1‑140 (MQ=255)                                                                                                                                 tttGTCTCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCa  >  2:365411/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                      ATCACTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCGGTTTGTCGCATTTCACCGTCAGGTCGAGGAAGTAAGCAGTTTCTTCCATCGTGTCGGTACCGTGGGTAATGACGAAGCCGTCGGTCTTATCGCAGTCGGTGTTAATTTTTTTCGCCAGTGTCAGCCAGACATTATCGTTCA  >  NC_000913/3100214‑3100481

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