Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,331,358	A→G	D104G (GAC→GGC)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,331,358	0	A	G	97.7%	143.5 / NA	44	D104G (GAC→GGC)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (26/17);  minor base C (1/0);  total (27/17)
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

CCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTAAAAACGCGATCCGCC  >  NC_000913/1331222‑1331492                                                                                                                                         |                                                                                                                                       ccTAAAAAGGATGAACGCGGCGGTCTCGTCAGCCGCATGTGGGTTGACTCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCca                                                                                                                                     <  2:344325/140‑1 (MQ=255)                         ctcGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCGGGCAACCTCTATCAGGCAACCGACCTTGa                                                                                                             >  1:325882/1‑140 (MQ=255)                              gacaggTATGGGGGTAGACCCGTGGCACACTTGGGAGCCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACcgcg                                                                                                        <  1:247133/135‑1 (MQ=255)                              cAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACcgcg                                                                                                        <  1:712634/140‑1 (MQ=255)                                     aTGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGa                                                                                                 <  2:809736/140‑1 (MQ=255)                                     aTGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGa                                                                                                 <  1:710383/140‑1 (MQ=255)                                      tGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTTAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGaa                                                                                                >  1:606818/1‑140 (MQ=255)                                           gTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCAt                                                                                           >  1:540138/1‑140 (MQ=255)                                             tGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGa                                                                                                 >  2:662243/1‑132 (MQ=255)                                             tGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGa                                                                                                 <  1:662243/132‑1 (MQ=255)                                              gACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTTTCTAGCATCCGACCTTTACCGAGAAGGGGAAGccctgtc                                                                                        >  2:247133/1‑135 (MQ=255)                                              gACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGc                                                                                        >  2:94199/1‑140 (MQ=255)                                              gACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGc                                                                                        >  1:101562/1‑140 (MQ=255)                                              gACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGc                                                                                        >  1:105727/1‑140 (MQ=255)                                                  cGTGGCACAATTGGGAGGCGCACTATGAAGGGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACATCTGGCTGAAAAAGCCGGCCACATCTATCTCGCACCCGACCTTGACCGCGAAGGGGAAGCCATTGCATgg                                                                                    <  2:193671/140‑1 (MQ=255)                                                              gggAGGCGCACTATGGAGTGTTGCCTGGTTAAGAGAATGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGaa                                                                        <  1:324334/140‑1 (MQ=255)                                                              gggAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGaa                                                                        <  1:454148/140‑1 (MQ=255)                                                              gggAGGCGCACTATGAAGTGTTCCCAGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGaa                                                                        <  1:342642/140‑1 (MQ=255)                                                                    cgcACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGGGAATGGGAAGCCATTGGATGGAACCTGCGGGAAGTGAtt                                                                  >  2:262419/1‑140 (MQ=255)                                                                      cACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCATGACCGCGGAGGAGAAGCCATTGCATGGCACCTGCGGGAAGTGATTgg                                                                >  2:324334/1‑140 (MQ=255)                                                                      cACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCATGACCGCGAAGGGGAATCCATTGCAGGGCATCTGCGGGAAGTGATTgg                                                                >  2:342642/1‑140 (MQ=255)                                                                      cACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCAAATCTATCTCGCAACCGACCTTGACCGCGAAGGAGAATCCATTGCACGGCAACTGAGGGAAGTGATTgg                                                                >  2:354644/1‑140 (MQ=255)                                                                             aaGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGCCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGCAGTGATTGGGGGtgat                                                         >  1:497673/1‑140 (MQ=255)                                                                                     ccTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATgcgcg                                                 >  1:547831/1‑140 (MQ=255)                                                                                     ccTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATgcgcg                                                 >  1:735679/1‑140 (MQ=255)                                                                                     ccTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATgcg                                                   >  1:521633/1‑138 (MQ=255)                                                                                     ccTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATgcg                                                   <  2:521633/138‑1 (MQ=255)                                                                                     ccTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGCCACCTGCGGGAACTGATACGGGGTGATGATgcgcg                                                 >  2:375560/1‑140 (MQ=255)                                                                                      cTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATgcgcgc                                                <  2:101562/140‑1 (MQ=255)                                                                                          tAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCtata                                            <  1:864748/140‑1 (MQ=255)                                                                                          tAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCtata                                            >  1:677747/1‑140 (MQ=255)                                                                                               agaAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGa                                       >  2:795697/1‑140 (MQ=255)                                                                                                gaAGGTCGTCTCTGAACTGAAATAACTGGTTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAg                                      <  1:262419/140‑1 (MQ=255)                                                                                                  aGGTCGTCTCTGAACTGAAACACCTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAgtg                                    <  1:375560/140‑1 (MQ=255)                                                                                                    gtcgtcTCTGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAgtggt                                  >  1:66539/1‑140 (MQ=255)                                                                                                         ctctGAACTGAAACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTa                             >  1:96445/1‑140 (MQ=255)                                                                                                          tctGAACTGAAACAAATGGCTGTAAAAGCTGGCCACACCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTaa                            <  1:354644/140‑1 (MQ=255)                                                                                                             gaactgaaACAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGa                         >  1:435201/1‑140 (MQ=255)                                                                                                                ctgaaACAACTGGCTGATAAAGTCGGCCACACCTACCTCCCAACCGACATTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAAt                      <  2:258366/140‑1 (MQ=255)                                                                                                                   aaaCAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTAc                   >  2:88572/1‑140 (MQ=255)                                                                                                                   aaaCAACTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATCATGCGCGGTATAGCGGAGTGGCGTTTAACGAAATTAc                   >  1:312450/1‑140 (MQ=255)                                                                                                                        acTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTaaaa              >  1:763029/1‑140 (MQ=255)                                                                                                                         cTGGCTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGTGCTATAGTCGAGTGATGTTTACCGAAATTACTaaaaa             >  1:258366/1‑140 (MQ=255)                                                                                                                             cTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTAAAAAcgcg         >  1:6682/1‑140 (MQ=255)                                                                                                                             cTGAAAAAGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTAAAAAcgcg         >  1:802071/1‑140 (MQ=255)                                                                                                                                    aGCCGGCCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTAAAAACGCGATCCGcc  <  2:435201/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCTTGACCGCGAAGGGGAAGCCATTGCATGGCACCTGCGGGAAGTGATTGGGGGTGATGATGCGCGCTATAGCCGAGTGGTGTTTAACGAAATTACTAAAAACGCGATCCGCC  >  NC_000913/1331222‑1331492

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