Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,405,230	A→C	10.8%	H76P (CAC→CCC)	PP_2108 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,405,230	0	A	C	10.8%	151.5 / 10.6	56	H76P (CAC→CCC)	PP_2108	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (2/4);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CGCCCCTCGAAGCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAC  >  NC_002947/2405083‑2405374                                                                                                                                                    |                                                                                                                                                 cGCCCCTCGAAGCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCAc                                                                                                                                                 >  1:183170/1‑149 (MQ=255)         gAAGCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAg                                                                                                                                          <  7:132154/148‑1 (MQ=255)             ccAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCTg                                                                                                                                     <  4:39837/149‑1 (MQ=255)              cAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCTgg                                                                                                                                    <  5:189260/149‑1 (MQ=255)               aGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTggc                                                                                                                                   >  4:258906/1‑149 (MQ=255)               aGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTggc                                                                                                                                   <  1:317435/149‑1 (MQ=255)               aGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGGCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTggc                                                                                                                                   <  2:190959/149‑1 (MQ=255)                 ggCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTggcgg                                                                                                                                 >  8:262491/1‑149 (MQ=255)                   ccAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGgcggcg                                                                                                                               >  8:220479/1‑149 (MQ=255)                     aaCCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGa                                                                                                                             >  2:64314/1‑149 (MQ=255)                      aCCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGAc                                                                                                                            <  5:3504/149‑1 (MQ=255)                       ccTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACt                                                                                                                           <  2:183170/149‑1 (MQ=255)                        cTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACtt                                                                                                                          >  6:303472/1‑149 (MQ=255)                         tGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCt                                                                                                                                      >  6:355122/1‑136 (MQ=255)                         tGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCt                                                                                                                                      <  5:355122/136‑1 (MQ=255)                            aaTTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      <  7:220479/149‑1 (MQ=255)                            aaTTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      <  4:320002/149‑1 (MQ=255)                             aTTACTGTGTGGTGCGCTGGAACAGGTAGGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACgg                                                                                                                     >  1:208314/1‑149 (MQ=255)                                    tgtgGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      >  3:227330/1‑141 (MQ=255)                                    tgtgGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      <  4:227330/141‑1 (MQ=255)                                      tggtgCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAg                                                                                                             <  4:184057/148‑1 (MQ=255)                                           cgcctcaaCAGGTACGTGCGTTGCTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGGCAATCGGTTGATCTGTCGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGc                                                                                                       <  8:140087/143‑1 (MQ=255)                                            gcTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGTAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCTggctatc                                                                                                                               <  1:112594/124‑5 (MQ=255)                                            gcTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCCCTTGCACGTGCTggctatc                                                                                                                               >  2:112594/1‑120 (MQ=255)                                                gAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCcgcg                                                                                                  >  4:187730/1‑149 (MQ=255)                                                   cAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCga                                                                                               <  3:187730/149‑1 (MQ=255)                                                       tACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGAcc                                                                                           <  8:98220/149‑1 (MQ=255)                                                          gTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTc                                                                                        >  3:342185/1‑149 (MQ=255)                                                              gttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGc                                                                                    >  8:223749/1‑149 (MQ=255)                                                               ttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATc                                                                                                                >  4:153992/1‑120 (MQ=255)                                                               ttgttTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATc                                                                                                                <  3:153992/120‑1 (MQ=255)                                                                     gACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAAt                                                                                <  1:65245/146‑1 (MQ=255)                                                                     gACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAAt                                                                                >  2:65245/1‑146 (MQ=255)                                                                          gcCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCgaga                                                                                             <  8:272516/128‑1 (MQ=255)                                                                          gcCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCgaga                                                                                             >  7:272516/1‑128 (MQ=255)                                                                             catcatcatCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTg                                                                     <  8:112907/149‑1 (MQ=255)                                                                                    atcaGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGcc                                                              <  2:208314/149‑1 (MQ=255)                                                                                     tcaGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCACCAAGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCAGGCGAATCAGCAAAAGTGCGGTGccc                                                             >  6:379677/1‑149 (MQ=255)                                                                                        gcagcGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTc                                                          >  5:76136/1‑149 (MQ=255)                                                                                            cGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGATGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATc                                                      <  7:262491/149‑1 (MQ=255)                                                                                               gTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGc                                                   <  1:64314/149‑1 (MQ=255)                                                                                                      aGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      >  8:265054/1‑75 (MQ=255)                                                                                                      aGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACg                                                                                                                      <  7:265054/75‑1 (MQ=255)                                                                                                          aaaGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGaa                                        <  5:303472/149‑1 (MQ=255)                                                                                                          aaaGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGaa                                        >  6:61652/1‑149 (MQ=255)                                                                                                               gTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGa                                   <  3:258906/149‑1 (MQ=255)                                                                                                                  aaTCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACa                                 <  4:337717/148‑1 (MQ=255)                                                                                                                   aTCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGa                                   >  3:2109/1‑145 (MQ=255)                                                                                                                    tCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGAt                              >  2:124238/1‑149 (MQ=255)                                                                                                                       gTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATggg                           >  6:54967/1‑149 (MQ=255)                                                                                                                        ttGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGt                          <  7:223749/149‑1 (MQ=255)                                                                                                                             cGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCAt                      <  4:2109/148‑1 (MQ=255)                                                                                                                                  gcgcCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGCACCCATCTGTc                 <  4:342185/148‑1 (MQ=255)                                                                                                                                    gcCTCCGATAGCCAGCACGTGCAAGGGCTGNCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGt              >  1:34153/1‑149 (MQ=255)                                                                                                                                         cGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGAcgc         >  2:196551/1‑149 (MQ=255)                                                                                                                                                cAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAc  >  4:347684/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGCCCCTCGAAGCCAGGGCCAACCTGCAATTACTGTGTGGTGCGCTGGAACAGGTACGTGCGTTGTTTGACGCGCCCATCATCATCAGCAGCGGGTATCGCAGCGAAAGGGTCAATCGGTTGATCGGTGGCGCCTCCGATAGCCAGCACGTGCAAGGGCTGGCGGCGGACTTCACGGTGATCGAGGTCAGCCCGCGCGAGACCGTCCGGCGAATCAGCAAAAGTGCGGTGCCCTTCGATCAGCTGATCCTGGAATTCGACAGATGGGTGCATCTGTCGGTGACGCGCGGCAC  >  NC_002947/2405083‑2405374

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