Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,204	T→C	G101G (GGT→GGC)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,204	0	T	C	98.2%	198.6 / NA	55	G101G (GGT→GGC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (27/27);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTG  >  NC_000913/1491069‑1491340                                                                                                                                        |                                                                                                                                         ggTCGTTCGTTTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggcttgg                                                                                                                                      <  2:231245/140‑1 (MQ=255) ggTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggtttgg                                                                                                                                      <  1:1613864/140‑1 (MQ=255) ggTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggcttgg                                                                                                                                      <  1:1011770/140‑1 (MQ=255) ggTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggcttgg                                                                                                                                      >  2:1388341/1‑140 (MQ=255)    ctttcgtCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtga                                                                                                                                   <  2:99914/138‑1 (MQ=255)    cgttcgtCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtga                                                                                                                                   <  2:341581/140‑1 (MQ=255)    cgttcgtCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtga                                                                                                                                   <  1:56606/140‑1 (MQ=255)    cgttcgtCTTTTGTAAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtga                                                                                                                                   <  2:463427/140‑1 (MQ=255)     gttcttCTTTTGTAAGGGCTGAAATCCCCACCGCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtgat                                                                                                                                  <  1:1129495/140‑1 (MQ=255)     gttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtgat                                                                                                                                  <  1:110012/140‑1 (MQ=255)      ttcgtCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGtgatg                                                                                                                                 <  1:113420/140‑1 (MQ=255)        cgtCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTgatgat                                                                                                                               <  1:217497/140‑1 (MQ=255)         gtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGatgatg                                                                                                                              >  1:1606657/1‑140 (MQ=255)          tCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGt                                                                                                                             <  2:1589767/140‑1 (MQ=255)              ttGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATaa                                                                                                                         >  1:204105/1‑140 (MQ=255)                   aaGGCTAAAAAACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCggg                                                                                                                    <  1:896455/140‑1 (MQ=255)                    aGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCgggg                                                                                                                   <  1:1544910/140‑1 (MQ=255)                    aGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCgggg                                                                                                                   <  1:1540283/140‑1 (MQ=255)                     ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGc                                                                                                                  >  2:126919/1‑140 (MQ=255)                         gAAATACCCAACCCCACCGATTATACCTAAGCCAAAACCGATGATGACGGGACTTACGCATTTGGGACATTTGGAGATTTATCTCCTGTTTATTGCTCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGCGGCAAcc                                                                                                              >  1:1562321/1‑140 (MQ=255)                              ccccACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTgg                                                                                                         <  2:483/139‑1 (MQ=255)                                ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                                                                       >  2:461045/1‑140 (MQ=255)                                ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                                                                       >  1:63145/1‑140 (MQ=255)                                ccAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGtt                                                                                                       >  1:647061/1‑140 (MQ=255)                                     cccGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgt                                                                                                  <  2:1638365/140‑1 (MQ=255)                                       cgccgATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgttt                                                                                                >  1:1176371/1‑140 (MQ=255)                                         ccgATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCgtttgg                                                                                              >  2:1418612/1‑140 (MQ=255)                                                 aCCTAAGCCAAAACCGATGATTACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGa                                                                                      >  2:205696/1‑140 (MQ=255)                                                    tAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGc                                                                                   >  2:1472086/1‑140 (MQ=255)                                                     aaGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGcc                                                                                  <  2:1415995/140‑1 (MQ=255)                                                      aGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCt                                                                                 >  1:1634146/1‑140 (MQ=255)                                                          aaaaCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACg                                                                             <  1:126172/140‑1 (MQ=255)                                                           aaaCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGc                                                                            >  1:1357912/1‑140 (MQ=255)                                                               cGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCa                                                                        <  1:780829/140‑1 (MQ=255)                                                                 atgatgACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCaga                                                                      >  2:626960/1‑140 (MQ=255)                                                                 atgatgACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCaga                                                                      >  1:1360066/1‑140 (MQ=255)                                                                  tgatgaCGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCagag                                                                     >  2:599647/1‑140 (MQ=255)                                                                           gACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTc                                                            >  1:1027421/1‑140 (MQ=255)                                                                              tGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGaa                                                         >  2:1183714/1‑140 (MQ=255)                                                                              tGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGaa                                                         <  2:1176371/140‑1 (MQ=255)                                                                                     tttGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAg                                                  >  1:436406/1‑140 (MQ=255)                                                                                      ttGGGACATTTGGTGATTTATCTTCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGa                                                 <  2:1633466/140‑1 (MQ=255)                                                                                      ttGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGa                                                 <  1:740923/140‑1 (MQ=255)                                                                                      ttGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGa                                                 <  2:92913/140‑1 (MQ=255)                                                                                        gggACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATa                                               >  2:788396/1‑140 (MQ=255)                                                                                               ttGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTaaa                                        >  2:791735/1‑140 (MQ=255)                                                                                               ttGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTaaa                                        >  2:407693/1‑140 (MQ=255)                                                                                                 ggTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGt                                      <  1:126919/140‑1 (MQ=255)                                                                                                  gTGATTTATCTGCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCATCCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTc                                     <  2:604474/140‑1 (MQ=255)                                                                                                        tATCTCCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACTCTTCAGAGGCCAATTTCGACCGGGTAGATAGCTTAAAGTCGTGGCa                               <  2:1322160/140‑1 (MQ=255)                                                                                                            tcCTGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGaa                           >  2:1096614/1‑140 (MQ=255)                                                                                                               tGTTTATTGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACg                        <  2:53137/140‑1 (MQ=255)                                                                                                                     ttGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACAATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCg                  <  1:1388341/140‑1 (MQ=255)                                                                                                                      tGCGCTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGa                 <  2:1634146/140‑1 (MQ=255)                                                                                                                        cgcTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAAt               >  2:1337889/1‑140 (MQ=255)                                                                                                                        cgcTGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAAt               >  1:444572/1‑140 (MQ=255)                                                                                                                           tGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTg            >  1:1432319/1‑140 (MQ=255)                                                                                                                           tGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTg            >  1:38699/1‑140 (MQ=255)                                                                                                                           tGCCAGTGATTGGCTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGCAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTg            >  1:1236712/1‑140 (MQ=255)                                                                                                                                  gATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATAtt     <  1:205696/140‑1 (MQ=255)                                                                                                                                     tggcttggtGCTGTTGTATAACCGGGGCATCCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTg  <  2:1059747/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCAACCCGTGGTTTGCGTTTGGTTTGACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTG  >  NC_000913/1491069‑1491340

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