Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,159	A→G	G86G (GGA→GGG)	cybB →	cytochrome b561

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

ATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGA  >  NC_000913/1491022‑1491293                                                                                                                                          |                                                                                                                                       atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCa                                                                                                                                      >  1:232662/1‑140 (MQ=255)    aaCATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCAttt                                                                                                                                   <  2:887759/140‑1 (MQ=255)         gATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGa                                                                                                                              >  2:222292/1‑140 (MQ=255)           ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGAtt                                                                                                                            >  2:99385/1‑140 (MQ=255)           ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGAtt                                                                                                                            >  1:436323/1‑140 (MQ=255)                   tCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCt                                                                                                                    >  2:484489/1‑140 (MQ=255)                     cTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGt                                                                                                                  >  1:137592/1‑140 (MQ=255)                      tgtgGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCAGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGtt                                                                                                                 >  2:833824/1‑140 (MQ=255)                          gCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTAtt                                                                                                             <  1:484489/140‑1 (MQ=255)                             tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTgcg                                                                                                          <  1:110979/140‑1 (MQ=255)                             tctcAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTgcg                                                                                                          <  2:137592/140‑1 (MQ=255)                                  aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGcc                                                                                                     <  2:46426/140‑1 (MQ=255)                                  aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTGGGGGCATTTGGGGATTTATCTCCTGTTTATTGCGCTGcc                                                                                                     >  2:95845/1‑140 (MQ=255)                                   tcctcGGGCTGATGGTTGTTCGTTTTTTGTAAAGGGTAAAACCCCCAACCCCGCCGATTATTCCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                                                                                    <  1:235272/138‑1 (MQ=255)                                        gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                                                                                    <  1:109334/135‑1 (MQ=255)                                        gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCa                                                                                                    >  2:109334/1‑135 (MQ=255)                                         tGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAtt                                                                                              >  2:173809/1‑140 (MQ=255)                                            tGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggt                                                                                           <  2:136890/140‑1 (MQ=255)                                             gATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtt                                                                                          <  1:818456/140‑1 (MQ=255)                                              aTGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGGGATTTATGTCCTGTTTATTGCGCGGCCAGTGAttggttt                                                                                         >  2:235272/1‑140 (MQ=255)                                                   cgttcgttTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                                                    <  2:232662/140‑1 (MQ=255)                                                   cgttcgtCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                                                    <  1:361674/140‑1 (MQ=255)                                                   cgttcgtCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                                                                                    <  2:64381/140‑1 (MQ=255)                                                                    ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                                                                                  <  1:761121/125‑1 (MQ=255)                                                                    ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                                                                                  >  2:761121/1‑125 (MQ=255)                                                                    ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCCCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc                                                                   <  1:605950/140‑1 (MQ=255)                                                                               ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                        >  1:495336/1‑140 (MQ=255)                                                                               ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                        >  1:847137/1‑140 (MQ=255)                                                                               ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                        >  1:782776/1‑140 (MQ=255)                                                                                   ccccGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCg                                                    >  1:208015/1‑140 (MQ=255)                                                                                        ccgATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtttgg                                               >  2:113310/1‑140 (MQ=255)                                                                                         cgATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtttggt                                              >  1:329085/1‑140 (MQ=255)                                                                                             tataCCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGa                                          <  2:329085/140‑1 (MQ=255)                                                                                               taCCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACg                                        <  1:222292/140‑1 (MQ=255)                                                                                                  cTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATg                                     <  1:289339/140‑1 (MQ=255)                                                                                                     aGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGa                                          >  1:514625/1‑132 (MQ=255)                                                                                                     aGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGa                                          <  2:514625/132‑1 (MQ=255)                                                                                                                atggtgaCGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg                              >  1:596090/5‑133 (MQ=255)                                                                                                                atggtgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg                              <  2:596090/133‑1 (MQ=255)                                                                                                                atgatgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGcc                                   <  2:457609/128‑1 (MQ=255)                                                                                                                atgatgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGcc                                   >  1:457609/1‑128 (MQ=255)                                                                                                                atgatgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg                              <  2:683295/133‑1 (MQ=255)                                                                                                                atgatgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg                              >  1:683295/1‑133 (MQ=255)                                                                                                                atgatgACGGGACTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCaga                       >  1:584784/1‑140 (MQ=255)                                                                                                                           aCTGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCg            >  2:306013/1‑140 (MQ=255)                                                                                                                             tGGCGCATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGAGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGaa          >  1:457471/1‑140 (MQ=255)                                                                                                                                cgcATTTGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACgg       >  1:462752/1‑140 (MQ=255)                                                                                                                                    tttGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAg   >  2:912153/1‑140 (MQ=255)                                                                                                                                     ttGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGa  <  1:99385/140‑1 (MQ=255)                                                                                                                                     ttGGGGCATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGa  <  1:398207/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       ATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGA  >  NC_000913/1491022‑1491293

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