Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,213,352	A→C	9.5%	L125V (TTG→GTG)	PP_3695 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,213,352	0	A	C	9.5%	129.1 / 2.5	49	L125V (TTG→GTG)	PP_3695	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base C (3/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGTATCTACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGCTCAC  >  NC_002947/4213210‑4213484                                                                                                                                                |                                                                                                                                     ttGGTATCTACGCCGTACCCGGACAGCAACCGCAGTGGCATTCCTTGGGCGAATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAAccgcc                                                                                                                                 >  1:153309/1‑148 (MQ=255) ttGGTATCTACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGTTGATCATCACccgccc                                                                                                                                >  7:75490/1‑148 (MQ=255)   ggTATCTACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGTTGATCATCACccgccc                                                                                                                                <  8:75490/147‑2 (MQ=255)         tACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTg                                                                                                                        >  6:324726/1‑149 (MQ=255)         tACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTg                                                                                                                        <  5:324726/149‑1 (MQ=255)            gCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAg                                                                                                                     <  8:85222/149‑1 (MQ=255)                   ccGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGgtg                                                                                                              >  7:257757/1‑149 (MQ=255)                    cGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGgtgt                                                                                                             <  5:151165/149‑1 (MQ=255)                      gACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGtt                                                                                                            >  7:79797/1‑148 (MQ=255)                          gCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGGTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAACCTGCAGGAAGGTGTTCAAcc                                                                                                       <  6:209740/149‑1 (MQ=255)                            aaCCGCAGTTGCATTACTTGGGCGCATCCGGGAGATAATCGCTGCGCTAATCCCGCGA‑CCTGATTGACCTCTTCTGAACGCTTATTCTTCCATAGTGACGGGCGGGTGATTATCAACCGCCGATAAACTGCAGGAAGGTGTTCATCCAg                                                                                                     <  2:153309/149‑1 (MQ=255)                            aaCCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAg                                                                                                     >  5:251292/1‑149 (MQ=255)                              ccGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGgt                                                                                                   <  6:278902/149‑1 (MQ=255)                                   gTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTACATAGGGACGGGCGGGTGATGATCAACCGCAGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTc                                                                                               >  6:77359/1‑148 (MQ=255)                                      gCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAg                                                                                            <  8:79797/148‑1 (MQ=255)                                      gCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGc                                                                                           >  2:31883/1‑149 (MQ=255)                                          tCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACa                                                                                       <  5:232152/149‑1 (MQ=255)                                             ttGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATg                                                                                     >  8:102549/1‑148 (MQ=255)                                                 gcgcATCCGGGAGATAATCGCTGCGCTAAACCCGCGATCCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTAtt                                                                                 >  1:56130/1‑149 (MQ=255)                                                  cgcATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTgg                                                                               <  7:81476/149‑1 (MQ=255)                                                     aTCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCa                                                                            <  3:290161/149‑1 (MQ=255)                                                      tCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCAc                                                                           <  6:251292/149‑1 (MQ=255)                                                         gggAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTa                                                                        >  8:152635/1‑149 (MQ=255)                                                             gaTAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTAtctc                                                                    <  8:257757/149‑1 (MQ=255)                                                                aaTCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTAATGGCCACTTATctcct                                                                  >  3:190281/1‑148 (MQ=255)                                                                    gCTGAGCTAACCCCGCGA‑CCAGAACGACCTCATCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCAtt                                                              <  8:250479/148‑1 (MQ=255)                                                                     aggcacTAAACCCTCTA‑CCTGAATGACCGCTTCAGACCGCGTATTCTTACATAGGGACGGGCGGGTTATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAg                                                            <  4:218554/144‑1 (MQ=255)                                                                             aaCCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCt                                                    <  7:152635/149‑1 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCTAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCg                                                      >  3:58223/1‑144 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCTAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCg                                                      <  4:58223/144‑1 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTg                                                                                >  6:102391/1‑118 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTg                                                                                <  5:102391/118‑1 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGt                                                  >  8:201142/1‑148 (MQ=255)                                                                                ccGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTa                                                 <  1:31883/149‑1 (MQ=255)                                                                                 cgcgA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTaa                                                >  2:70015/1‑149 (MQ=255)                                                                                           aaTGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTc                                         <  3:112213/147‑1 (MQ=255)                                                                                                    ttCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACt                                <  5:234639/147‑1 (MQ=255)                                                                                                    ttCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACAGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTgc                              >  5:181738/1‑149 (MQ=255)                                                                                                             cTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACAcc                     <  7:102549/149‑1 (MQ=255)                                                                                                              ttATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCa                    >  5:81739/1‑149 (MQ=255)                                                                                                               tATTCTTGCATCGGGATGGGCGTGGGATGATCACCCGCCGATAAATTTCAGGTAGGTGTTGAACCATGTGTTCCAGCTACATGTATTGGTCACTTATCTCCTCATTAGACATCGGTTTAATCCACTCGGTTGGAGTGCGCAAACTCCa                    <  4:190281/148‑1 (MQ=255)                                                                                                                 ttcttcCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCATGAAGGTGTTCAACCAGTTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGTACTGCGCAAACACCAtt                  <  5:77359/148‑1 (MQ=255)                                                                                                                    ttcCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGt              >  6:50775/1‑149 (MQ=255)                                                                                                                           tttatcggcggtTGATCATCACCCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATg       >  5:248792/7‑149 (MQ=255)                                                                                                                            ggACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGc      >  6:293414/1‑149 (MQ=255)                                                                                                                            ggACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGc      >  5:323279/1‑149 (MQ=255)                                                                                                                              aCGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGCTc    <  6:181738/149‑1 (MQ=255)                                                                                                                                gggcgggTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGCTCAc  <  5:50775/149‑1 (MQ=255)                                                                                                                                cggcggtTGATCATCACCCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTACAGCGACATGTATTGGGAACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGCTCAc  >  7:70071/2‑149 (MQ=255)                                                                                                                                                |                                                                                                                                     TTGGTATCTACGCCGTACCCGGACAGCAACCGCAGTGGCAGTCCTTGGGCGCATCCGGGAGATAATCGCTGCGCTAAACCCGCGA‑CCTGAATGACCTCTTCAGAACGCTTATTCTTCCATAGGGACGGGCGGGTGATGATCAACCGCCGATAAACTGCAGGAAGGTGTTCAACCAGGTGTTCCAGCGACATGTATTGGGCACTTATCTCCTCATTAGACATCGCTGTAATCCACTCGGTGGGACTGCGCAAACACCATTGAGTGGGTATGCTCAC  >  NC_002947/4213210‑4213484

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