Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,238	A→T	7.1%	L42H (CTC→CAC)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,238	0	A	T	7.1%	204.2 / 5.0	70	L42H (CTC→CAC)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (38/27);  new base T (3/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297091‑5297385                                                                                                                                                    |                                                                                                                                                     gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                                                                                                                      <  8:158376/148‑1 (MQ=255)    tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCTGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                  >  2:234184/1‑149 (MQ=255)      ggCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                                                                                                                                  <  7:225608/147‑1 (MQ=255)          aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                                                                                                                            >  4:131202/1‑149 (MQ=255)             tcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagcag                                                                                                                                         >  2:179895/1‑149 (MQ=255)                tcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGt                                                                                                                                      >  8:284980/1‑149 (MQ=255)                 cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                                                                        <  1:35737/146‑1 (MQ=255)                 cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGt                                                                                                                                      <  3:131202/148‑1 (MQ=255)                   gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                   >  5:126640/1‑149 (MQ=255)                      cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAggg                                                                                                                                 <  6:161063/148‑1 (MQ=255)                         gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgaga                                                                                                                              >  6:98656/1‑148 (MQ=255)                         gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCACAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGAt                                                                                                                             <  5:128415/149‑1 (MQ=255)                           aaTCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAc                                                                                                                           >  8:272469/1‑149 (MQ=255)                                   ccAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑gg                                                                                                                  >  8:192063/1‑149 (MQ=255)                                         aTTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACg                                                                                                            <  8:67445/149‑1 (MQ=255)                                              ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGgcaag                                                                                                       <  8:228555/149‑1 (MQ=255)                                                             tCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  6:126640/149‑1 (MQ=255)                                                               aTCGAACACCACAGCCGCTTCAACGCTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTAATGGATGGCAAAGACACGCAGGGCGAGCACGCCAAGCAGGGTCAGGCAGATACCAGCCA‑CGACAACGCCCAGCAAGGCGACGCTGAGGt                                                                                      >  6:254853/1‑149 (MQ=255)                                                                      accacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACg                                                                               >  8:36138/1‑149 (MQ=255)                                                                        cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                                                                             >  2:115253/1‑149 (MQ=255)                                                                          caGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGcc                                                                            >  3:208327/1‑148 (MQ=255)                                                                           aGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCa                                                                          <  4:208327/149‑1 (MQ=255)                                                                               gCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGccc                                                                           >  8:54582/1‑144 (MQ=255)                                                                               gCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGccc                                                                           <  7:54582/144‑1 (MQ=255)                                                                                   cTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCTGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCg                                                                    <  1:234184/147‑1 (MQ=255)                                                                                   cTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAg                                                                  >  4:20916/1‑149 (MQ=255)                                                                                     aacctgcaGTTGATCGCCCGCTCGATCAACTACAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAgcg                                                                >  7:158393/8‑149 (MQ=255)                                                                                        gtgtAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGGATAGACTCGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAg                                                             >  3:21853/1‑149 (MQ=255)                                                                                           caGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGt                                                          >  4:17747/2‑149 (MQ=255)                                                                                           caGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGt                                                          >  5:239877/2‑149 (MQ=255)                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTc                                                         <  7:272469/148‑1 (MQ=255)                                                                                             gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTc                                                         <  1:115253/148‑1 (MQ=255)                                                                                               tGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg                                                       <  7:284980/148‑1 (MQ=255)                                                                                                gATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc                                                     >  5:230492/1‑149 (MQ=255)                                                                                                 aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgcg                                                    >  8:54382/1‑149 (MQ=255)                                                                                                  tCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgcg                                                    >  7:240637/1‑148 (MQ=255)                                                                                                          gCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACg                                           >  7:99596/1‑149 (MQ=255)                                                                                                            gATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCt                                         >  3:139604/1‑149 (MQ=255)                                                                                                              tCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTg                                       >  3:278969/1‑149 (MQ=255)                                                                                                                aaCTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGc                                      <  8:99596/148‑1 (MQ=255)                                                                                                                aaCGGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGc                                      <  4:12672/148‑1 (MQ=255)                                                                                                                 aCTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCgg                                    <  2:42705/149‑1 (MQ=255)                                                                                                                  cTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACgg                                               >  5:67779/1‑137 (MQ=255)                                                                                                                  cTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACgg                                               <  6:67779/137‑1 (MQ=255)                                                                                                                        gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGGTCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGc                              >  1:246101/1‑149 (MQ=255)                                                                                                                        gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACg                                                <  8:257284/130‑1 (MQ=255)                                                                                                                        gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACg                                                >  7:257284/1‑130 (MQ=255)                                                                                                                        gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCt                             >  2:250808/1‑149 (MQ=255)                                                                                                                         ttATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  2:124459/1‑149 (MQ=255)                                                                                                                           aTGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTc                           <  7:36138/148‑1 (MQ=255)                                                                                                                            tGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCAt                         <  6:239877/149‑1 (MQ=255)                                                                                                                             ggTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATg                        >  1:153365/1‑149 (MQ=255)                                                                                                                             ggTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATg                        >  7:100393/1‑149 (MQ=255)                                                                                                                                  gCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgctgc                   >  3:110021/1‑149 (MQ=255)                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGTGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            <  4:18257/138‑1 (MQ=255)                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  7:45149/1‑138 (MQ=255)                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            <  8:45149/138‑1 (MQ=255)                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  3:18257/1‑138 (MQ=255)                                                                                                                                    aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTc                 <  8:158393/149‑1 (MQ=255)                                                                                                                                      aGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCAt                         <  7:131135/139‑1 (MQ=255)                                                                                                                                      aGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACCGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCAt                         >  8:131135/1‑139 (MQ=255)                                                                                                                                       gACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCg               <  4:139604/148‑1 (MQ=255)                                                                                                                                         cacGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGcc             <  4:21853/148‑1 (MQ=255)                                                                                                                                          ctgctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCa                          <  5:122573/129‑1 (MQ=255)                                                                                                                                          ctgctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCa                          >  6:122573/6‑134 (MQ=255)                                                                                                                                          ctgctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCg           >  7:96267/6‑149 (MQ=255)                                                                                                                                           tgctgggCGTGCTCGCCCTGCAGGGTGAGGGAGATACCCGCCC‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCgg          >  3:302464/5‑149 (MQ=255)                                                                                                                                            gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         >  2:237572/1‑149 (MQ=255)                                                                                                                                            gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         >  4:35052/1‑149 (MQ=255)                                                                                                                                             ctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         <  8:96267/146‑1 (MQ=255)                                                                                                                                                   gAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGa  >  4:158802/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297091‑5297385

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