Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,605,721	Δ1 bp	7.1%	coding (231/2346 nt)	PP_2284 →	tail tubular protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,605,721	0	G	.	7.1%	176.5 / 7.0	56	coding (231/2346 nt)	PP_2284	tail tubular protein B
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base . (2/2);  total (28/28)
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.06e-01

GACCGTCTCTGAGGGACTGATGAAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCTTCATC  >  NC_002947/2605598‑2605865                                                                                                                            |                                                                                                                                                 gACCGTCTCTGAGGGACTGATGAAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCgtg                                                                                                                         >  1:38786/1‑149 (MQ=255)   ccGTCTCTGAGGGACTGATGAAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTa                                                                                                                       <  3:324397/149‑1 (MQ=255)          tGAGGGACTGATGAAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGc                                                                                                                  <  2:38786/147‑1 (MQ=255)                  tgatgaAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGACCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACgg                                                                                                         >  5:96661/1‑148 (MQ=255)                    atgaAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCAc                                                                                                      >  1:29882/1‑149 (MQ=255)                      gaAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACt                                                                                                     >  4:240274/1‑148 (MQ=255)                            ccccccAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCa                                                                                              >  2:163961/1‑149 (MQ=255)                              ccccAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACg                                                                                            >  3:231948/1‑149 (MQ=255)                                 cAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTg                                                                                                                   >  4:171727/1‑123 (MQ=255)                                 cAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTg                                                                                                                   <  3:171727/123‑1 (MQ=255)                                    ccacccacCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCAc                                                                                                                                          >  8:59198/1‑97 (MQ=255)                                    ccacccacCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCAc                                                                                                                                          <  7:59198/97‑1 (MQ=255)                                    ccacccacCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCt                                                                                      <  2:29882/149‑1 (MQ=255)                                        ccaccTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTcc                                                                                  <  3:240274/149‑1 (MQ=255)                                         caccTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTccc                                                                                 >  3:193072/1‑149 (MQ=255)                                                        aCTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCgt                                                                                                                          >  8:106411/1‑93 (MQ=255)                                                        aCTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCgt                                                                                                                          <  7:106411/93‑1 (MQ=255)                                                               ccccATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGc                                                            <  1:163961/148‑1 (MQ=255)                                                                      ggAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGc                                                     <  6:96661/148‑1 (MQ=255)                                                                         gTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGAcc                                                 >  6:43782/1‑149 (MQ=255)                                                                            ggCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCgg                                              >  5:14921/1‑149 (MQ=255)                                                                                    tCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGT‑GCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACt                                     >  6:7556/1‑149 (MQ=255)                                                                                    tCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGT‑GCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACt                                     >  6:2582/1‑149 (MQ=255)                                                                                     ccACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGt                                                                                                               <  2:86441/75‑1 (MQ=255)                                                                                     ccACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGt                                                                                                               >  1:86441/1‑75 (MQ=255)                                                                                       acacCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGAc                                      >  2:94562/1‑146 (MQ=255)                                                                                       acacCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGAc                                      <  1:94562/146‑1 (MQ=255)                                                                                       acacCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACtt                                    >  4:268307/1‑148 (MQ=255)                                                                                        cacCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAAcc                                                                       >  1:333067/1‑112 (MQ=255)                                                                                        cacCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAAcc                                                                       <  2:333067/112‑1 (MQ=255)                                                                                         acCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCAc                                 <  4:231948/149‑1 (MQ=255)                                                                                          ccATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGc                                         <  1:262616/140‑1 (MQ=255)                                                                                          ccATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGc                                         >  2:262616/1‑140 (MQ=255)                                                                                           cATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATg                                                                            <  6:46976/104‑1 (MQ=255)                                                                                           cATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATg                                                                            >  5:46976/1‑104 (MQ=255)                                                                                            aTCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACcgc                              >  5:152685/1‑149 (MQ=255)                                                                                               aaCCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCACGCGACTTCACCGCGAt                           >  2:343856/1‑149 (MQ=255)                                                                                                      aTGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGTGACTTCACCGCGATCACCGTg                    <  5:43782/149‑1 (MQ=255)                                                                                                         ccAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACcg                               <  8:183580/135‑1 (MQ=255)                                                                                                         ccAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACcg                               >  7:183580/1‑135 (MQ=255)                                                                                                            gcgagcgCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCt                                                                   >  5:71682/1‑96 (MQ=255)                                                                                                            gcgagcgCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCt                                                                   <  6:71682/96‑1 (MQ=255)                                                                                                             cgagcgCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGAc              >  3:173633/1‑148 (MQ=255)                                                                                                              gagcgCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTa            <  4:193072/149‑1 (MQ=255)                                                                                                              gagcgCTACCAAGT‑GCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTac           <  5:2582/149‑1 (MQ=255)                                                                                                              gagcgCTACCAAGT‑GCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTac           <  5:7556/149‑1 (MQ=255)                                                                                                               agcgCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTac           <  4:173633/149‑1 (MQ=255)                                                                                                                gcgcTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGa               <  7:106085/144‑1 (MQ=255)                                                                                                                gcgcTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGa               >  8:106085/1‑144 (MQ=255)                                                                                                                gcgcTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTaca          <  3:268307/149‑1 (MQ=255)                                                                                                                   cTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCt       <  6:152685/149‑1 (MQ=255)                                                                                                                    tACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCAcc                                <  4:231572/123‑1 (MQ=255)                                                                                                                    tACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCAcc                                >  3:231572/1‑123 (MQ=255)                                                                                                                       cAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCTTCAt   >  7:101226/1‑149 (MQ=255)                                                                                                                       cAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCTTCAt   <  6:14921/149‑1 (MQ=255)                                                                                                                        aaGTGGCTATCAGCAACGGCGGACTCCGTGTATTTGACGTATACGGCACTGAGCGCACGGTGTCTTTCCCGGATTGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCTTCATc  <  3:235198/149‑1 (MQ=255)                                                                                                                            |                                                                                                                                                 GACCGTCTCTGAGGGACTGATGAAGCGCCCCCCAACCACCCACCTGGCCCGAGTGACTGCTTCCCCATTGGAGTCGGCATTCGTCCACACCATCAACCGCGATGCCAGCGAGCGCTACCAAGTGGCGATCACCAACGGCGGACTCCGTGTATTTGCCGTAGACGGCACTGAGCGCACGGTGTCTTTCCCGGATGGAACCGGCTACCTGGCCGCAAGCGACCCGGCAAGCGACTTCACCGCGATCACCGTGGCGGACTACACCTTCATC  >  NC_002947/2605598‑2605865

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