Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,129	A→G	K76K (AAA→AAG)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,129	0	A	G	100.0%	175.6 / NA	49	K76K (AAA→AAG)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/23);  total (26/23)

TTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACG  >  NC_000913/1490995‑1491265                                                                                                                                       |                                                                                                                                         ttcttcCCACGTAGCGATCGGCCACTCATAAACATTATTCATGTTTCGTGGGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCgat                                                                                                                                     <  2:1057488/140‑1 (MQ=255)       ccACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGAc                                                                                                                               <  2:79157/140‑1 (MQ=255)          cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACggg                                                                                                                            >  2:1065381/1‑140 (MQ=255)          cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACggg                                                                                                                            >  1:317966/1‑140 (MQ=255)          cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACggg                                                                                                                            >  1:550768/1‑140 (MQ=255)                 aTCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGgcg                                                                                                                     <  1:457352/140‑1 (MQ=255)                   cGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCa                                                                                                                   >  2:773114/1‑140 (MQ=255)                       cACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTg                                                                                                               >  1:36779/1‑140 (MQ=255)                         cTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCC‑CCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGa                                                                                                            >  1:288893/1‑140 (MQ=255)                                      ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                                                                >  1:859835/1‑140 (MQ=255)                                      ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                                                                >  1:722652/1‑140 (MQ=255)                                                    ggCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTAt                                                                                  >  1:888314/1‑140 (MQ=255)                                                        tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgcg                                                                              <  2:888314/140‑1 (MQ=255)                                                            aaTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                                                          <  2:859835/140‑1 (MQ=255)                                                             aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                                         >  2:949066/1‑140 (MQ=255)                                                              ttCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                                                        <  1:1065381/140‑1 (MQ=255)                                                                   gTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                                                   <  2:722652/140‑1 (MQ=255)                                                                   gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                                                   <  2:60886/140‑1 (MQ=255)                                                                   gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                                                   <  1:986396/140‑1 (MQ=255)                                                                   gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                                                                   >  1:518510/1‑140 (MQ=255)                                                                    tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAtt                                                                  >  1:433254/1‑140 (MQ=255)                                                                    tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAtt                                                                  >  2:3153/1‑140 (MQ=255)                                                                     gCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgcg                                                                              <  1:439094/127‑1 (MQ=255)                                                                     gCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgcg                                                                              >  2:439094/1‑127 (MQ=255)                                                                      cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttgg                                                                <  1:109352/140‑1 (MQ=255)                                                                           ggTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAGCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                                                           >  2:969926/1‑140 (MQ=255)                                                                              cgttcgtCTTTTGTTAAGGCTGTAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                        <  1:773114/140‑1 (MQ=255)                                                                              cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                        >  2:897419/1‑140 (MQ=255)                                                                              cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                        <  1:949066/140‑1 (MQ=255)                                                                               gtccgtCTTTTGTTAAGGCGGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                                       <  2:189271/137‑1 (MQ=255)                                                                                ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                                                      >  1:189271/1‑140 (MQ=255)                                                                                ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                                       >  2:565347/1‑139 (MQ=255)                                                                                ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                                       <  1:565347/139‑1 (MQ=255)                                                                                               ggCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc                                       <  1:910453/140‑1 (MQ=255)                                                                                               ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc                                       >  2:817444/1‑140 (MQ=255)                                                                                                gCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                                       >  2:701167/1‑123 (MQ=255)                                                                                                gCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgat                                                       <  1:701167/123‑1 (MQ=255)                                                                                                    aaaTACCCACCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAAccc                                  <  1:969926/140‑1 (MQ=255)                                                                                                    aaaTACCCAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAAccc                                  >  1:707207/1‑140 (MQ=255)                                                                                                         cccAACCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGt                             <  1:817444/140‑1 (MQ=255)                                                                                                          ccACCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                            <  2:504885/140‑1 (MQ=255)                                                                                                              ccccGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCg                        >  2:521724/1‑140 (MQ=255)                                                                                                                         ataCCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGAc             <  1:1058178/140‑1 (MQ=255)                                                                                                                             cTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATg         <  2:707207/140‑1 (MQ=255)                                                                                                                                aGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCt      <  1:106493/140‑1 (MQ=255)                                                                                                                                aGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCt      >  1:236323/1‑140 (MQ=255)                                                                                                                                aGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTGGCGTTTGGTTTGACGATGCCt      >  2:996042/1‑140 (MQ=255)                                                                                                                                 gCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCtt     >  1:489734/1‑140 (MQ=255)                                                                                                                                 gCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCtt     >  1:1071054/1‑140 (MQ=255)                                                                                                                                    aaagCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg  <  1:897419/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         TTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACG  >  NC_000913/1490995‑1491265

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