Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	772,209	C→A	T251N (ACC→AAC)	cydA →	cytochrome d terminal oxidase, subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	772,209	0	C	A	100.0%	210.7 / NA	60	T251N (ACC→AAC)	cydA	cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (37/23);  total (37/23)

CTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAC  >  NC_000913/772114‑772303                                                                                                |                                                                                               cTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACaaa                                                                                            >  1:444307/1‑100 (MQ=255)  tttGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGTAATGGGCGACGTGCAGAAaaacaaac                                                                                           <  1:964941/100‑1 (MQ=255)  tttGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAaaacaaac                                                                                           <  1:170717/100‑1 (MQ=255)    tGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTg                                                                                         >  1:716180/1‑100 (MQ=255)    tGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTg                                                                                         <  2:334174/100‑1 (MQ=255)     gCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTgg                                                                                        >  1:567321/1‑100 (MQ=255)        aTCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGgctg                                                                                     >  2:488472/1‑100 (MQ=255)        aTCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGgctg                                                                                     >  2:487678/1‑100 (MQ=255)           gCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTa                                                                                  >  2:744211/1‑100 (MQ=255)             tGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTAtt                                                                                >  1:99992/1‑100 (MQ=255)              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                               >  1:782792/1‑100 (MQ=255)              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                               >  1:71882/1‑100 (MQ=255)              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                               <  1:1292018/100‑1 (MQ=255)               ccAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGa                                                                              >  2:1382220/1‑100 (MQ=255)                cAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGaa                                                                             >  2:208598/1‑100 (MQ=255)                        gTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTg                                                                     >  1:1330645/1‑100 (MQ=255)                         tATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTgg                                                                    >  1:798996/1‑100 (MQ=255)                            ngctgctGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGaa                                                                 >  1:554969/2‑100 (MQ=255)                            ggCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGaa                                                                 >  1:989585/1‑100 (MQ=255)                               tgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACg                                                              <  1:298338/100‑1 (MQ=255)                                gctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGATAAACAAACTGGCTGCTATTGAAGCCGCGTGGGAAACGc                                                             >  2:1002680/1‑100 (MQ=255)                                gctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGc                                                             <  2:268683/100‑1 (MQ=255)                                  tgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCaa                                                           <  1:808815/100‑1 (MQ=255)                                      ctgtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCaa                                                           >  2:1298002/1‑96 (MQ=255)                                      ctgtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTg                                                       <  1:416132/100‑1 (MQ=255)                                         tctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTgcgc                                                    >  2:580101/1‑100 (MQ=255)                                         tctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTgcgc                                                    <  1:1042791/100‑1 (MQ=255)                                             ttATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGGGGGAAACGCAACCTGCGCctgc                                                >  1:1568515/1‑100 (MQ=255)                                              tATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgct                                               >  2:490350/1‑100 (MQ=255)                                               aTTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgctg                                              >  1:1421405/1‑100 (MQ=255)                                               aTTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgctg                                              >  2:53873/1‑100 (MQ=255)                                                ttGTTCTGGATGATGAATCCGGCTACGAAATGGGCGACGGGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgctg                                              <  2:630349/99‑1 (MQ=255)                                                   ttCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCtt                                          >  2:526302/1‑100 (MQ=255)                                                   ttCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCtt                                          <  2:1316335/100‑1 (MQ=255)                                                    tCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCttt                                         <  2:444306/100‑1 (MQ=255)                                                     cTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTa                                        >  1:419601/1‑100 (MQ=255)                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                       <  2:887594/100‑1 (MQ=255)                                                       gggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAct                                      >  1:1518062/1‑100 (MQ=255)                                                         gTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActct                                    >  1:563606/1‑100 (MQ=255)                                                             tgaATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActct                                    >  1:732413/1‑96 (MQ=255)                                                               aagccGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCgg                              <  1:559020/97‑1 (MQ=255)                                                               aaTCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCgg                              <  1:1323333/100‑1 (MQ=255)                                                                aTCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTCACTCTGTTCGGc                             <  1:1038510/100‑1 (MQ=255)                                                                 tCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGc                             <  1:1055474/99‑1 (MQ=255)                                                                  ccGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCAt                           <  1:199992/100‑1 (MQ=255)                                                                     gCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTcc                        >  1:1232288/1‑100 (MQ=255)                                                                       tACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTcccg                      <  1:1069822/100‑3 (MQ=255)                                                                        aCGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGa                     >  2:1439827/1‑100 (MQ=255)                                                                        aCGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGa                     >  2:709603/1‑100 (MQ=255)                                                                          gAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATc                   <  2:99992/100‑1 (MQ=255)                                                                          gAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATc                   <  2:716178/100‑1 (MQ=255)                                                                                 gCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAag             >  1:634144/1‑99 (MQ=255)                                                                                  cGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAagag           >  2:161241/1‑100 (MQ=255)                                                                                   gACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagg          >  2:1431663/1‑100 (MQ=255)                                                                                   gACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagg          >  2:1140024/1‑100 (MQ=255)                                                                                   gACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagg          <  1:561112/100‑1 (MQ=255)                                                                                   gACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagg          <  1:709605/100‑1 (MQ=255)                                                                                     cGTGCAGAAAAACAAACTGGCTGCTACTGCAGCCAAGGGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgaggag        <  2:267611/100‑1 (MQ=255)                                                                                      gTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGgaga       <  2:83489/100‑1 (MQ=255)                                                                                       tGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGgcga       >  2:84034/1‑99 (MQ=255)                                                                                        gCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACg     >  1:866258/1‑100 (MQ=255)                                                                                           gAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAc  >  2:318663/1‑100 (MQ=255)                                                                                                |                                                                                               CTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAC  >  NC_000913/772114‑772303

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