Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,238	A→T	7.4%	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.4%	195.9 / 4.6	67	L42H (CTC→CAC)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base T (2/3);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297093‑5297385                                                                                                                                                   |                                                                                                                                                    cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                                                                                  >  2:243145/1‑149 (MQ=255) cTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCa                                                                                                                                                  >  6:212924/1‑149 (MQ=255)  tCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCTGGGCGTGctc                                                                                                                                                 >  2:2815/1‑147 (MQ=255)     gCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                                                                                                                              >  5:83368/1‑149 (MQ=255)       cACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                                                                                                                                            >  2:451/1‑149 (MQ=255)         cttcttCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGTGCTCGCCctgc                                                                                                                                          <  7:51579/149‑1 (MQ=255)            cttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagca                                                                                                                                         <  5:204703/147‑1 (MQ=255)             ttcGATACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAggg                                                                                                                                      >  2:285885/1‑149 (MQ=255)                ggTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCTGGGCGTGCTCGCCctgcgtgtct                                                                                                                                    <  1:2815/148‑7 (MQ=255)                 gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                  >  3:253742/1‑149 (MQ=255)                  tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 <  6:231645/149‑1 (MQ=255)                  tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 <  5:212924/149‑1 (MQ=255)                    cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGga                                                                                                                               <  5:353302/149‑1 (MQ=255)                     tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 >  5:195802/1‑146 (MQ=255)                     tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 >  4:55235/1‑146 (MQ=255)                     tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                                                                 >  2:68908/1‑146 (MQ=255)                       gCAATCACCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCAG‑CTTCTACGGTGTAGTTGATCGCGCGGGCGATCAACTGCAGGTGATGGTTGGCAGAGACACGCAGGGCGCGCAGGCGCAGCAGGGTGAGGgaga                                                                                                                             >  8:167397/1‑148 (MQ=255)                         aaTCAACGCCAGCGATTCGTTGGCGCCGGCGCAGTCATCGAACACCACAGCCGCCTTCTACGGTCTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGATCACGACCAGCAGGGAGAGGGAGATa                                                                                                                           >  1:290243/1‑149 (MQ=255)                               cGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCa                                                                                                                    <  1:68908/149‑1 (MQ=255)                                gCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg                                                                                                                   <  5:369865/149‑1 (MQ=255)                                  cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                                                                 <  3:55235/149‑1 (MQ=255)                                      gATTCGTTGGCGGCGTCGCTGGCATCGAACACGACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAAc                                                                                                             <  6:251282/149‑1 (MQ=255)                                       aTTCGTTGGCGGCGGCGCTGGCCGCGAACCCCACAGCCG‑CGTCTACGGGGTAGGTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                                                                            <  6:131407/149‑1 (MQ=255)                                        ttCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACgg                                                                                                           >  1:241944/1‑149 (MQ=255)                                          cGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgc                                                                                                          <  4:253742/148‑1 (MQ=255)                                               gcggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc                                                                                                                      >  1:228084/1‑131 (MQ=255)                                               gcggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc                                                                                                                      <  2:228084/131‑1 (MQ=255)                                                 ggcggcgCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAgg                                                                                                  <  1:285885/149‑1 (MQ=255)                                                         tGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCt                                                                                           >  3:2208/1‑148 (MQ=255)                                                           tCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        >  4:208647/1‑149 (MQ=255)                                                           tCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAg                                                                                        <  1:243145/149‑1 (MQ=255)                                                             aTCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                                                                      <  3:178885/149‑1 (MQ=255)                                                              tCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTg                                                                                     <  4:2208/149‑1 (MQ=255)                                                                gAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTgg                                                                                    <  1:451/148‑1 (MQ=255)                                                                      cacaGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGCTACCCGCCAGGc                                                                                                                 >  3:136984/1‑113 (MQ=255)                                                                        caGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                                                                       <  6:266218/89‑1 (MQ=255)                                                                        caGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAAGCCCAGCAgg                                                                                                                                       >  5:266218/1‑89 (MQ=255)                                                                                  cTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGa                                                                   <  6:195802/148‑1 (MQ=255)                                                                                        tgtAGTTGATCGAGCGGGGGATCAACTGCTTGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCGCAGCAGGGTGAGGGAGATACCCGCCTGGCCAACGGCAAGCAAGGCGACGCTGAGTTGGCCACGGCCCTGCACCGAGCGCAgg                                                            <  2:290243/149‑1 (MQ=255)                                                                                        tgcaGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAgg                                                            >  6:153701/4‑149 (MQ=255)                                                                                            gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                                                                   >  1:47529/1‑106 (MQ=255)                                                                                            gTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGcaagcaag                                                                                                   <  2:47529/106‑1 (MQ=255)                                                                                             ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                  >  5:22409/1‑74 (MQ=38)                                                                                             ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                                                                                                                                  <  6:22409/74‑1 (MQ=38)                                                                                                aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc                                                     <  6:83368/148‑1 (MQ=255)                                                                                                         gCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACg                                           <  6:190376/149‑1 (MQ=255)                                                                                                             tCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTg                                       >  2:145227/1‑149 (MQ=255)                                                                                                              cAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGc                                      >  7:134039/1‑149 (MQ=255)                                                                                                                    cAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTc                                 <  3:208647/148‑1 (MQ=255)                                                                                                                      ggTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGc                              >  2:286431/1‑149 (MQ=255)                                                                                                                       gTTATGGTTGGCATAGACACGCAGGGCGAGGACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGGAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCgg                               <  6:10691/147‑1 (MQ=255)                                                                                                                       gTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCt                             >  5:77147/1‑149 (MQ=255)                                                                                                                        ttATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCtt                            >  3:154928/1‑149 (MQ=255)                                                                                                                              ttGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgc                      <  2:241944/149‑1 (MQ=255)                                                                                                                                ggCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgct                     <  5:153701/148‑1 (MQ=255)                                                                                                                                   aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTc                 >  7:2524/1‑149 (MQ=255)                                                                                                                                   aTAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTc                 >  6:277652/1‑149 (MQ=255)                                                                                                                                         ctgctgggCGTGCTCGCCCTGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgc                      >  7:346219/6‑138 (MQ=255)                                                                                                                                         ctgctgggCGTGCTCGCCCTGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgc                      <  8:346219/133‑1 (MQ=255)                                                                                                                                         acGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCg           >  6:218275/1‑149 (MQ=255)                                                                                                                                           gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgctgct                  <  4:135925/140‑1 (MQ=255)                                                                                                                                           gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGgctgct                  >  3:135925/1‑140 (MQ=255)                                                                                                                                           gCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgc         >  6:177749/1‑149 (MQ=255)                                                                                                                                             aGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGgcgc       >  5:215788/1‑149 (MQ=255)                                                                                                                                                gCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCgg          <  6:152928/143‑1 (MQ=255)                                                                                                                                                gCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCgg          >  5:152928/1‑143 (MQ=255)                                                                                                                                                  gAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGa  >  5:322401/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCG‑CTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGA  >  NC_002947/5297093‑5297385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑