Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,520,084	T→A	16.3%	intergenic (+703/‑177)	PP_5580 → / → PP_5582	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,520,084	0	T	A	16.3%	97.0 / 12.7	43	intergenic (+703/‑177)	PP_5580/PP_5582	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/13);  new base A (5/2);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

TGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAG  >  NC_002947/3519937‑3520232                                                                                                                                                    |                                                                                                                                                     tGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTAt                                                                                                                                                      <  8:141646/148‑1 (MQ=255)   aCGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTTTgc                                                                                                                                                    <  6:58617/148‑1 (MQ=255)            gAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGATTATGCAAGCAATAGt                                                                                                                                          >  3:177498/1‑149 (MQ=255)                                    tCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAtt                                                                                                                  >  3:54181/1‑149 (MQ=255)                                    tCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAtt                                                                                                                  >  3:131131/1‑149 (MQ=255)                                    tCATATTTCAGAATCGAGTCACAGATATATAAAGACGCATGCAGGACGATGGTGAATATTGGAGCTATGAATATATGTGTTAAGCAGCACGGATAGTATCAATTCTTGCTTATGCAAGCAAAAATGCTACTCTTAAAGAAATAATCaat                                                                                                                  >  8:51457/1‑147 (MQ=255)                                      atatTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGa                                                                                                                <  5:94324/149‑1 (MQ=255)                                       tatTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAg                                                                                                               <  4:54181/149‑1 (MQ=255)                                       tatTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAg                                                                                                               <  4:131131/149‑1 (MQ=255)                                         tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCg                                                                                                             >  3:248171/1‑149 (MQ=255)                                                  cGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCg                                                                                                     >  4:271507/1‑148 (MQ=255)                                                        ccAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGAAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGt                                                                                              >  1:310731/1‑149 (MQ=255)                                                        ccAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGAGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAAAATTCATTGAGCGAGATGGCGAAAGCGt                                                                                              >  6:321720/1‑149 (MQ=255)                                                                   aTGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAAGATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGGTTGt                                                                                   >  4:248011/1‑149 (MQ=255)                                                                                gaCGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACaa                                                                      <  3:271507/149‑1 (MQ=255)                                                                                      ggTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGAc                                                                        <  8:5513/141‑1 (MQ=255)                                                                                      ggTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACaacatcg                                                                 >  7:5513/1‑143 (MQ=255)                                                                                        tGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCa                                                               >  1:292633/1‑148 (MQ=255)                                                                                        tGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGGTGGCTTGTATTTTTACGACAAAAATTCa                                                               >  4:148463/1‑148 (MQ=255)                                                                                             atTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTg                                                         >  3:233816/1‑149 (MQ=255)                                                                                                ggTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCa                                                       <  7:51457/148‑1 (MQ=255)                                                                                                    cTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGATGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCtt                                                   >  5:326112/1‑148 (MQ=255)                                                                                                              cggTGTGAAGCCCCCCGGATAGACACACTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTa                                        <  4:177498/148‑1 (MQ=255)                                                                                                                  gtgAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTg                                    >  4:57437/1‑149 (MQ=255)                                                                                                                         agcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAAtt                             <  3:57437/149‑1 (MQ=255)                                                                                                                          gcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       <  8:54982/138‑1 (MQ=255)                                                                                                                          gcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       >  7:54982/1‑138 (MQ=255)                                                                                                                              ggATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGc                                      <  6:211959/135‑1 (MQ=255)                                                                                                                              ggATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGc                                      >  5:211959/1‑135 (MQ=255)                                                                                                                                   gTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGa                    <  5:321720/148‑1 (MQ=255)                                                                                                                                    tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       >  1:47616/1‑128 (MQ=255)                                                                                                                                    tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAg                                       <  2:47616/128‑1 (MQ=255)                                                                                                                                    tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAtt                  >  4:250736/1‑149 (MQ=255)                                                                                                                                    tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAtt                  >  2:108757/1‑149 (MQ=255)                                                                                                                                    tATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAtt                  >  3:22070/1‑149 (MQ=255)                                                                                                                                      ctattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgt                >  6:177506/5‑149 (MQ=255)                                                                                                                                       tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCtt                                                                                     <  4:283448/76‑1 (MQ=39)                                                                                                                                       tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCtt                                                                                     >  3:283448/4‑79 (MQ=39)                                                                                                                                       tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGaa                               <  4:34187/130‑1 (MQ=255)                                                                                                                                       tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGaa                               >  3:34187/4‑133 (MQ=255)                                                                                                                                        attgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGGGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgttt              >  5:228671/3‑149 (MQ=255)                                                                                                                                         ttgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgtttg             >  8:69659/2‑149 (MQ=255)                                                                                                                                           gcttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATtgtttgtt           >  5:84028/1‑149 (MQ=255)                                                                                                                                               gcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTtgtg        >  8:41811/1‑148 (MQ=255)                                                                                                                                                    tGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAg  <  4:248171/149‑1 (MQ=255)                                                                                                                                                    tGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAg  <  4:22070/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGTCTCAG  >  NC_002947/3519937‑3520232

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