Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,660,183	G→A	7.8%	F51F (TTC→TTT)	PP_2331 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,660,183	0	G	A	7.8%	148.2 / 2.6	51	F51F (TTC→TTT)	PP_2331	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (2/2);  total (27/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 = 1.00e+00

GACCAGGTAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTAAGAGAGC  >  NC_002947/2660037‑2660331                                                                                                                                                   |                                                                                                                                                      gACCAGGTAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGCCCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGaa                                                                                                                                                     >  3:227175/1‑149 (MQ=255)   ccAGGTAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTATTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTTGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACt                                                                                                                                                   <  7:86259/149‑1 (MQ=255)    cAGGTAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑c                                                                                                                                                 >  2:176238/1‑149 (MQ=255)        tAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGAcc                                                                                                                                             <  5:172477/149‑1 (MQ=255)            cTTGCGGTAGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTAGGAGATGAAGCTCTCCTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTCCACCTTGACCGTGTCCATGCTGCTGGCCACGAACA‑CGACCttcg                                                                                                                                         >  4:207977/1‑147 (MQ=255)               gCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGATCTCGTTGAGGATGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTAGGCGt                                                                                                                                      >  6:164186/1‑149 (MQ=255)                   ttGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGAt                                                                                                                                  >  6:209870/1‑149 (MQ=255)                       cGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTc                                                                                                                              <  8:209174/149‑1 (MQ=255)                        gTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTc                                                                                                                              >  4:287570/1‑148 (MQ=255)                         tCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCca                                                                                                                            <  1:176238/149‑1 (MQ=255)                            gCACGGTTGATGCTGACCCGGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACTTCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAg                                                                                                                         <  3:207977/149‑1 (MQ=255)                                    gATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGtcgtcg                                                                                                                  <  1:149460/148‑1 (MQ=255)                                       gCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGttt                                                                                                              >  6:42874/1‑149 (MQ=255)                                       gCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTACCCAGGGCGTCGGtgt                                                                                                              >  6:251398/1‑147 (MQ=255)                                        cTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTc                                                                                                             <  3:287570/149‑1 (MQ=255)                                           aCCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATc                                                                                                          >  1:186737/1‑149 (MQ=255)                                             cccGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGTTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGCTTCATCCt                                                                                                        >  4:266342/1‑149 (MQ=255)                                              ccGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTc                                                                                                       <  5:240876/149‑1 (MQ=255)                                               cGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCg                                                                                                      >  5:319209/1‑149 (MQ=255)                                                               gAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑AGACATTGGCGTCGATTTTCCACAGGTCGTCGGTTTCATCCTCGGTCAACGGGATGACGa                                                                                      >  6:347036/1‑149 (MQ=255)                                                                 aGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGc                                                                                    <  5:42874/149‑1 (MQ=255)                                                                  gCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGACCT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGACCa                                                                                   <  5:251398/149‑1 (MQ=255)                                                                  gCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCa                                                                                   <  4:227175/149‑1 (MQ=255)                                                                      tcGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACg                                                                                 <  2:186737/147‑1 (MQ=255)                                                                         ttGAGGCTGACGTAGTCGCGGGTCAGCGGGGGCACGCACATCTGCACCTTGACCGGGGCCTTGCTGCTGGCCCCGAACT‑CGACCTGGGCGGCGAGGTTCCCCAGGTCGTCGGTTTCCTCCGCGGGCACCGGGCTGCCGAGCCCGCAGag                                                                            >  2:326312/1‑148 (MQ=255)                                                                            aGGGTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACAAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCt                                                                          >  2:311906/1‑148 (MQ=255)                                                                            aGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACAAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCt                                                                          >  2:317128/1‑148 (MQ=255)                                                                                tGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGAt                                                                       >  6:364068/1‑147 (MQ=255)                                                                                tGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGAt                                                                       >  5:305889/1‑147 (MQ=255)                                                                                tGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATACTCGGTCACCGGGATGCCGAGCACGAAGATCTGATa                                                                      >  4:276926/1‑148 (MQ=255)                                                                                     tAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCAt                                                                                                           <  6:272707/106‑1 (MQ=255)                                                                                     tAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGt                                                                <  3:276926/149‑1 (MQ=255)                                                                                              ttCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCAc                                                       <  5:364068/149‑1 (MQ=255)                                                                                              ttCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCAc                                                       <  6:319209/149‑1 (MQ=255)                                                                                                     gTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCaa                                                 >  3:331575/1‑148 (MQ=255)                                                                                                     gTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGACCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGt                                                                                                    >  6:348729/1‑97 (MQ=255)                                                                                                     gTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGACCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGt                                                                                                    <  5:348729/97‑1 (MQ=255)                                                                                                      tGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATa                                               <  5:209870/149‑1 (MQ=255)                                                                                                        gCACGAACATCTGCACCTTGACCGGGTCCTTGCTGTTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGCTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAg                                              <  3:266342/148‑1 (MQ=255)                                                                                                          aCGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCAc                                           >  3:326275/1‑149 (MQ=255)                                                                                                           cGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACAAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCAcc                                          <  1:317128/149‑1 (MQ=255)                                                                                                           cGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACAAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCAcc                                          <  1:311906/149‑1 (MQ=255)                                                                                                                      aCCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCTTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGAt                               >  6:310756/1‑149 (MQ=255)                                                                                                                       ccTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATc                              <  4:331575/149‑1 (MQ=255)                                                                                                                        cTTGACCGGGTCCTTGCTGCTGGCCACGAACTGCGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGAt                               >  1:226838/1‑148 (MQ=255)                                                                                                                        cTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCa                             <  5:347036/149‑1 (MQ=255)                                                                                                                           gACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAgg                           <  6:305889/148‑1 (MQ=255)                                                                                                                              cGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACtt                       >  4:109282/1‑149 (MQ=255)                                                                                                                                          gctgGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAggg           >  4:123303/1‑149 (MQ=255)                                                                                                                                                cACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTaa        >  1:272357/1‑146 (MQ=255)                                                                                                                                                cACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTaa        >  4:137989/1‑146 (MQ=255)                                                                                                                                                   gAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTAAGAGAgc  <  3:109282/149‑1 (MQ=255)                                                                                                                                                   gAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGACACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTAAGAGAgc  <  5:164186/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GACCAGGTAACCTTGCGGTTGCCGTCGGCACGGTTGATGCTGACCCCGTAGTTGTTGGAGATGAAGCTCTCGTTGAGGCTGACGTAGTCGCGGTTCAGCGGTGGCACGAACATCTGCACCTTGACCGGGTCCTTGCTGCTGGCCACGAACT‑CGACCTTGGCGTCGATGTTCCACAGGTCGTCGGTTTCATCCTCGGTCACCGGGATGCCGAGCACGAAGATCTGATAAGCCGTGACCGCCACGCCCAATAGCACCAGCACGGTGATCAGGACTTTCAGATGGAGGGTAAGAGAGC  >  NC_002947/2660037‑2660331

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