Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,039,704	G→T	8.0%	P33H (CCT→CAT) ‡	PP_1812 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,039,704	0	G	T	8.0%	156.0 / 4.1	50	P33H (CCT→CAT) ‡	PP_1812	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base T (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTG  >  NC_002947/2039578‑2039849                                                                                                                                |                                                                                                                                                   ctgaCCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                              <  7:109353/146‑1 (MQ=255)  ttACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAgggccggg                                                                                                                             >  7:79476/1‑149 (MQ=255)      cGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgggcc                                                                                                                           >  4:15804/1‑147 (MQ=255)            ggggCGGGCGGTTTACCGGGCGGGCCTGGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCgggg                                                                                                                   <  3:142200/149‑1 (MQ=255)            ggggAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCgggg                                                                                                                   >  4:118289/1‑149 (MQ=255)              ggAGGGCGTTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGATTGGCCGGGCAGTGGTGGGCAGTGTGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCg                                                                                                                 <  2:225587/149‑1 (MQ=255)               gAGGCCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgggcc                                                                                                                           <  5:313945/138‑1 (MQ=255)               gAGGACGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgggcc                                                                                                                           >  6:313945/1‑138 (MQ=255)                 gggCGGTTGACGGGGCGGGGCTTTTTCGGCGGTCGTTGTTGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑Gaa                                                                                                             <  5:144925/149‑1 (MQ=255)                  ggCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGTGTTCAGGGCCGGGCCGAGCGGGGCG‑Gaaa                                                                                                            >  1:142851/1‑149 (MQ=255)                    cGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGTGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACg                                                                                                          <  2:142851/149‑1 (MQ=255)                      gTTTACCGGGCGGGCCTTGGTCGGCGGGCGTTGTATGAGAATGGCCGGGCAGGGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACATGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACgg                                                                                                         <  1:89593/149‑1 (MQ=255)                        ttACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCa                                                                                                      >  8:305437/1‑149 (MQ=255)                                 gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGttt                                                                                                >  6:176620/1‑146 (MQ=255)                                 gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGttt                                                                                                <  5:176620/146‑1 (MQ=255)                                   gcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTgtggg                                                                                           >  2:84906/1‑149 (MQ=255)                                        gTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGtgggtggg                                                                                       <  1:106612/148‑1 (MQ=255)                                                 gcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAg                                                                             <  3:118289/149‑1 (MQ=255)                                                  cgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGc                                                                            <  7:305437/149‑1 (MQ=255)                                                          gagaATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                  <  6:221683/104‑1 (MQ=255)                                                          gagaATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                  >  5:221683/1‑104 (MQ=255)                                                              aTGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGAcc                                                                <  8:79476/149‑1 (MQ=255)                                                               tGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGAcc                                                                >  3:305322/1‑148 (MQ=255)                                                                     ggCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GCAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCgg                                                         <  1:84906/149‑1 (MQ=255)                                                                        aGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTaa                                                      >  5:99622/1‑149 (MQ=255)                                                                            gtggGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCAt                                                  >  8:181405/1‑149 (MQ=255)                                                                                  aGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCg                                            >  3:266944/1‑149 (MQ=255)                                                                                    tGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCggg                                          <  4:305322/149‑1 (MQ=255)                                                                                    tGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCggg                                          <  4:266944/149‑1 (MQ=255)                                                                                          tCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGtt                                    >  5:259795/1‑149 (MQ=255)                                                                                             gACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                  <  7:186415/148‑1 (MQ=255)                                                                                             gACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                  >  8:186415/1‑148 (MQ=255)                                                                                                  aCGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGTGAAACGGCCACGGTTTGTGGGTGGGTGAAGATGAGCGTTATGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTAtt                             >  1:97067/1‑149 (MQ=255)                                                                                                  aCGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATtg                            >  3:87153/1‑149 (MQ=255)                                                                                                   cGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTgtgg                                                                                            <  8:281323/84‑1 (MQ=255)                                                                                                   cGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTgtgg                                                                                            >  7:281323/1‑84 (MQ=255)                                                                                                               tACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTa               >  5:67482/1‑149 (MQ=255)                                                                                                                    cgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGt                                        <  6:298081/119‑1 (MQ=255)                                                                                                                    cgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGt                                        >  5:298081/2‑120 (MQ=255)                                                                                                                     gCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTAGAGCGTTTTGCTGACCGATTCGGTAAGCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGc          >  5:228006/1‑149 (MQ=255)                                                                                                                      cacaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTg                                                                      <  2:202913/88‑1 (MQ=255)                                                                                                                      cacaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTg                                                                      >  1:202913/1‑88 (MQ=255)                                                                                                                      cacaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGc          >  3:67577/1‑148 (MQ=255)                                                                                                                       acaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTa        <  4:87153/149‑1 (MQ=255)                                                                                                                       acaAGGCCAGGGCCTAGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACCGTTTCTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGc          >  7:245074/1‑147 (MQ=255)                                                                                                                        caAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                  <  3:64715/122‑1 (MQ=255)                                                                                                                        caAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                  >  4:64715/1‑122 (MQ=255)                                                                                                                           ggCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggt    >  5:243390/1‑149 (MQ=255)                                                                                                                             ccAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttg  <  6:259795/149‑1 (MQ=255)                                                                                                                             ccAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttg  <  6:243390/149‑1 (MQ=255)                                                                                                                                |                                                                                                                                                   CTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACA‑GGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCG‑GAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTG  >  NC_002947/2039578‑2039849

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