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	100.0%	162.8 / NA	44	H158H (CAT→CAC)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/20);  total (24/20)

GGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCT  >  NC_000913/3026107‑3026370                                                                                                                                     |                                                                                                                                    ggATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                               <  2:126581/139‑1 (MQ=255)   aTTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAatat                                                                                                                             >  2:915211/1‑139 (MQ=255)   aTTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAatat                                                                                                                             >  1:354329/1‑139 (MQ=255)            aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGTACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    <  2:354329/139‑1 (MQ=255)            aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGACGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    >  2:476349/1‑139 (MQ=255)            aCGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    >  2:400473/1‑139 (MQ=255)             cGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAtt                                                                                                                   >  2:814311/1‑139 (MQ=255)             cGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAtt                                                                                                                   >  1:165078/1‑139 (MQ=255)                  gcgAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCgg                                                                                                              <  1:842931/139‑1 (MQ=255)                       aTGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGTACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAgg                                                                                                         <  1:683643/139‑1 (MQ=255)                        tGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    <  2:339115/127‑1 (MQ=255)                        tGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    >  1:339115/1‑127 (MQ=255)                            ggCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    >  2:801622/1‑123 (MQ=255)                            ggCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAt                                                                                                                    <  1:801622/123‑1 (MQ=255)                                      atcaAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCa                                                                                          <  1:325002/139‑1 (MQ=255)                                          agcagcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGc                                                                                      <  2:584817/139‑1 (MQ=255)                                                              accaccACGGCGCTTGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCAGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGaca                                                                  <  2:806310/139‑1 (MQ=255)                                                              accaccACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGaca                                                                  <  1:915211/139‑1 (MQ=255)                                                                 accacGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTg                                                               >  2:500192/1‑139 (MQ=255)                                                                     cGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGa                                                           >  2:230162/1‑139 (MQ=255)                                                                        cgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAag                                                        <  1:541484/139‑1 (MQ=255)                                                                        cgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAag                                                        <  1:687888/139‑1 (MQ=255)                                                                        cgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAag                                                        <  2:165078/139‑1 (MQ=255)                                                                            ggTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAagcagc                                                    <  1:814311/139‑1 (MQ=255)                                                                             gtgtTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCt                                                   <  1:349149/139‑1 (MQ=255)                                                                                  tGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGG‑AAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCAcc                                             >  2:205099/1‑139 (MQ=255)                                                                                    gCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACca                                            <  1:407678/139‑1 (MQ=255)                                                                                     cACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaa                                           >  2:11571/1‑139 (MQ=255)                                                                                        tGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagc                                        >  1:904645/1‑139 (MQ=255)                                                                                          ttCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaa                                      >  2:152726/1‑139 (MQ=255)                                                                                              tCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGaa                                  >  2:92664/1‑139 (MQ=255)                                                                                                    aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGt                                                                                                             >  2:144452/1‑58 (MQ=255)                                                                                                    aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGt                                                                                                             <  1:144452/58‑1 (MQ=255)                                                                                                    aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATc                            >  1:303578/1‑139 (MQ=255)                                                                                                    aTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGAATGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATc                            >  1:843921/1‑139 (MQ=255)                                                                                                       tGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                               >  1:911556/1‑37 (MQ=255)                                                                                                       tGTTGATGCGCTGTTTGAAGCCGCCAGTCACATCAat                                                                                                                               <  2:911556/37‑1 (MQ=255)                                                                                                                 cTGTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGcaca               >  2:544906/1‑139 (MQ=255)                                                                                                                   gTTTGAAGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATAGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCAAACGCTGGCACaaa             >  1:1141272/1‑139 (MQ=255)                                                                                                                         aGCCGTCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgt       >  2:432994/1‑139 (MQ=255)                                                                                                                         aGCCGCCAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgt       <  2:509372/139‑1 (MQ=255)                                                                                                                           ccgtcAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcg     <  1:432994/139‑1 (MQ=255)                                                                                                                              ccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCt  >  1:275765/1‑139 (MQ=255)                                                                                                                              ccAGTCACATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCt  >  1:1058260/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                    GGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCT  >  NC_000913/3026107‑3026370

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