Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,291	A→G	V130V (GTA→GTG)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,291	0	A	G	100.0%	167.9 / NA	50	V130V (GTA→GTG)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/23);  total (27/23)

TTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACG  >  NC_000913/1491162‑1491425                                                                                                                                  |                                                                                                                                       tttGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTaa                                                                                                                              >  2:697761/1‑140 (MQ=255) tttGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTaa                                                                                                                              <  1:468299/140‑1 (MQ=255)       gATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGt                                                                                                                        >  2:582786/1‑140 (MQ=255)        aTTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTg                                                                                                                       <  1:582786/140‑1 (MQ=255)          ttATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGc                                                                                                                     >  1:635583/1‑140 (MQ=255)                         tGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGa                                                                                                      <  1:697761/140‑1 (MQ=255)                           cgcTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAAt                                                                                                    >  2:567708/1‑140 (MQ=255)                              tGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTg                                                                                                 <  2:127931/140‑1 (MQ=255)                                   gTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGata                                                                                            <  2:554802/140‑1 (MQ=255)                                      aTTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATAttt                                                                                         >  1:543698/1‑140 (MQ=255)                                             tggtGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATAATTTGTCATc                                                                                  <  2:94970/140‑1 (MQ=255)                                              ggtGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCg                                                                                 >  2:205204/1‑140 (MQ=255)                                              ggtGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCg                                                                                 >  1:573450/1‑140 (MQ=255)                                              ggtGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCg                                                                                 >  1:667999/1‑140 (MQ=255)                                              ggtGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGCGGCTTGAAACGCTGGCAAATCGGGGATATTTTGTCATCg                                                                                 >  2:318436/1‑140 (MQ=255)                                                tgatgaTGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCggg                                                                               >  1:860315/1‑140 (MQ=255)                                                    gatgTATAACCGGGGCAACCCGTGGTTTGCGTTTGGGTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGACACGCTGGCGAAACTCGGATATTTTGTCATCGGGTTGc                                                                           >  1:356430/1‑140 (MQ=255)                                                          taACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTg                                                                     >  2:566426/1‑140 (MQ=255)                                                             ccGGGGCAACCCGTGGTTTGCGCTTGGTTTGTCGATGCCTTACGCTTCGGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCgg                                                                  <  2:356430/140‑1 (MQ=255)                                                              cGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGc                                                                 <  1:567708/140‑1 (MQ=255)                                                               ggggCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCt                                                                >  1:86017/1‑140 (MQ=255)                                                                gggCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTg                                                               >  2:139561/1‑140 (MQ=255)                                                                gggCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTg                                                               <  1:139561/140‑1 (MQ=255)                                                                  gCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCa                                                             >  1:582079/1‑140 (MQ=255)                                                                  gCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCa                                                             <  2:582079/140‑1 (MQ=255)                                                                  gCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCa                                                             <  2:833812/140‑1 (MQ=255)                                                                  gCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCa                                                             >  1:463602/1‑140 (MQ=255)                                                                     aCCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTg                                                          >  1:747698/1‑140 (MQ=255)                                                                      cccGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTgg                                                         <  2:635583/140‑1 (MQ=255)                                                                      cccGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTgg                                                         <  2:543698/140‑1 (MQ=255)                                                                        cGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGca                                                       <  2:86017/140‑1 (MQ=255)                                                                            gTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCAcacc                                                   <  1:566426/140‑1 (MQ=255)                                                                               tGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACt                                                <  2:860315/140‑1 (MQ=255)                                                                                   tttggtttgACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTAttt                                            >  2:57615/1‑140 (MQ=255)                                                                                      ggtttgACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTg                                         >  2:452871/1‑140 (MQ=255)                                                                                       gtttgACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTgg                                        <  2:449072/140‑1 (MQ=255)                                                                                       gtttgACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTgg                                        <  2:167966/140‑1 (MQ=255)                                                                                       gtttgACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTgg                                        <  2:590585/140‑1 (MQ=255)                                                                                            aCGATGCCTTACGCTTCAGAGGTCAATTTAGAACGGGCGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGa                                   <  1:279114/140‑1 (MQ=255)                                                                                            aCGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGa                                   <  1:205204/140‑1 (MQ=255)                                                                                                 gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACAGGCACAGCACTATGTCTGGAAGGacaaca                              >  1:267283/1‑140 (MQ=255)                                                                                                 gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGAATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGacaaca                              <  2:463602/140‑1 (MQ=255)                                                                                                            cAGAGGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGc                   >  1:52674/1‑140 (MQ=255)                                                                                                             agagGCCAATTTCGAACGGGTGGATCGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCa                  <  1:318436/140‑1 (MQ=255)                                                                                                               agGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCatg                >  2:465438/1‑140 (MQ=255)                                                                                                               agGCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCCGCACTGGCGCACGACTGTTTCTGGAGGGACGACACACTTCTACGCatg                >  2:248868/1‑140 (MQ=255)                                                                                                                 gCCAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCatgat              >  2:125324/1‑140 (MQ=255)                                                                                                                  ccAATTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCatgatg             <  1:125324/140‑1 (MQ=255)                                                                                                                    aaTTTCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGcc           >  1:545031/1‑140 (MQ=255)                                                                                                                      tttCGAACGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCcgc         >  1:613200/1‑140 (MQ=255)                                                                                                                           aaCGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACTACAAACTTCTACGCATGATGCCGCGTaaa    >  1:328621/1‑140 (MQ=255)                                                                                                                             cGGGTGGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACg  <  1:452871/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                       TTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACG  >  NC_000913/1491162‑1491425

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: ‑