Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,538,573	A→T	10.1%	intergenic (+63/+58)	PP_4871 → / ← PP_4872	lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,538,573	0	A	T	10.1%	119.5 / 2.4	46	intergenic (+63/+58)	PP_4871/PP_4872	lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/22);  new base T (2/3);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTG  >  NC_002947/5538428‑5538714                                                                                                                                                  |                                                                                                                                              agcagcCGGACGAGCTGCTGGAAGGGATGGAGAGCTTTGTGGCGCCGGTCGGGCCTATGTGGGGCGACAAGCAGCCTGAGTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACAGAAGGCCGCTCCCACAAAGGCAATg                                                                                                                                            <  3:147623/149‑1 (MQ=255)  gcagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgg                                                                                                                                           <  7:180222/149‑1 (MQ=255)           cGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGc                                                                                                                                  >  3:21828/1‑149 (MQ=255)              gctgctGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTACGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                                                               <  3:173053/149‑1 (MQ=255)              gctgcgGGAAGGGATGGACAGCTTGGTGGCGCCGGGCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCtt                                                                                                                               <  4:34057/149‑1 (MQ=37)               ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGc                                                                                                                                  <  8:143708/145‑1 (MQ=38)               ctgctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGc                                                                                                                                  >  7:143708/1‑145 (MQ=38)                 gctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtgtg                                                                                                                            >  3:315781/1‑149 (MQ=255)                  ctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTgtgt                                                                                                                           <  6:170911/149‑1 (MQ=255)                     gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                        <  8:99021/149‑1 (MQ=255)                     gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTc                                                                                                                          <  1:163047/147‑1 (MQ=255)                               aCAGTTTTGTGGCGGGGTTCGGCCTTATGTGGGTCTACATGCTGCTTGACTTATGGTTAGAGTTTTGGGACGGTTTCGGGGTCTTTTCGCGACACAAGGCCGCTCCCACAAAGGAAATGGGGTAGGGGCCTTGTGTCGCGATGGGCTGc                                                                                                              <  4:21828/149‑1 (MQ=255)                                           cgcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATggg                                                                                                                                          >  4:20779/1‑109 (MQ=255)                                           cgcCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATggg                                                                                                                                          <  3:20779/109‑1 (MQ=255)                                                     gcCTAAGTGGGTCGACAAGCAGCCTGACTTACGCGTAGACTTTGGGGACCGCTTCGCGGTCCTTTCGCGACACACGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATc                                                                                        <  4:215978/149‑1 (MQ=38)                                                     gcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTcgc                                                                                                                        <  2:73516/117‑1 (MQ=21)                                                     gcCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTcgc                                                                                                                        >  1:73516/1‑117 (MQ=21)                                                           gTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTaa                                                                                  <  3:83544/149‑1 (MQ=255)                                                               gTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACAGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAACGCAATGGGGGACCGGCCTTGTGTCCCGATGGGCTGCAAAGCAGCACCAAAAATCGAACAGTTAAcgag                                                                              >  4:256400/1‑147 (MQ=38)                                                                tCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAcgcg                                                                              >  2:201329/1‑148 (MQ=255)                                                                     aaGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGa                                                                        <  4:315781/149‑1 (MQ=255)                                                                      agcagcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCACCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAAAAGTCGATCAGTTAACGCGGAATGa                                                                        >  1:146929/1‑148 (MQ=255)                                                                           ccTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCt                                                                  >  4:183800/1‑149 (MQ=255)                                                                               actgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCCAAGCAGCCCCAATAGTCGATCAGTTTACGCGGCATGACAGGCTGAcg                                                              >  7:106407/1‑149 (MQ=255)                                                                                   aaGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGc                                                          <  1:201329/149‑1 (MQ=255)                                                                                     gCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGc                                                          <  8:186533/147‑1 (MQ=255)                                                                                     gCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGc                                                          >  7:186533/1‑147 (MQ=255)                                                                                      cttcgaCTTTTGGGACCGCTTCGCGGTCCTTTCGCCACACAAGGCCGCTCCCACAAAGGCAATGGGGGCGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttc                                                       <  5:280541/145‑1 (MQ=255)                                                                                        gagaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGAAGGGCTGCAAAGCAGCCCCAAAAGTCGATCAGTTAAACCGGAATGACAGGCTGACGCGGcttctt                                                     >  3:72176/1‑149 (MQ=255)                                                                                         agtcTTTTGGGTCGTATTCGTGGTCGTTTCGGGAATCATGGCCGCTCCCACAATGGGAATGGGTGAGGGGCGTTGTGTCTCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGcttctt                                                     <  3:256400/145‑1 (MQ=255)                                                                                         agaCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACCGGCTGACGCGGCTTCttg                                                    >  7:138545/1‑149 (MQ=255)                                                                                               ttGCGACCGCTTCGCGGTCCTTGCGCGACACAAGCCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGcc                                               <  8:127535/148‑1 (MQ=255)                                                                                                  ggACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGCATGACAGGCTGACGCGGCTGCTTGTTGCCCTTg                                           >  5:93567/1‑149 (MQ=255)                                                                                                            gcgGTCCTTTCGCGACACAAGGCCGCTCCCACAACGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgttg                                                 <  6:252269/133‑1 (MQ=255)                                                                                                            gcgGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCttgttg                                                 >  5:252269/1‑133 (MQ=255)                                                                                                                   tttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAAt                                                                          >  3:261686/1‑101 (MQ=255)                                                                                                                   tttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAAt                                                                          <  4:261686/101‑1 (MQ=255)                                                                                                                    ttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAcc                                      >  2:227480/1‑136 (MQ=255)                                                                                                                    ttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCAcc                                      <  1:227480/136‑1 (MQ=255)                                                                                                                        cgACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCATGTTGCCCTTGCCACCCTTGCTCGCCTCTTTgc                     >  3:25336/1‑149 (MQ=255)                                                                                                                            acacGGCCGCTCCCACAATGGAATGGGGGGTGCGGCGTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTc                 <  2:146929/149‑1 (MQ=39)                                                                                                                                   cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg            >  1:326850/1‑147 (MQ=255)                                                                                                                                    gCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAg         >  1:154585/1‑149 (MQ=255)                                                                                                                                           caAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgctg  <  2:326850/149‑1 (MQ=255)                                                                                                                                                gCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGc           >  3:247180/1‑135 (MQ=255)                                                                                                                                                gCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGc           <  4:247180/135‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              AGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTG  >  NC_002947/5538428‑5538714

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