Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,203,154	G→A	9.3%	intergenic (+377/‑96)	PP_2809 → / → PP_2810	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,203,154	0	G	A	9.3%	118.6 / 2.8	43	intergenic (+377/‑96)	PP_2809/PP_2810	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base A (2/2);  total (20/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 = 9.95e-01

ACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTT  >  NC_002947/3203008‑3203302                                                                                                                                                   |                                                                                                                                                     acacATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                                     <  2:184802/148‑1 (MQ=255) acacATGTTTCCTTACTTCCCAATAATGAGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGt                                                                                                                                                     <  4:100137/148‑1 (MQ=255)  cacaTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTg                                                                                                                                                   <  7:112224/149‑1 (MQ=255)     aTGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGAt                                                                                                                                                 >  1:335013/1‑148 (MQ=255)          tCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                         >  3:350677/1‑147 (MQ=255)           ccTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCC‑‑TGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcat                                                                                                                                         <  4:350677/148‑3 (MQ=255)                    ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                                 >  5:13169/1‑149 (MQ=255)                    ccAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAg                                                                                                                                 >  4:246474/1‑149 (MQ=255)                        taatGTGCAGGCGCGATGCACTCATTAGTCGATAAATAAATTCAACTCTCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGTCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                             <  6:141153/149‑1 (MQ=255)                        taatGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                             >  7:151764/1‑149 (MQ=255)                                aGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCa                                                                                                                                  <  3:263868/136‑1 (MQ=255)                                aGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCa                                                                                                                                  >  4:263868/1‑136 (MQ=255)                                     ccAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGa                                                                                                                <  5:135925/149‑1 (MQ=255)                                              cATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAc                                                                                                       >  8:148273/1‑149 (MQ=255)                                              cATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATAATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAc                                                                                                       <  1:22721/149‑1 (MQ=255)                                                      gATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCAcccttcgg                                                                                                <  6:13169/148‑6 (MQ=255)                                                        taaataaaTTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTc                                                                                             <  7:148273/149‑1 (MQ=255)                                                           ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                                        >  5:163075/1‑100 (MQ=255)                                                           ataaatTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGATCAACTGACCATGTcatg                                                                                                                                        <  6:163075/103‑4 (MQ=255)                                                               atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                             >  1:17257/1‑110 (MQ=255)                                                               atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATg                                                                                                                             <  2:17257/110‑1 (MQ=255)                                                               atTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGagcagc                                                                                      <  3:246474/149‑1 (MQ=255)                                                                      tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAgg                                                                                                >  6:213437/1‑132 (MQ=255)                                                                      tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAgg                                                                                                <  5:213437/132‑1 (MQ=255)                                                                      tCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTcc                                                                                >  3:154770/1‑148 (MQ=255)                                                                        gCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCaa                                                                              <  2:335013/148‑1 (MQ=255)                                                                        gCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCaa                                                                              <  8:151764/148‑1 (MQ=255)                                                                          tGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGcc                                                                           >  7:160291/1‑149 (MQ=255)                                                                                 tctAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGccc                                                                          >  6:177891/1‑143 (MQ=255)                                                                                 tctAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGccc                                                                          <  5:177891/143‑1 (MQ=255)                                                                                          aTCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGa                                                                                           >  5:31064/1‑117 (MQ=255)                                                                                          aTCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGa                                                                                           <  6:31064/117‑1 (MQ=255)                                                                                                cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                                          <  6:154852/112‑1 (MQ=255)                                                                                                cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGag                                                                                          >  5:154852/1‑112 (MQ=255)                                                                                                cTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAAt                                                     <  8:160291/149‑1 (MQ=255)                                                                                                     aCGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGaaa                                                 >  6:65586/1‑148 (MQ=255)                                                                                                      cGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAAc                                                >  8:68528/1‑148 (MQ=255)                                                                                                       gCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAAc                                                >  5:226141/1‑147 (MQ=255)                                                                                                       gCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACaa                                              <  5:65586/149‑1 (MQ=255)                                                                                                                             tgtTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAgtg                        <  6:226141/149‑1 (MQ=255)                                                                                                                                       atgacatgGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATccgc              >  1:228662/1‑149 (MQ=255)                                                                                                                                                  aGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc                      >  6:181450/1‑130 (MQ=255)                                                                                                                                                  aGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTc                      <  5:181450/130‑1 (MQ=255)                                                                                                                                                   gTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTAtttt  >  6:200144/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ACACATGTTTCCTTACTTCCCAATAATGTGCAGGCGCCAGGCACGCATTAGTCGATAAATAAATTCAACTCGCTGATTTCTCTAGCGTAATCCGCCTATGACGCATTCTGGCATGTGGCTTGCGTGTTGATCACATGACATGGTCAGTTGATCATGTCGCTTGATCAGGATGTAGTAGGCGCCGATCTTTCCACCCGAAGGGTCGAGCAGCCGGTCCAAGCCCTTCGTTTGAGGATCCTCCAATGAAACAAGTCACTGCACGGTCGGAAAGTGTCGTATCCGCCCGCCGTATTTT  >  NC_002947/3203008‑3203302

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