Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,241	A→T	7.2%	V41E (GTG→GAG)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,241	0	A	T	7.2%	198.9 / 5.4	70	V41E (GTG→GAG)	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.03e-01

TCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297094‑5297385                                                                                                                                                    |                                                                                                                                                  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                                                                           <  7:54582/144‑1 (MQ=255)                                                                            gCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGccc                                                                           >  8:54582/1‑144 (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                                                >  7:257284/1‑130 (MQ=255)                                                                                                                     gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACg                                                <  8:257284/130‑1 (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                        >  7:100393/1‑149 (MQ=255)                                                                                                                          ggTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATg                        >  1:153365/1‑149 (MQ=255)                                                                                                                               gCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgctgc                   >  3:110021/1‑149 (MQ=255)                                                                                                                                 aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGTGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            <  4:18257/138‑1 (MQ=255)                                                                                                                                 aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  3:18257/1‑138 (MQ=255)                                                                                                                                 aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            <  8:45149/138‑1 (MQ=255)                                                                                                                                 aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  7:45149/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         >  4:35052/1‑149 (MQ=255)                                                                                                                                         gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         >  2:237572/1‑149 (MQ=255)                                                                                                                                          ctgggCGTGCTCGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         <  8:96267/146‑1 (MQ=255)                                                                                                                                                gAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGa  >  4:158802/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCA‑GGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297094‑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: ‑