Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,399,798	T→A	8.5%	V204E (GTA→GAA)	PP_4743 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,399,798	0	T	A	8.5%	135.1 / 3.5	47	V204E (GTA→GAA)	PP_4743	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/26);  new base A (2/2);  total (19/28)
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.96e-01

ATTACTAATTTTAGTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTCGAGCTCCAAGACGACAAGATTGTAGTCAATAAAGTCAAT  >  NC_002947/5399653‑5399940                                                                                                                                                   |                                                                                                                                               aTTACTAATTTTAGTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTAc                                                                                                                                              >  2:246769/1‑148 (MQ=255)     cTAATTTTACTGATGGTCTGATCCAGTCGTGCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTaga                                                                                                                                          <  8:71425/148‑1 (MQ=255)        aTTTTAGTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGCAAGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGAt                                                                                                                                       <  1:280935/149‑1 (MQ=255)            tAGTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACt                                                                                                                                   <  2:348297/148‑1 (MQ=255)              gTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGa                                                                                                                                 <  2:26075/148‑1 (MQ=255)                    gTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCg                                                                                                                          >  4:297320/1‑149 (MQ=255)                             aGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTa                                                                                                                                               <  8:197346/119‑1 (MQ=255)                              gTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAgaga                                                                                                                                        <  2:4646/125‑1 (MQ=255)                              gTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAgaga                                                                                                                                        >  1:4646/1‑125 (MQ=255)                                cgttCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAAc                                                                                                              <  4:241072/146‑1 (MQ=255)                                 ggTCTACATCTATTACATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCACGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAAc                                                                                                              >  5:131051/1‑148 (MQ=255)                                   tCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACca                                                                                                            >  3:176019/1‑148 (MQ=255)                                      aCATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACAt                                                                                                         >  5:134985/1‑148 (MQ=255)                                         tCTATTCCATCGTCGAGTTATTCAAAGAGCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACt                                                                                                      >  1:201652/1‑148 (MQ=255)                                             ttCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTaa                                                                                                    <  6:134985/146‑1 (MQ=255)                                                cATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAAc                                                                                                              >  8:69596/1‑133 (MQ=255)                                                cATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAAc                                                                                                              <  7:69596/133‑1 (MQ=255)                                                cATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAg                                                                                                >  1:260347/1‑147 (MQ=255)                                                cATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTa                                                                                              <  2:85065/149‑1 (MQ=255)                                                 aTCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAg                                                                                             <  1:246769/149‑1 (MQ=255)                                                    gtcgAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGcc                                                                                           <  2:260347/148‑1 (MQ=255)                                                                aaaGCCCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGAACTAGAGATCACTGATGCAGCGGCTaaaa                                                                                                                   >  7:328744/1‑112 (MQ=255)                                                                aaaGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGa                                                                               >  7:313606/1‑148 (MQ=255)                                                                aaaGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGAc                                                                              >  7:274747/1‑149 (MQ=255)                                                                aaaGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGAc                                                                              >  6:156636/1‑149 (MQ=255)                                                                aaaGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGAACTAGAGATCACTGATGCAGCGGCTaaaa                                                                                                                   <  8:328744/112‑1 (MQ=255)                                                                  aGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCa                                                                            <  6:131051/149‑1 (MQ=255)                                                                      cGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCAc                                                                        <  8:274747/149‑1 (MQ=255)                                                                         aaTTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGa                                                                      <  5:156636/148‑1 (MQ=255)                                                                         aaTTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGa                                                                      <  8:313606/148‑1 (MQ=255)                                                                               cTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAAt                                                                                                   <  3:335344/113‑1 (MQ=255)                                                                               cTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAAt                                                                                                   >  4:335344/1‑113 (MQ=255)                                                                                    tttGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGaa                                                           >  3:134/1‑148 (MQ=255)                                                                                    tttGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGAACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGaa                                                           >  4:62985/1‑148 (MQ=255)                                                                                           gACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCtca                                                    <  4:134/148‑1 (MQ=255)                                                                                           gACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGAACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCtca                                                    <  3:62985/148‑1 (MQ=255)                                                                                             cGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCtcatc                                                  <  3:297320/148‑1 (MQ=255)                                                                                               aGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCAc                                                                        >  5:282936/1‑124 (MQ=255)                                                                                               aGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCAc                                                                        <  6:282936/124‑1 (MQ=255)                                                                                                     gATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGa                                                                      <  1:153553/121‑1 (MQ=255)                                                                                                     gATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGa                                                                      >  2:153553/1‑121 (MQ=255)                                                                                                      attcattcCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTc                                         <  4:176019/149‑1 (MQ=255)                                                                                                                   tCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTc                                         >  8:243791/1‑136 (MQ=255)                                                                                                                   tCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTc                                         <  7:243791/136‑1 (MQ=255)                                                                                                                         aTGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTCGAGCTCCAAGACGACAAGa                      <  2:201652/149‑1 (MQ=255)                                                                                                                            gCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTCGAGCTCCAAGACGACAAGATTg                   <  1:315988/149‑1 (MQ=255)                                                                                                                                             aGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTCGAGCTCCAAGACGACAAGATTGTagtcaataaagtcaat  <  1:162401/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               ATTACTAATTTTAGTGATGGTCAGTTCCAGTCGGTCTACATCTATTCCATCGTCGAGTTATTCAAAGACCGAAATTGTCTAACTTTGCTAGACGAGC‑CGGATTCATTCCTACATCCCGAATGGCAGTTTGAATTCCCGAAGCAAGTACTAGAGATCACTGATGCAGCGGCTAAAAGCAACCACATACTAATGAGTAGCCATAGTGCTTCGACCATCACGACTGCAAACGAAAACCTCATCAACTTATTCGAGCTCCAAGACGACAAGATTGTAGTCAATAAAGTCAAT  >  NC_002947/5399653‑5399940

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