Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,438,210	A→T	5.9%	K457M (AAG→ATG) ‡	PP_3932 →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,438,210	0	A	T	5.9%	194.2 / 2.8	68	K457M (AAG→ATG) ‡	PP_3932	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base T (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTT  >  NC_002947/4438064‑4438353                                                                                                                                                    |                                                                                                                                                cgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAgg                                                                                                                                               <  7:73110/149‑1 (MQ=255)     aTACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCg                                                                                                                                           >  1:149101/1‑149 (MQ=255)       aCAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCggg                                                                                                                                         >  7:60960/1‑149 (MQ=255)        cAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGGTTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCgggg                                                                                                                                        <  2:13723/149‑1 (MQ=255)        cAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCTTGGCCCggta                                                                                                                                        <  3:358074/149‑3 (MQ=255)         aGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGc                                                                                                                                       >  6:212710/1‑149 (MQ=255)           gACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCg                                                                                                                                     <  2:149101/149‑1 (MQ=255)           gACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCg                                                                                                                                     >  8:219782/1‑149 (MQ=255)                    gTGACACAGGCCGCCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCACCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGATGTA‑CCGGGCCAAGGCCCCAGGCCGCACCCGCGt                                                                                                                            >  5:16058/1‑149 (MQ=255)                    gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCTTGGCCCgg                                                                                                                                          >  5:129796/1‑135 (MQ=255)                    gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCTTGGCCCgg                                                                                                                                          <  6:129796/135‑1 (MQ=255)                         acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGaa                                                                                                                        >  8:34497/1‑148 (MQ=255)                         acaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAAc                                                                                                                       <  8:210403/149‑1 (MQ=255)                          caGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                      <  2:79713/149‑1 (MQ=255)                           aGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAg                                                                                                                     <  7:219782/149‑1 (MQ=255)                            ggCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAg                                                                                                                     <  1:63006/149‑1 (MQ=255)                              ccGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                  >  7:325524/1‑149 (MQ=255)                              ccGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                  >  2:83944/1‑149 (MQ=255)                              ccGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CAGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                  >  8:210230/1‑149 (MQ=255)                               cGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                  >  5:356512/1‑148 (MQ=255)                               cGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGcc                                                                                                                  >  3:362962/1‑148 (MQ=255)                                gACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAAc                                                                                                                       <  2:354615/142‑1 (MQ=255)                                gACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAAc                                                                                                                       >  1:354615/1‑142 (MQ=255)                                 aCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCTTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGaa                                                                                                               <  5:296606/149‑1 (MQ=255)                                   cAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                              >  5:284312/1‑148 (MQ=255)                                           gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAAcc                                                                                                                                >  5:48533/1‑122 (MQ=255)                                           gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAAcc                                                                                                                                <  6:48533/122‑1 (MQ=255)                                           gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGAAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCCACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCTTCGAACAGGCCGAAGCACAccgc                                                                                                      >  8:146215/1‑148 (MQ=255)                                                ttACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTg                                                                                                <  5:101403/149‑1 (MQ=255)                                                   ccGTGAGTTTCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATc                                                                                             <  8:2950/149‑1 (MQ=255)                                                    cGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCGAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTgg                                                                                                  >  1:387221/1‑143 (MQ=255)                                                    cGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCAAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTgg                                                                                                  <  2:387221/143‑1 (MQ=255)                                                     gTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTg                                                                                           <  2:39248/149‑1 (MQ=255)                                                        aGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGaa                                                                                         <  8:325524/148‑1 (MQ=255)                                                         gTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAAc                                                                                        <  3:273444/148‑1 (MQ=255)                                                           aTCGGAGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGAAGAACTGCTGGACAGTGACGACCAGGCGCTGTA‑CCGGGCCCAGGCCCGGCGCCGCAACCGCGACGAACAGGCAGAAGACAACCGCCTGGTGATCTGAACTAg                                                                                     >  1:306635/1‑149 (MQ=255)                                                            tCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAg                                                                                     <  4:362962/148‑1 (MQ=255)                                                              ggTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGcc                                                                                   <  6:356512/148‑1 (MQ=255)                                                                  tgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCca                                                                              <  4:70817/149‑1 (MQ=255)                                                                    gCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCcag                                                                             >  6:306849/1‑148 (MQ=255)                                                                        gCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGt                                                                        >  8:9456/1‑149 (MQ=255)                                                                             gACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCACCCGCGTCGAACAGGCCGAAGCACACCGCCTGGGGATCTGAACTAGCCGGCCAGCCGGTCCCAc                                                                   >  8:195597/1‑149 (MQ=255)                                                                               cGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                 >  4:181689/1‑149 (MQ=255)                                                                                gCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGTCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                 >  2:69936/1‑148 (MQ=255)                                                                                 cAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACa                                                               <  7:9456/149‑1 (MQ=255)                                                                                  aGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAg                                                              <  7:34497/149‑1 (MQ=255)                                                                                    aCATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCt                                                            <  3:74594/149‑1 (MQ=255)                                                                                     cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCtt                                                           <  4:336375/149‑1 (MQ=255)                                                                                      aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                          <  1:251715/149‑1 (MQ=255)                                                                                           ggCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGCCCAGCAGGTCCCACGACAGCTTGGCGAt                                                     >  4:169498/1‑149 (MQ=255)                                                                                            gCCCCCAGCCTCGACGAACTGCTGGCCAGTGCGGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATc                                                    <  8:60960/149‑1 (MQ=255)                                                                                                   gCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcacac                                             <  5:212710/149‑1 (MQ=255)                                                                                                    aatCGACGAACTGCTGGCCAGTGCCGACCAGGCTCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCTCGTCTAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcacaca                                            <  1:260680/147‑1 (MQ=255)                                                                                                        gacgaACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACacaac                                        <  7:210230/149‑1 (MQ=255)                                                                                                          cgaACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAg                                      >  3:270144/1‑149 (MQ=255)                                                                                                           gaACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGc                                     <  3:169498/149‑1 (MQ=255)                                                                                                            aaCTGCTGGCCAGTGCAGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGTGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCa                                    <  8:339311/149‑1 (MQ=255)                                                                                                                 ctgGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGTGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                <  7:146215/148‑1 (MQ=255)                                                                                                                 ctgGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                <  5:306849/148‑1 (MQ=255)                                                                                                                      cAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCa                                 <  8:28841/142‑1 (MQ=255)                                                                                                                      cAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCa                                 >  7:28841/1‑142 (MQ=255)                                                                                                                       aGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGc                         <  8:303704/149‑1 (MQ=255)                                                                                                                            cGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACCCACAACAGCACCAGGAACAGCCcgcg                    >  5:108502/1‑149 (MQ=255)                                                                                                                               ccAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGa                  >  5:18922/1‑148 (MQ=255)                                                                                                                                         tA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCa       >  1:287777/1‑149 (MQ=255)                                                                                                                                          ccccGGGCCTTGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCa       >  5:363280/3‑149 (MQ=255)                                                                                                                                              gggCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCt   >  7:344240/1‑149 (MQ=255)                                                                                                                                               ggCCAAGGCGCGGGGCCGCTACCGCGTCGAACATGGCGATGCACAGCGCCTTGTGATCTGAACAAGCCGGCCAGCAGGTCCGTCGACAGCTTGGCGATCAGCACACACAACAGCACCAGTAACAGCCCGCGTACAAAACCTGCACCCtt  <  2:306635/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTA‑CCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGCCTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTT  >  NC_002947/4438064‑4438353

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