Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,973,572	A→G	10.6%	L585L (TTA→CTA)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,572	0	A	G	10.6%	79.5 / 3.1	39	L585L (TTA→CTA)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base G (2/2);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

CACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCT  >  NC_000913/1973432‑1973708                                                                                                                                             |                                                                                                                                         cACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGCGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACCTTGTAAGa                                                                                                                                        >  1:102980/1‑143 (MQ=255)    gCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATca                                                                                                                                     >  1:316117/1‑143 (MQ=255)             gCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGCTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTGGTAAGATCACCACAACTc                                                                                                                            >  2:230526/1‑143 (MQ=255)              ccAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTcc                                                                                                                           >  2:809476/1‑143 (MQ=255)               cAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACACCCTGTTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAGGATCACCACAATTccc                                                                                                                          >  2:844122/1‑143 (MQ=255)               cAAGAATGGTCGCAGAAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAGGATCACCACAATTccc                                                                                                                          >  2:770865/1‑143 (MQ=255)                    aTGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGt                                                                                                                     <  1:510124/143‑1 (MQ=255)                          gcagcaGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTCAGATCACCAGAATTCCCTGGGTGGCtta                                                                                                               >  1:279939/1‑142 (MQ=255)                            agcagAAATGCCGGCCACTCTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCgg                                                                                                             <  1:758524/143‑1 (MQ=255)                             gcagAAATGCGTGGCACTGGGCGATAGTTCCTTTCAAGGTCTTTAACTACAACCTGGTGTTGACCAATTAATTTATCCACCAGCAAGGCGGAGCGGCGACCGCCACTGTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGt                                                                                                            <  2:729160/143‑1 (MQ=255)                              cagAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGtt                                                                                                           <  2:676366/143‑1 (MQ=255)                                 aaaTGCCGGGGACTTTGCGATAGTTACTTTCAAGCTTTTTATCCACAACCTGGTGTTGACCAATCAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTc                                                                                                        <  2:736377/143‑1 (MQ=255)                                 aaaTGCCGGGGACTTTGAGATAGTTACTTGCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTc                                                                                                        <  2:836251/143‑1 (MQ=255)                                              tttGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGt                                                                                           <  1:700859/143‑1 (MQ=255)                                               ttGCGAGAGTTACTTTCAAGGTTTTTAACGAGAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGtt                                                                                          <  1:570464/143‑1 (MQ=255)                                                    aTAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGcaca                                                                                     >  1:282617/1‑143 (MQ=255)                                                        ttACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGAGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACttt                                                                                 >  1:779420/1‑143 (MQ=255)                                                         gaCTTTCATGGTTTTCAAGCACAACCTGGTTTGGACCAGTTAATTGATCCACCAGCAATGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTc                                                                                <  2:442688/142‑1 (MQ=255)                                                            tttCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATAGATCCACCAGCAAGGCGTAGCAGCGACCGCCACTATGTAAGATCACCACAATTCCCTGGGTGGCTTCAGGCGTAGCGGGCGGGAAGTTCAACACTTTCcac                                                                             >  2:742243/1‑143 (MQ=255)                                                                 atgTTTTTAAACATAACCTTGTGTTGACCGATTTATTGATCCACCATCAAGGCGTAGCGGCGACTGCCATTTTGTAAGATCTCCACATTTCCCTGGGTGGCTTCGGTTTTCGCGCCTGCGACGTTGAACACTTTCCACAGTTc                                                                        <  2:36883/141‑1 (MQ=255)                                                                 aGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACACTTCCCTGGGTGGCTTCGGTATTCGCGCCCGCGCCGTTGAACACTTCCCACGGTTc                                                                        >  2:572495/1‑143 (MQ=255)                                                                       ttACCCACGTCCTGGTGTTTGCCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGCGGCTTCTGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGAt                                                                  <  2:744629/143‑1 (MQ=255)                                                                       ttAACCCCAGCCTGGTGTTGACCGATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAACGTCCCCACAATTCCCTGGGTGGCGTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGAt                                                                  <  2:72850/143‑1 (MQ=255)                                                                                   tGGAGTTGACCAATTAATTGATCCACCAGCAAGCCGTCGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTc                                                      <  2:734453/143‑1 (MQ=255)                                                                                    ggTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAGGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCa                                                     <  1:770865/143‑1 (MQ=255)                                                                                    ggTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAGGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCa                                                     <  1:844122/143‑1 (MQ=255)                                                                                               attaattGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTc                                          <  2:488737/143‑1 (MQ=255)                                                                                                          caccaGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACCTCCAGCACCCGcc                               >  2:829433/1‑142 (MQ=255)                                                                                                                cAATGCGCAGCGGCGTCCGCCACTTTGTACGATCACCACAATCCCCTGGGTGGCTTCGGTTTGCGCGCCCGCTACGTTGAACACTTTCCACAGTTCGACCATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTcgccgc                         <  2:418692/143‑1 (MQ=255)                                                                                                                  aGGCGTAGCGGCGACCGCCACTTTGTAAGATTACCTTTATTTCTTTGGTGGCTTCGGTGTTCGCGCCCGCGACGATGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGccgccg                       <  1:102859/143‑1 (MQ=255)                                                                                                                  aGGCGTAGCGGCGACCGCCACTTTGTAAGATCGCCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTGCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGccgccg                       <  2:876545/143‑1 (MQ=255)                                                                                                                    gCGTAGCGGCGACCGCCACTCTTTTAGATCACCACCATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGCTGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGccgccgac                     <  2:279939/143‑3 (MQ=255)                                                                                                                     cgcggcggccACGCCCACGTTGTAAAATCACCACAATTCCCTGGGTGGCTTCGGTGTTCGCGCCCGCGACGTTGAACACCTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGcc                    <  1:714612/139‑1 (MQ=255)                                                                                                                           cgcgACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGtgga              <  2:909035/142‑1 (MQ=255)                                                                                                                            gcgACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGtggat             >  1:353699/1‑143 (MQ=255)                                                                                                                                ccGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGcag         >  1:400610/1‑143 (MQ=255)                                                                                                                                  gCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCTTGGACCCGCTCTCCGCCGGCCTGTGGGTGgagac       >  1:377717/1‑142 (MQ=255)                                                                                                                                   ccACTTTTTTAGATCACCACATTTCCCTGGGTGGCTTCGGTTTCCGAGCAATCGACGTTCAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCTCTCGCCGCCGGCCAGTGGATGGAGATc      <  2:316117/143‑1 (MQ=255)                                                                                                                                       tttGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGCTATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGCCt  >  2:595739/1‑143 (MQ=255)                                                                                                                                             |                                                                                                                                         CACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATTCACCCCGCACTTCCAGCACCCGCTCGCCGCCGGCCAGTGGATGGAGATCGGCT  >  NC_000913/1973432‑1973708

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