Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,110,515	C→G	10.2%	A317P (GCG→CCG)	PP_5355 ←	sodium/proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,110,515	0	C	G	10.2%	123.7 / 6.7	40	A317P (GCG→CCG)	PP_5355	sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (17/18);  new base G (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

ACCGGGTATAGGAAGTCGCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAATATGGCGCT  >  NC_002947/6110382‑6110641                                                                                                                                      |                                                                                                                               aCCGGGTATATGAAGTCGCTGGTCACCTGCATGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGccc                                                                                                                              <  1:127075/136‑1 (MQ=255) aCCGGGTATAGGAAGTCGCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGccc                                                                                                                              <  2:599324/136‑1 (MQ=255) aCCGGGTATAGGAAGTCGCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGccc                                                                                                                              <  2:114022/136‑1 (MQ=255)      gTATAGGAAGTCGCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATg                                                                                                                         <  2:440006/136‑1 (MQ=255)                  gCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGc                                                                                                             >  1:13404/1‑136 (MQ=255)                  gCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGc                                                                                                             <  1:321518/136‑1 (MQ=255)                  gCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGACCCCATGCCGCAAGACAGc                                                                                                             <  1:217745/136‑1 (MQ=255)                   cTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCTAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCa                                                                                                            >  1:395025/1‑136 (MQ=255)                   cTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCa                                                                                                            >  1:160166/1‑136 (MQ=255)                   cTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCa                                                                                                            >  1:125015/1‑136 (MQ=255)                       tCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCtt                                                                                                        <  2:13404/136‑1 (MQ=255)                            tGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCa                                                                                                   >  1:455610/1‑136 (MQ=255)                             gCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAg                                                                                                  <  2:141521/136‑1 (MQ=255)                                   tCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCACCATGCCGCAAGACAGCATCTTGCCCAGTaccac                                                                                            >  2:404310/1‑136 (MQ=255)                                        cccAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAg                                                                                                                  >  2:662533/1‑109 (MQ=255)                                        cccAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAg                                                                                                                  <  1:662533/109‑1 (MQ=255)                                        cccAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCa                                                                                       >  2:455473/1‑136 (MQ=255)                                         ccAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAg                                                                                      >  2:254394/1‑136 (MQ=255)                                              gcggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGGCCCCATGCCGCAAGACAGCATCTTGCCCAGTaccacc                                                                                           >  1:417152/1‑126 (MQ=255)                                              gcggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGGCCCCATGCCGCAAGACAGCATCTTGCCCAGTaccacc                                                                                           >  1:418702/1‑126 (MQ=255)                                              gcggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGGCCCCATGCCGCAAGACAGCATCTTGCCCAGTaccacc                                                                                           <  2:418702/126‑1 (MQ=255)                                              gcggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGGCCCCATGCCGCAAGACAGCATCTTGCCCAGTaccacc                                                                                           <  2:417152/126‑1 (MQ=255)                                                ggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGtgatg                                                                               >  1:229384/1‑136 (MQ=255)                                                ggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGtgagg                                                                               >  1:496658/1‑134 (MQ=255)                                                ggcgATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGCGtgatg                                                                               >  1:493531/1‑136 (MQ=255)                                                     tgatgaAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTAc                                                                          <  1:184574/136‑1 (MQ=255)                                                     tgatgaAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTAc                                                                          <  2:580113/136‑1 (MQ=255)                                                     tgatgaAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTAc                                                                          <  1:238645/136‑1 (MQ=255)                                                         gaAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATc                                                                      <  1:39615/136‑1 (MQ=255)                                                               gaACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAg                                                                >  1:234610/1‑136 (MQ=255)                                                               gaACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAg                                                                <  2:455610/136‑1 (MQ=255)                                                               gaACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAg                                                                <  1:484651/136‑1 (MQ=255)                                                                          tCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGAt                                                     <  2:366169/136‑1 (MQ=255)                                                                             gCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTAc                                                                          <  2:375373/112‑1 (MQ=255)                                                                             gCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTAc                                                                          >  1:375373/1‑112 (MQ=255)                                                                                       cATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGa                                        <  1:376870/136‑1 (MQ=255)                                                                                                  gTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATcatcat                             >  1:492530/1‑136 (MQ=255)                                                                                                    aaaGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATAc                           <  2:160166/136‑1 (MQ=255)                                                                                                                   tCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACgaaga            <  2:125015/136‑1 (MQ=255)                                                                                                                     aTTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAat          >  1:278750/1‑136 (MQ=255)                                                                                                                      ttGCCTGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAata         >  2:161452/1‑136 (MQ=255)                                                                                                                             cTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAATATGGCGCt  <  2:21274/136‑1 (MQ=255)                                                                                                                             cTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAATATGGCGCt  <  1:644723/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                               ACCGGGTATAGGAAGTCGCTGGTCACCTGCAGGGTCATGCCCAGCGCGGCGATGATGAAAGAGAACTCGCCAATCTGCGAAAGCCCCATGCCCACGCGTAAAGACGTACGCCCGTCATTGCCGGCTATAAACGCCCCCATGCCGCAAGACAGCATCTTGCCCAGTACCACCGCCAGGGTGATGACTACGATCGGCCAGGCATAGTCGATGAGCACCTCAGGATCAATCATCATACCGATGGCGACGAAGAATATGGCGCT  >  NC_002947/6110382‑6110641

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