Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,543,060	A→G	10.2%	A280A (GCA→GCG)	PP_1353 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,543,060	0	A	G	10.2%	98.0 / 3.5	39	A280A (GCA→GCG)	PP_1353	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base G (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

TTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGCGTGACCCCGCCCGCCACACCCGATACGCC  >  NC_002947/1542929‑1543193                                                                                                                                    |                                                                                                                                      ttGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCATGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGt                                                                                                                      >  3:77728/1‑149 (MQ=255)         cAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCg                                                                                                              <  4:309702/149‑1 (MQ=255)             aCCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGgc                                                                                                          >  5:199735/1‑149 (MQ=255)                tCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGgcgc                                                                                                        >  6:151255/1‑148 (MQ=255)                        cTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGcc                                                                                               >  6:126581/1‑149 (MQ=255)                         tGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGccc                                                                                              >  1:326421/1‑149 (MQ=255)                                 gATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCtttt                                                                                      >  1:107753/1‑149 (MQ=255)                                     ttCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCt                                                                                         >  4:248537/1‑142 (MQ=255)                                     ttCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCt                                                                                         <  3:248537/142‑1 (MQ=255)                                       cATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGAc                                                                                >  4:67621/1‑149 (MQ=255)                                       cATGATCAACGAGTATCCGCGTGACCGGTTGAAGCACCAAGGGATCGACATTAAGTTTCCTCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAAAGCAGTAATGCCGCGTTCAAGCCTGGGCGCAAAGCGCCCCCAGCTTTTTCCCAc                                                                                >  6:254685/1‑149 (MQ=255)                                            tCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCggg                                                                                                            >  1:169432/1‑116 (MQ=255)                                            tCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCggg                                                                                                            <  2:169432/116‑1 (MQ=255)                                             cAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCgg                                                                                                             >  3:126207/1‑114 (MQ=255)                                             cAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCgg                                                                                                             <  4:126207/114‑1 (MQ=255)                                                  aGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGTTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGTTTacac                                                                     <  2:270090/149‑1 (MQ=255)                                                  aGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTacac                                                                     <  2:326421/149‑1 (MQ=255)                                                  aGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTacac                                                                     <  2:107753/149‑1 (MQ=255)                                                      tGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTg                                                                  <  8:48397/148‑1 (MQ=255)                                                         gcatGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTca                                                              <  4:77728/146‑1 (MQ=255)                                                           gTGACCGGTTGAAGGCCGAAGTGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTcaca                                                            <  6:199735/149‑1 (MQ=255)                                                           gTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCGACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGAcc                                                                               <  4:306648/130‑1 (MQ=255)                                                           gTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCGACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGAcc                                                                               >  3:306648/1‑130 (MQ=255)                                                                   ttGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTc                                                    <  4:244683/149‑1 (MQ=255)                                                                   ttGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTc                                                    <  5:126581/149‑1 (MQ=255)                                                                   ttGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACACTGATTACACTTGCTCACAGCAAttt                                                     >  3:104087/1‑148 (MQ=255)                                                                      aaGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCa                                                  <  5:151255/148‑1 (MQ=255)                                                                      aaGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCa                                                  <  6:160761/148‑1 (MQ=255)                                                                       aGGCCGAATGGATCGACATTCCGTTTCCGCATCGGGTTATTCGCGTTGTGCATGATACTGCAACGCATTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGc                                                <  4:104087/149‑1 (MQ=255)                                                                            gAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGTCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCAttt                                           <  7:217581/149‑1 (MQ=255)                                                                                  aTCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCGACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCt                                                                >  8:271562/1‑122 (MQ=255)                                                                                  aTCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCGACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCt                                                                <  7:271562/122‑1 (MQ=255)                                                                                  aTCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGc                                      <  1:269390/148‑1 (MQ=255)                                                                                  aTCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGc                                      >  2:269390/1‑148 (MQ=255)                                                                                         ttCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGCg                              >  3:264680/1‑149 (MQ=255)                                                                                            cGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGc                               >  1:291229/1‑145 (MQ=255)                                                                                            cGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGc                               <  2:291229/145‑1 (MQ=255)                                                                                                        gggTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGCGTGACccccccccccc               >  5:204868/1‑148 (MQ=255)                                                                                                                      tctcAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCTGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATTGCGTGACCCCGCCCGCCACACCCGATACGcc  <  5:184933/145‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      TTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTCGCGTGGTGCAGGAGACTGCAACGCAGTAATGCCGCGTTCAAGCCGGGGCGCAAAGCGGCCCCAGCTTTTTCCGACCCTGATTACACTTGCTCACAGCAATTTCCAGCCATTTCTGGCATATAGCGTGACCCCGCCCGCCACACCCGATACGCC  >  NC_002947/1542929‑1543193

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