Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,019,896	G→T	A72A (GCC→GCA)	PP_4423 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,019,896	0	G	T	97.5%	129.3 / ‑3.8	40	A72A (GCC→GCA)	PP_4423	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (19/20);  total (20/20)
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.21e-01

ACTGTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGG  >  NC_002947/5019750‑5020039                                                                                                                                                   |                                                                                                                                                aCTGTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGc                                                                                                                                               <  1:139490/149‑1 (MQ=255) aCTGTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGc                                                                                                                                               <  5:190222/149‑1 (MQ=255)     tCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCa                                                                                                                                             <  2:361921/147‑1 (MQ=255)       ccGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCtt                                                                                                                                         >  1:187961/1‑149 (MQ=255)       ccGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCtt                                                                                                                                         <  7:238180/149‑1 (MQ=255)          tCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGt                                                                                                                                      <  8:78392/149‑1 (MQ=255)              ccTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAATCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGt                                                                                                                                  <  2:300238/149‑1 (MQ=255)              ccTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCATGGTGAGt                                                                                                                                  <  5:337710/149‑1 (MQ=255)                              aCGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTc                                                                                                                    >  2:15161/1‑147 (MQ=255)                              aCGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTc                                                                                                                    <  1:15161/147‑1 (MQ=255)                                          cccGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCa                                                                                                                   <  2:30903/136‑1 (MQ=255)                                          cccGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTc                                                                                                              >  1:2138/1‑141 (MQ=255)                                          cccGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTc                                                                                                              <  2:2138/141‑1 (MQ=255)                                          cccGCGATCAATCGGAGAGGTTCTGGAGGCCTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCa                                                                                                                   >  1:30903/1‑136 (MQ=255)                                                    aTCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAGGGGCAAAATAATGATTAGCCCTCTGACTGTTTAGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccgc                                                                                             >  1:86452/1‑148 (MQ=255)                                                      cGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccg                                                                                              >  8:238261/1‑145 (MQ=255)                                                      cGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccg                                                                                              >  8:232278/1‑145 (MQ=255)                                                      cGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccg                                                                                              <  7:238261/145‑1 (MQ=255)                                                      cGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccg                                                                                              <  7:232278/145‑1 (MQ=255)                                                        gagaGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGt                                                                                         >  8:21238/1‑148 (MQ=255)                                                                      gCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGACAGCTTGGTGAGTGCAATCGGTCCTTCATAATCGTCTCCATGATTGCCGCAAGTGAGGAGCACAGTctg                                                                          >  5:266049/1‑147 (MQ=255)                                                                       cGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCggg                                                                         <  3:85077/149‑1 (MQ=255)                                                                        gTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGc                                                                        <  2:187961/149‑1 (MQ=255)                                                                           tttAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGa                                                                     >  6:319577/1‑149 (MQ=255)                                                                           tttAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGAGTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCGTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGa                                                                     >  2:102284/1‑149 (MQ=255)                                                                              aaaaGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATcc                                                                  >  7:177719/1‑149 (MQ=255)                                                                                 aGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCAAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGt                                                               >  4:121370/1‑149 (MQ=255)                                                                                       aTAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCa                                                         <  5:319577/149‑1 (MQ=255)                                                                                             cATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACg                                                    <  7:21238/148‑1 (MQ=255)                                                                                                 aaGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGa                                                <  3:121370/148‑1 (MQ=255)                                                                                                              agATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATc                                  <  5:284924/148‑1 (MQ=255)                                                                                                                      ccTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCa                          >  7:106937/1‑149 (MQ=255)                                                                                                                                 ttCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGtca                >  4:255161/1‑148 (MQ=255)                                                                                                                                         ccAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATATGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAg       >  7:329753/1‑149 (MQ=255)                                                                                                                                          cAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAg       >  6:393937/1‑148 (MQ=255)                                                                                                                                           aGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt          >  1:221294/1‑144 (MQ=255)                                                                                                                                           aGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGt          <  2:221294/144‑1 (MQ=255)                                                                                                                                            gTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTaa    <  2:86452/149‑1 (MQ=255)                                                                                                                                            gTTGATTTGCCAGCTTGGTGAGTCCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCCTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGt      >  6:182716/1‑147 (MQ=255)                                                                                                                                                aTTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGCTCCCGTTTTTCATGACGGTGATAGGAATGAGGATCGCGCCCCCCGCAGAGTCGTCAGGTGAGTAggg  >  7:267727/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                ACTGTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAGTAAGG  >  NC_002947/5019750‑5020039

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