Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,413,705	G→A	11.6%	D126N (GAC→AAC)	PP_5643 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,413,705	0	G	A	11.6%	87.2 / 2.0	35	D126N (GAC→AAC)	PP_5643	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/13);  new base A (2/2);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TTAACCGTAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCATGCAGTTCAAGGATGAGTACCTAATTCCGCACGCAGCTGG  >  NC_002947/4413565‑4413844                                                                                                                                             |                                                                                                                                            ttaaCCGTAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTAt                                                                                                                                     >  8:59591/1‑149 (MQ=255)   aaCCGTAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTAtt                                                                                                                                    <  1:131677/148‑1 (MQ=255)    aCCGTAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTAc                                                                                                                                  <  8:235229/149‑1 (MQ=255)        tAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTAt                                                                                                                              <  3:182192/149‑1 (MQ=255)         aGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTAt                                                                                                                              <  7:59591/148‑1 (MQ=255)          gCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATg                                                                                                                             >  8:72650/1‑148 (MQ=255)            gccccGACGGCTCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCTAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGt                                                                                                                            <  4:15649/145‑1 (MQ=255)            gACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGtaa                                                                                                                          >  3:15649/1‑148 (MQ=255)               ccGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGTCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAAt                                                                                                                       >  2:221603/1‑149 (MQ=255)                  aCGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAAtt                                                                                                                      >  5:5844/1‑147 (MQ=255)                    ggCGCCCGATTCGTTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCTCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTgaga                                                                                                                  <  5:58201/149‑1 (MQ=255)                     gcgcCCGCTTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGAc                                                                                                                 <  4:193473/149‑1 (MQ=255)                         ccGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGtt                                                                                                             >  4:263404/1‑149 (MQ=255)                                             cccGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGc                                                                                                >  3:15074/1‑142 (MQ=255)                                                          gCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGc                                                                            <  5:260941/149‑1 (MQ=255)                                                                  gtaCACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCt                                                                    >  3:202982/1‑149 (MQ=255)                                                                     cacaGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTaa                                                                 <  7:72650/149‑1 (MQ=255)                                                                          aaaGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTAAACCTTAGCAGTATCAGCGATTATCTTAAATCAg                                                            >  2:158447/1‑149 (MQ=255)                                                                             gCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGa                                                         >  3:131030/1‑149 (MQ=255)                                                                              cccGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGAt                                                        <  6:5844/149‑1 (MQ=255)                                                                               ccGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGCCCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGAt                                                        >  1:174731/1‑148 (MQ=255)                                                                                  ccAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCAACAACTATTACTTa                                                                                                                               >  7:75867/1‑74 (MQ=255)                                                                                  ccAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCAACAACTATTACTTa                                                                                                                               >  7:79640/1‑74 (MQ=255)                                                                                  ccAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCAACAACTATTACTTa                                                                                                                               <  8:75867/74‑1 (MQ=255)                                                                                  ccAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCAACAACTATTACTTa                                                                                                                               <  8:79640/74‑1 (MQ=255)                                                                                        aCGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCtt                                              <  4:202982/149‑1 (MQ=255)                                                                                        aCGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCtt                                              <  4:131030/149‑1 (MQ=255)                                                                                          ggAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCa                                                                                    >  8:105857/1‑109 (MQ=255)                                                                                             aGCTAATTTTCCGCAAACGTTTTAAATTTCATATACCCTAATATAGGCTAAATTTATTACTTCTGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCa                                          <  1:221603/148‑1 (MQ=255)                                                                                               cGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGCGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGAATg                                        >  6:242998/1‑148 (MQ=255)                                                                                                          cacaCGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCATGCAGTTCAAGGAt                            >  8:68301/1‑149 (MQ=255)                                                                                                                         gCATAGACCCGACTATAGGCGACAACTATTACTTATGt                                                                                                                            >  2:218694/1‑38 (MQ=255)                                                                                                                         gCATAGACCCGACTATAGGCGACAACTATTACTTATGt                                                                                                                            <  1:218694/38‑1 (MQ=255)                                                                                                                               aCCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCATGCAGTTCAAGGATGAGTACCTAATTCcgcacg         >  4:131167/1‑147 (MQ=255)                                                                                                                                    aCTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCATGCAGTTCAAGGATGAGTACCTAATTCCGCACGCAGCTgg  >  7:184239/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                            TTAACCGTAGCGACCCGACGGCGCCCGATTCGGTATCCTCAGACCCCGACAACAAAAGCCGCGTAGTACACAGAAAGCCCGCCAACCACGGAAGCGAATATTCCGCACACGTAGTAATAGGCATAGACCCGACTATAGGCGACAACTATTACTTATGTATAATTGAGACTGTTTTCGGCTCCGGGCTACACCTTAGCAGTATCAGCGATTATCTTAAATCAGTGATCAAGTTCCTTCGCATGCAGTTCAAGGATGAGTACCTAATTCCGCACGCAGCTGG  >  NC_002947/4413565‑4413844

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