Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,490,865	T→C	intergenic (‑152/‑37)	gapC ← / → cybB	glyceraldehyde‑3‑phosphate dehydrogenase/cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,490,865	0	T	C	100.0%	181.6 / NA	49	intergenic (‑152/‑37)	gapC/cybB	glyceraldehyde‑3‑phosphate dehydrogenase/cytochrome b561
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/19);  total (30/19)

TTCGGATGAAAATTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTC  >  NC_000913/1490735‑1491000                                                                                                                                   |                                                                                                                                        ttCGGATGAAAATTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACg                                                                                                                                >  1:358734/1‑140 (MQ=255)  tCGGATGAAAATTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACgg                                                                                                                               >  1:175384/1‑140 (MQ=255)  tCGGATGAAAATTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACgg                                                                                                                               >  1:762945/1‑140 (MQ=255)           aaTTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATg                                                                                                                      <  2:358734/140‑1 (MQ=255)                 tCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAg                                                                                                                <  1:10271/140‑1 (MQ=255)                     ttCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGaaaa                                                                                                                           >  1:1018434/1‑125 (MQ=255)                     ttCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGaaaa                                                                                                                           <  2:1018434/125‑1 (MQ=255)                      tCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGa                                                                                                           >  2:256360/1‑140 (MQ=255)                      tCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGa                                                                                                           >  2:1132423/1‑140 (MQ=255)                       ccAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGaa                                                                                                          <  2:762945/140‑1 (MQ=255)                       ccAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGaa                                                                                                          <  1:1132423/140‑1 (MQ=255)                       ccAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGaa                                                                                                          <  1:256360/140‑1 (MQ=255)                       ccAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGaa                                                                                                          >  1:713710/1‑140 (MQ=255)                               aGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATg                                                                                                  >  2:983773/1‑140 (MQ=255)                                     ggATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAaaat                                                                                            >  2:604097/1‑140 (MQ=255)                                      gATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAaaata                                                                                           >  1:42377/1‑140 (MQ=255)                                      gATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAaaata                                                                                           <  1:983773/140‑1 (MQ=255)                                      gATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAaaata                                                                                           >  1:1021531/1‑140 (MQ=255)                                          gCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAaataaata                                                                                       <  2:754380/140‑1 (MQ=255)                                             ggATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTGACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTc                                                                                    <  1:1164844/140‑1 (MQ=255)                                                tGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAg                                                                                 >  1:1008935/1‑140 (MQ=255)                                                  ggCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGt                                                                               <  1:758830/140‑1 (MQ=255)                                                     aaTCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGATATGTTATGGAAAATAGATAATCAAGGATAc                                                                            >  2:758830/1‑140 (MQ=255)                                                        ctcCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACaaa                                                                         <  1:604097/140‑1 (MQ=255)                                                         tcCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAAt                                                                        <  2:1173136/140‑1 (MQ=255)                                                             ttCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGc                                                                    >  1:71377/1‑140 (MQ=255)                                                                 cacaGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGc                                                                    >  1:380565/1‑136 (MQ=255)                                                                 cacaGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGc                                                                    <  2:380565/136‑1 (MQ=255)                                                                   caGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCAc                                                              >  2:758912/1‑140 (MQ=255)                                                                     gTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCActg                                                            >  1:337120/1‑140 (MQ=255)                                                                            gtaacgtaCTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAATTATTCAAGGTTACAAATCAGCATTCACTGGCTGGTc                                                     >  2:504036/1‑140 (MQ=255)                                                                                      gAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTgg                                           >  2:664671/1‑140 (MQ=255)                                                                                      gAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTgg                                           >  1:711924/1‑140 (MQ=255)                                                                                       aaaaCGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGt                                          <  2:697565/140‑1 (MQ=255)                                                                                       aaaaCGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGt                                          >  1:697565/1‑140 (MQ=255)                                                                                         aaCGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTa                                        <  2:970492/140‑1 (MQ=255)                                                                                         aaCGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTa                                        <  2:337120/140‑1 (MQ=255)                                                                                            gggTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCg                                     >  2:178696/1‑140 (MQ=255)                                                                                            gggTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCg                                     >  2:1214583/1‑140 (MQ=255)                                                                                               tGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAg                                  >  2:1167627/1‑140 (MQ=255)                                                                                                      tatTTAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTg                           >  1:1218372/1‑140 (MQ=255)                                                                                                        tttAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTgcg                         >  1:16552/1‑140 (MQ=255)                                                                                                        tttAATGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTgcg                         >  2:558652/1‑140 (MQ=255)                                                                                                            aTGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAAt                     >  2:706027/1‑140 (MQ=255)                                                                                                             tGAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATg                    >  1:251729/1‑140 (MQ=255)                                                                                                              gAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCa                                   <  1:517985/124‑1 (MQ=255)                                                                                                              gAAATTTTGAGAAAAGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCa                                   >  2:517985/1‑124 (MQ=255)                                                                                                                            aGCCCGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTttc     <  2:251729/140‑1 (MQ=255)                                                                                                                               ccGTCATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTttcttc  <  1:558652/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        TTCGGATGAAAATTGATCCTTTCCAAGCTTAGACCAGGATGGCGGGATGGGCAATCTCCATTCTCACAGTGAAACGTAACGTACTGAAAACGGGTGAACAATATTTAATGAAATTTTGAGAAAAGCCCGTTATGTTAACGGAAAATTATGTTAAAGCAGGAAATGTTATGGAAAATAAATATTCAAGGTTACAAATCAGCATTCACTGGCTGGTCTTTTTACTGGTTATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTC  >  NC_000913/1490735‑1491000

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