Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,004,526	A→T	6.8%	Q184L (CAG→CTG)	PP_0866 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,004,526	0	A	T	6.8%	218.5 / 5.4	73	Q184L (CAG→CTG)	PP_0866	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (3/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

GATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004381‑1004674                                                                                                                                                  |                                                                                                                                                      gATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑ggg                                                                                                                                                   <  4:212439/149‑1 (MQ=255)  aTGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑ggg                                                                                                                                                   >  7:331325/1‑148 (MQ=255)      cTCGCATCTGGAGGTGAATGGCATCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAgtg                                                                                                                                              <  8:349129/149‑1 (MQ=255)       tCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAgtgt                                                                                                                                             >  6:31854/1‑149 (MQ=255)         gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCCGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCT‑GGGCAgtgt                                                                                                                                             <  3:219778/147‑1 (MQ=255)           aTCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATc                                                                                                                                         <  1:197306/149‑1 (MQ=255)                      aaTGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTca                                                                                                                                     <  8:255958/142‑1 (MQ=255)                      aaTGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTca                                                                                                                                     >  7:255958/1‑142 (MQ=255)                      aaTGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GCGCAGTGTGATCGTCACAGGTAc                                                                                                                              >  7:334987/1‑149 (MQ=255)                         ggCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCg                                                                                                                           >  2:272832/1‑149 (MQ=255)                               acGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGa                                                                                                                     <  4:251557/149‑1 (MQ=255)                               acGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGa                                                                                                                     >  1:245953/1‑149 (MQ=255)                                     ttGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAg                                                                                                                                  <  6:374611/130‑1 (MQ=255)                                     ttGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAg                                                                                                                                  >  5:374611/1‑130 (MQ=255)                                       gCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTg                                                                                                                      <  4:369792/140‑1 (MQ=255)                                       gCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTg                                                                                                                      >  3:369792/1‑140 (MQ=255)                                       gCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGcc                                                                                                             >  6:293927/1‑149 (MQ=255)                                       gCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGc                                                                                                              <  5:31854/148‑1 (MQ=255)                                        cAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGcc                                                                                                             <  2:245953/148‑1 (MQ=255)                                         aGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTg                                                                                                           >  4:352090/1‑149 (MQ=255)                                            tCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGt                                                                                                                                       >  7:124293/1‑118 (MQ=255)                                            tCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGt                                                                                                                                       <  8:124293/118‑1 (MQ=255)                                                  tCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACg                                                                                                                   <  1:65442/132‑1 (MQ=255)                                                  tCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACg                                                                                                                   >  2:65442/1‑132 (MQ=255)                                                       gcCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCT‑GGGCAgtgtcga                                                                                                                                          >  6:193481/1‑101 (MQ=255)                                                       gcCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCT‑GGGCAgtgtcga                                                                                                                                          <  5:193481/104‑4 (MQ=255)                                                       gcCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCAc                                                                                             >  1:199186/1‑149 (MQ=255)                                                        ccGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCAcc                                                                                            >  3:341938/1‑149 (MQ=255)                                                                 ggCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCg                                                                                   <  4:236776/149‑1 (MQ=255)                                                                 ggCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCg                                                                                   <  6:290955/149‑1 (MQ=255)                                                                 ggCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCg                                                                                   <  2:203604/149‑1 (MQ=255)                                                                          gcCTGAACGCTGGTGAAGTGATGTTCAGCGGCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAg                                                                          <  8:331325/149‑1 (MQ=255)                                                                            cTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 >  1:290080/1‑140 (MQ=255)                                                                            cTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 <  2:290080/140‑1 (MQ=255)                                                                               aaCGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGtgtg                                                                                      <  2:110285/132‑1 (MQ=255)                                                                               aaCGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTACCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGtgtg                                                                                      >  1:110285/1‑132 (MQ=255)                                                                                      gTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGCTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  5:319068/149‑1 (MQ=255)                                                                                      gTGAAGTGATGATCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTTCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGc                                                                <  1:272832/147‑1 (MQ=255)                                                                                       tGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCggg                                                             >  8:283906/1‑149 (MQ=255)                                                                                            tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAg                                                        >  4:28670/1‑149 (MQ=255)                                                                                              aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAggg                                                      >  7:194920/1‑149 (MQ=255)                                                                                                gTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCTCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGc                                                     >  4:242173/1‑148 (MQ=255)                                                                                                 ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGc                                                     >  2:172170/1‑147 (MQ=255)                                                                                                   cAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                  <  1:172170/148‑1 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACcca                                                                                                     >  4:207501/1‑95 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACcca                                                                                                     <  3:207501/95‑1 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  7:110492/134‑1 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  8:110492/1‑134 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGc                                               <  7:283906/149‑1 (MQ=255)                                                                                                     gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGc                                               >  6:284423/1‑149 (MQ=255)                                                                                                           cccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGAcccc                                         >  1:111677/1‑146 (MQ=255)                                                                                                            ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagc                                                                                                   >  3:268848/1‑90 (MQ=255)                                                                                                            ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagc                                                                                                   <  4:268848/90‑1 (MQ=255)                                                                                                              aTGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGcc                                       >  1:65774/1‑149 (MQ=255)                                                                                                                    gTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCaa                                <  3:28670/149‑1 (MQ=255)                                                                                                                         cGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGAc                           >  1:26100/1‑149 (MQ=255)                                                                                                                         cGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGAc                           <  2:199186/149‑1 (MQ=255)                                                                                                                         cGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGAc                           >  5:121491/1‑149 (MQ=255)                                                                                                                              gTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  7:115984/1‑109 (MQ=255)                                                                                                                              gTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  8:115984/109‑1 (MQ=255)                                                                                                                              gTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCg                      >  4:342624/1‑149 (MQ=255)                                                                                                                                 tgtgGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCggg                    <  6:121491/148‑1 (MQ=255)                                                                                                                                     gTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCg                      <  2:366212/142‑1 (MQ=255)                                                                                                                                     gTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCg                      >  1:366212/1‑142 (MQ=255)                                                                                                                                     gTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAggg               <  8:334987/149‑1 (MQ=255)                                                                                                                                     gTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAg                 <  7:337077/147‑1 (MQ=255)                                                                                                                                     gTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAg                 >  8:337077/1‑147 (MQ=255)                                                                                                                                       tcacacTGCCCT‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGt             >  5:316768/3‑149 (MQ=255)                                                                                                                                        gACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGAGACTCGCCGGCGCAACTGCCGGCCGGGCAGGGGTc            >  8:151423/1‑149 (MQ=255)                                                                                                                                          cacTGCCCT‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTAGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCg           >  2:261821/1‑148 (MQ=255)                                                                                                                                           acTTCCCA‑GTGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGTCGACTCGCCGGCGCAACTGACGCCCGGGCAGGGGTCGCa         >  6:103099/1‑149 (MQ=255)                                                                                                                                               cccA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata       >  4:305511/1‑147 (MQ=255)                                                                                                                                               cccA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata       >  5:9417/1‑147 (MQ=255)                                                                                                                                                  a‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  <  5:293927/149‑1 (MQ=255)                                                                                                                                                  a‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  <  4:341938/149‑1 (MQ=255)                                                                                                                                                  a‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACATCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  >  2:131171/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GATGGCTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCA‑GGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004381‑1004674

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