Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,995,504	T→A	7.3%	intergenic (‑756/+163)	PP_1781 ← / ← PP_mr26	acyltransferase/C4 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,995,504	0	T	A	7.3%	154.5 / 2.9	55	intergenic (‑756/+163)	PP_1781/PP_mr26	acyltransferase/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base A (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GTTAAGCCTCAATGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTGAGCTGCT  >  NC_002947/1995371‑1995648                                                                                                                                       |                                                                                                                                                 gTTAAGCCTCAATGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGc                                                                                                                                     >  5:99513/1‑147 (MQ=255)       ccTCAATGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCa                                                                                                                             >  5:160322/1‑149 (MQ=255)          cAATGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACgg                                                                                                                          >  1:5748/1‑149 (MQ=255)           aaTGAGACGTTTTTGGCTAAGAGGGCATTATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTAATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACgg                                                                                                                          >  1:20206/1‑149 (MQ=255)           aaTGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACggg                                                                                                                         >  3:142901/1‑149 (MQ=255)           aaTGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACggg                                                                                                                         >  2:66383/1‑149 (MQ=255)           aaTGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACACCGCTTGATGCACgg                                                                                                                          <  7:273586/148‑1 (MQ=255)             tGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTAATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTAGATgacatcgg                                                                                                                        >  3:296564/1‑141 (MQ=255)                   gTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAgcgg                                                                                                                 >  2:272445/1‑149 (MQ=255)                        tGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCa                                                                                                             >  4:129446/1‑148 (MQ=255)                              agagGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAgcgc                                                                                                       <  6:297161/148‑1 (MQ=255)                              agagGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCAGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGcgcg                                                                                                      >  5:19919/1‑149 (MQ=255)                              agagGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCAGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGcgcg                                                                                                      >  5:24487/1‑149 (MQ=255)                                         tGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATtgtg                                                                                             <  6:24075/148‑1 (MQ=255)                                         tGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTgg                                                                                            <  2:41847/149‑1 (MQ=255)                                          gTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGc                                                                                           >  3:183527/1‑149 (MQ=255)                                           tAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCt                                                                                          >  3:67468/1‑149 (MQ=255)                                             aCTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTgg                                                                                        <  1:66383/149‑1 (MQ=255)                                                 cagcagCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGAc                                                                                      <  1:108562/147‑1 (MQ=255)                                                 cagcagCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGAc                                                                                      <  4:142901/147‑1 (MQ=255)                                                   gcagcTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAat                                                                                  <  5:278513/149‑1 (MQ=255)                                                             cccAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGAGCGTAAATTGTGGCTGGACCATTATCTTCAAGa                                                                        >  7:232739/1‑149 (MQ=255)                                                                  aGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAAt                                                                   <  6:99513/149‑1 (MQ=255)                                                                  aGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAAt                                                                   <  2:5748/149‑1 (MQ=255)                                                                   gCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGAAGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATg                                                                  >  4:250991/1‑149 (MQ=255)                                                                       gggCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGc                                                              <  3:214968/149‑1 (MQ=255)                                                                                 gTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTgg                                                    <  6:120866/149‑1 (MQ=255)                                                                                      cTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCa                                                <  5:266890/148‑1 (MQ=255)                                                                                       tGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAg                                               >  6:45442/1‑148 (MQ=255)                                                                                         ttCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAgg                                              <  3:129446/147‑1 (MQ=255)                                                                                            aTACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTa                                         >  7:64014/1‑149 (MQ=255)                                                                                              aCATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCa                                                                                    <  3:79590/104‑1 (MQ=255)                                                                                              aCATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCa                                                                                    >  4:79590/1‑104 (MQ=255)                                                                                                aTGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGAAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTAttt                                      >  1:125822/1‑148 (MQ=255)                                                                                                   gACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCaa                                  <  6:160322/149‑1 (MQ=255)                                                                                                         atcatcaCAATCGACGCTGGGGCCTTGGGCTGTAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATgg                             <  2:125822/148‑1 (MQ=255)                                                                                                         atcatcaCAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATgg                             <  8:64014/148‑1 (MQ=255)                                                                                                         atcatcaCAATCGACGCTGGGGCCTTGGGCAGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATgg                             <  6:19919/148‑1 (MQ=255)                                                                                                         atcatcaCAATCGACGCTGGGGCCTTGGGCAGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATgg                             <  6:24487/148‑1 (MQ=255)                                                                                                            atcaCAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTa                          >  8:36765/1‑148 (MQ=255)                                                                                                              cacaATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCGGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTacac                       >  5:83419/1‑149 (MQ=255)                                                                                                                 aaTCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGAATAGCGAGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTc                     >  8:138621/1‑148 (MQ=255)                                                                                                                  aTCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCa                   <  1:239957/149‑1 (MQ=255)                                                                                                                    cGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAg                  >  4:131168/1‑148 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGTTGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTAtttt                                     >  1:27402/1‑123 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCtt                                                                             <  4:262269/83‑1 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCtt                                                                             >  3:262269/1‑83 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTAtttt                                     <  2:27402/123‑1 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGtt                           >  7:238074/1‑133 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGtt                           <  8:238074/133‑1 (MQ=255)                                                                                                                          tGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGt           >  4:26627/1‑149 (MQ=255)                                                                                                                            tggCCTTGTGCTGCAGCAACATCGCTTTATGCACGGTGATAGCTGCGCATAGCGCGTAAATTGTGGCTGGACAAATATCTTCAAGACAAATGTTGCTGAGATGTTGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTg         <  7:138621/148‑1 (MQ=255)                                                                                                                               ccTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTGAgc      >  1:66839/1‑149 (MQ=255)                                                                                                                               ccTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTGAgc      <  6:179125/149‑1 (MQ=255)                                                                                                                                    ggCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTGAgctgct  <  4:67468/148‑1 (MQ=255)                                                                                                                                       |                                                                                                                                                 GTTAAGCCTCAATGAGACGTTTTTGGCTAAGAGGGCA‑TATGTAACTGCAGCAGCTGCGGCCCAAAGCAAGGGCAGCAGCGTCTCCTGTTCATACATGGACGTCATCATCACAATCGACGCTGGGGCCTTGGGCTGCAGCAACATCGCTTGATGCACGGGGATAGCGGCGCATAGCGCGTAAATTGTGGCTGGACCAATATCTTCAAGACAAATGGTGCTGAGAGGTGGATCAGGAGCTATTTTCAACATGGTTACACTCCAGCAATGGTTGAGCTGCT  >  NC_002947/1995371‑1995648

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