Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,301,102	A→G	6.7%	T75A (ACG→GCG)	PP_3774 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,301,102	0	A	G	6.7%	175.9 / 2.4	60	T75A (ACG→GCG)	PP_3774	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base G (2/2);  total (31/29)
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.40e-01

TGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATC  >  NC_002947/4300966‑4301248                                                                                                                                         |                                                                                                                                                   tGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAgcgc                                                                                                                                       >  1:991906/1‑150 (MQ=255)  gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAgcgtc                                                                                                                                      <  1:644842/150‑3 (MQ=255)  gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGc                                                                                                                                             <  1:1029110/143‑1 (MQ=255)  gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGc                                                                                                                                             >  2:1029110/1‑143 (MQ=255)  gAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGcc                                                                                                                                      >  2:111868/1‑150 (MQ=255)   aaGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGCTGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGccc                                                                                                                                     >  2:302378/1‑150 (MQ=255)     gCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAg                                                                                                                                   <  2:991906/150‑1 (MQ=255)           gCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCAc                                                                                                                             <  2:1676/150‑1 (MQ=255)            cccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                            <  1:111868/150‑1 (MQ=255)            cccGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTCCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACt                                                                                                                            <  1:1161318/150‑1 (MQ=255)                 tcgctgaGCACCGTCACGATGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAAt                                                                                                                       >  2:258641/7‑150 (MQ=255)                  ttCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATc                                                                                                                      >  2:1637277/1‑150 (MQ=255)                    cccAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGc                                                                                                                    <  1:1637277/150‑1 (MQ=255)                                                 cTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTc                                                                                       >  2:850691/1‑150 (MQ=255)                                                  tCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCaa                                                                                          >  2:406435/1‑146 (MQ=255)                                                  tCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCaa                                                                                          <  1:406435/146‑1 (MQ=255)                                                  tCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCg                                                                                      <  2:1498266/150‑1 (MQ=255)                                                    gATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgc                                                                                    >  1:476597/1‑150 (MQ=255)                                                     aTCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGgcg                                                                                   >  1:33605/1‑150 (MQ=255)                                                         aGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTg                                                                                                                           <  1:14438/106‑1 (MQ=255)                                                         aGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTg                                                                                                                           >  2:14438/1‑106 (MQ=255)                                                          gTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGaa                                                                              >  2:1543464/1‑150 (MQ=255)                                                                 gCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCAgtct                                                                                                                                <  1:1052268/93‑4 (MQ=39)                                                                 gCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGGCGCTGGCCAGCGTCCAgtct                                                                                                                                >  2:1052268/1‑90 (MQ=39)                                                                        ccaaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGccc                                                                >  1:1162368/1‑150 (MQ=255)                                                                         caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                               >  1:1263025/1‑150 (MQ=255)                                                                         caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                               >  1:450243/1‑150 (MQ=255)                                                                         caaAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCa                                                               >  2:1451250/1‑150 (MQ=255)                                                                            aGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGAt                                                            <  2:1609946/150‑1 (MQ=255)                                                                                      aGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCCGCGCGCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTAcgc                                                  <  2:465848/150‑1 (MQ=255)                                                                                        cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCATCTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCgg                                                <  2:788913/150‑1 (MQ=255)                                                                                        cAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCgg                                                <  1:288517/150‑1 (MQ=255)                                                                                                 ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCNCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCaan                                                                                         <  1:232868/100‑2 (MQ=255)                                                                                                 ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGt                                                                                             >  2:755943/1‑96 (MQ=255)                                                                                                 ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGt                                                                                             <  1:755943/96‑1 (MQ=255)                                                                                                 ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAAc                                                                                         >  2:232868/1‑100 (MQ=255)                                                                                                 ccACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACgg                                       >  1:1500388/1‑150 (MQ=255)                                                                                                  cacaGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGc                                      >  2:1344877/1‑150 (MQ=255)                                                                                                            gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCAc                                                        >  1:147431/1‑122 (MQ=255)                                                                                                            gCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCAc                                                        <  2:147431/122‑1 (MQ=255)                                                                                                                  aaTGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGgtgt                      <  2:33605/150‑1 (MQ=255)                                                                                                                        gggAATCGAAAGACTGGCCGCTGGCCAGCGAACAGGATCACTGGACTCGCCAGACAATGACGCTGCCCAAGGTCACCTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGa                <  2:1162368/149‑1 (MQ=255)                                                                                                                         ggAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGaa               >  1:442639/1‑150 (MQ=255)                                                                                                                           aaaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacga                                          <  2:111575/121‑1 (MQ=255)                                                                                                                           aaaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGgacga                                          >  1:111575/1‑121 (MQ=255)                                                                                                                            aaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACaa            >  1:348012/1‑150 (MQ=255)                                                                                                                            aaCGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACaa            <  2:1500388/150‑1 (MQ=255)                                                                                                                              cGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAgg          <  2:825080/150‑1 (MQ=255)                                                                                                                                 aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCggg                                               <  1:433276/110‑1 (MQ=255)                                                                                                                                 aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCggg                                               >  2:433276/1‑110 (MQ=255)                                                                                                                                 aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCt                   >  1:1575994/1‑138 (MQ=255)                                                                                                                                 aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCt                   <  2:1575994/138‑1 (MQ=255)                                                                                                                                 aaGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCa       >  1:71414/1‑150 (MQ=255)                                                                                                                                   gACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc     <  1:302378/150‑1 (MQ=255)                                                                                                                                     cTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTc                    >  1:1401256/1‑133 (MQ=255)                                                                                                                                     cTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTc                    <  2:1401256/133‑1 (MQ=255)                                                                                                                                     cTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc     >  1:764029/1‑148 (MQ=255)                                                                                                                                     cTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAAc     <  2:764029/148‑1 (MQ=255)                                                                                                                                      tGGGCGCTGGCCAGCGTCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTAGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATc  >  2:101637/1‑150 (MQ=255)                                                                                                                                      tGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATc  <  1:581342/150‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                   TGAAGCGATTGCCCGCTTTCCCAGCATCGTGACGGTGCTCAGCGATGGCTCGATCCAGTTCTCAGCCCCAACCAAAGGGGCTTCGAGCAAAAGCACCCACAGGACACGCTGCGAATGGACGGAAACGCAAGACTGGACGCTGGCCAGCGCCCAGGATCACTGGAATCGCCAGACAATGACGCTGACCAAGGTCAACTCGGCGCAGAAGGTGGTCATCGCCCAGATGCACGTACGCGGGGACGACGGCCCGCCGGTAAAGGTGTTCTGGAACAAGGGCAACATC  >  NC_002947/4300966‑4301248

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