Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,791,313	Δ1 bp	coding (33/1026 nt)	tdh ←	threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,791,311	0	C	.	100.0%	187.5 / NA	44	coding (35/1026 nt)	tdh	threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (18/26);  total (18/26)

AGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATAC  >  NC_000913/3791234‑3791386                                                                              |                                                                            aGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGcctc                                                                          >  1:678468/1‑79 (MQ=255) aGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGcctc                                                                          >  1:678478/1‑79 (MQ=255) aGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGcctc                                                                          <  2:678399/81‑3 (MQ=255)    tGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTc                                                                      >  1:678456/1‑81 (MQ=255)    tGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTc                                                                      <  1:678401/81‑1 (MQ=255)       cTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGc                                                                   <  1:678402/81‑1 (MQ=255)          ttttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCttt                                                                >  2:678471/1‑81 (MQ=255)          ttttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCttt                                                                >  1:678467/1‑81 (MQ=255)          ttttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCt                                                                  <  1:678403/79‑1 (MQ=255)            ttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCttt                                                                >  1:678458/1‑79 (MQ=255)            ttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTc                                                               <  2:678405/80‑1 (MQ=255)            ttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTc                                                               <  2:678404/80‑1 (MQ=255)            ttACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTc                                                               <  2:678403/80‑1 (MQ=255)             tACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAg                                                             >  2:678476/1‑81 (MQ=255)             tACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAg                                                             >  1:678459/1‑81 (MQ=255)                ggATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGttt                                                          <  1:678408/81‑1 (MQ=255)                ggATTTTAATCAGCAGATCGTTATG‑CCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTg                                                         <  1:678407/81‑1 (MQ=255)                 gATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTg                                                         >  2:678484/1‑81 (MQ=255)                     ttAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGAt                                                      >  1:678462/1‑80 (MQ=255)                      tAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATa                                                     <  2:678408/80‑1 (MQ=255)                        aTCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAAc                                                   <  2:678409/80‑1 (MQ=255)                        aTCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACg                                                  >  1:678473/1‑81 (MQ=255)                         tCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGc                                                 >  2:678464/1‑81 (MQ=255)                           agcagATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCt                                                <  2:678410/80‑1 (MQ=255)                            gcagATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCttt                                              <  2:678411/81‑1 (MQ=255)                             cagATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTc                                             <  2:678412/81‑1 (MQ=255)                              agATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCa                                            >  1:678471/1‑81 (MQ=255)                                  cGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCAtct                                         >  2:678468/1‑80 (MQ=255)                                    ttATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTcac                                      <  1:678414/81‑1 (MQ=255)                                         cccGAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATcc                                  <  2:678414/80‑1 (MQ=255)                                            gAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAg                              <  1:678416/81‑1 (MQ=255)                                            gAGTTCCGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAg                              >  1:678482/1‑81 (MQ=255)                                                 ccGGTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGa                          <  1:678417/80‑1 (MQ=255)                                                    gTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAAc                      <  1:678419/81‑1 (MQ=255)                                                    gTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAAc                      <  2:678417/81‑1 (MQ=255)                                                    gTACAGGAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAAc                      <  1:678420/81‑1 (MQ=255)                                                         ggAACGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCa                 <  2:678420/81‑1 (MQ=255)                                                           aaCGTCGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAgt               >  1:678474/1‑81 (MQ=255)                                                           aaCGTCGGTCATCCAGATG‑CCTCTTCCGCTTACAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAgt               <  2:678421/81‑1 (MQ=255)                                                                cGGTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAgttgtt           >  2:678469/1‑80 (MQ=255)                                                                 ggTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGttt          <  2:678422/80‑1 (MQ=255)                                                                 ggTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGttt          <  2:678423/80‑1 (MQ=255)                                                                  gTCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTAc        >  1:678479/1‑81 (MQ=255)                                                                   tCATCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTAc        >  2:678475/1‑80 (MQ=255)                                                                     aTCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCaa     <  1:678425/81‑1 (MQ=255)                                                                      tCCAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAAt    <  1:678426/81‑1 (MQ=255)                                                                        cAGATG‑CCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATAc  >  2:678479/1‑81 (MQ=255)                                                                              |                                                                            AGATGGCTGTTTTACGGATTTTAATCAGCAGATCGTTATGCCCGAGTTCCGGTACAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATAC  >  NC_000913/3791234‑3791386

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