Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,239	T→C	H158H (CAT→CAC)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,239	0	T	C	98.0%	191.7 / NA	51	H158H (CAT→CAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/28);  minor base A (1/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCT  >  NC_000913/3026105‑3026373                                                                                                                                       |                                                                                                                                       gAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCa                                                                                                                                    <  1:503395/139‑1 (MQ=255)   ggATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                                  <  1:822812/139‑1 (MQ=255)   ggATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                                  <  2:427688/139‑1 (MQ=255)       ttAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGc                                                                                                                              >  2:270829/1‑139 (MQ=255)             tACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGa                                                                                                                        >  1:1060867/1‑139 (MQ=255)              aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                       <  2:846953/139‑1 (MQ=255)              aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                       >  1:846953/1‑139 (MQ=255)              aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                       <  1:398721/139‑1 (MQ=255)               cGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAtt                                                                                                                      >  1:192962/1‑139 (MQ=255)               cGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAtt                                                                                                                      >  2:856784/1‑139 (MQ=255)                  ccGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGcc                                                                                                                   >  1:443258/1‑139 (MQ=255)                     cgAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGt                                                                                                                <  2:809873/139‑1 (MQ=255)                     cgAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGt                                                                                                                <  1:820927/139‑1 (MQ=255)                      gAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTa                                                                                                               <  1:4297/139‑1 (MQ=255)                                    cttcATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGAc                                                                                                 >  1:28213/1‑139 (MQ=255)                                        atcaAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCa                                                                                             <  1:852655/139‑1 (MQ=255)                                            agcagcTTTTTCGTAACGGAACCACCACGGCGGTGGTGAATGGCACTGTTCATCCGCAGTCTGTTGATGGGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGGTTGCCGGTAGGGTGATGATGGACCGCAACGc                                                                                         <  1:1090675/139‑1 (MQ=255)                                            agcagcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGc                                                                                         <  2:934959/139‑1 (MQ=255)                                             gcagcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCa                                                                                        >  2:769957/1‑139 (MQ=255)                                               agcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCAcc                                                                                      >  2:105761/1‑139 (MQ=255)                                                 cTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                                  <  2:697921/93‑1 (MQ=255)                                                 cTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                                  >  1:697921/1‑93 (MQ=255)                                                       cGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATc                                                                              >  1:814485/1‑139 (MQ=255)                                                            cGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTc                                                                         >  2:392568/1‑139 (MQ=255)                                                              gAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGa                                                                       <  1:341344/139‑1 (MQ=255)                                                                         gcgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACACACCGGATTATCTGCTCGACACTGCCGaaa                                                            >  2:1104778/1‑139 (MQ=255)                                                                          cgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAag                                                           <  2:547609/139‑1 (MQ=255)                                                                                 gtTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTAt                                                    <  1:856784/139‑1 (MQ=255)                                                                                      gCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACca                                               <  1:306694/139‑1 (MQ=255)                                                                                          tGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagc                                           >  1:246183/1‑139 (MQ=255)                                                                                            ttCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGTGTATGATTGCCGGTAAGGTGACGATGGACCGCAACGCACCGGATTATCTACTCGACACTGCCGAAAGCAGCTATCACcaaagcaa                                         >  2:669912/1‑139 (MQ=255)                                                                                            ttCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaa                                         >  2:418957/1‑139 (MQ=255)                                                                                                      aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCCGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATc                               >  1:243772/1‑139 (MQ=255)                                                                                                      aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCAAATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGCACGGATc                               >  1:274666/1‑139 (MQ=255)                                                                                                           ttGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACg                          >  1:166099/1‑139 (MQ=255)                                                                                                           ttGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACg                          <  1:392568/139‑1 (MQ=255)                                                                                                            tGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCa                                                                                             >  2:360982/1‑71 (MQ=255)                                                                                                            tGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCa                                                                                             <  1:360982/71‑1 (MQ=255)                                                                                                               tGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTgg                      <  1:893467/139‑1 (MQ=255)                                                                                                                 cgcTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGca                    >  1:360630/1‑139 (MQ=255)                                                                                                                         gAAGCCGCCAGTCACATCAATATGTGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATg            <  1:669912/139‑1 (MQ=255)                                                                                                                         gAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATg            <  2:548571/139‑1 (MQ=255)                                                                                                                          aaGCCGCCAGTCACATCAATATGCGTATGGTTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGGCACTGCCCAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATgg           <  2:360630/139‑1 (MQ=255)                                                                                                                           aGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCCGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgt          <  2:243772/139‑1 (MQ=255)                                                                                                                             ccgccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcg        <  1:769957/139‑1 (MQ=255)                                                                                                                              cgccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcgt       >  1:411823/1‑139 (MQ=255)                                                                                                                              cgccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGAGGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcgt       <  2:984110/139‑1 (MQ=255)                                                                                                                               gccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcgtc      <  1:86766/139‑1 (MQ=255)                                                                                                                               gccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcgtc      <  1:509916/139‑1 (MQ=255)                                                                                                                                ccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCt     >  1:689537/1‑139 (MQ=255)                                                                                                                                   gTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCt  <  2:166099/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       GAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCT  >  NC_000913/3026105‑3026373

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