Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,057,006	Δ1 bp	9.8%	coding (385/990 nt)	PP_2669 ←	outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,057,006	0	T	.	9.8%	236.3 / 21.3	82	coding (385/990 nt)	PP_2669	outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/33);  new base . (4/4);  total (45/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 = 3.96e-01

TCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGC  >  NC_002947/3056862‑3057155                                                                                                                                                  |                                                                                                                                                       tCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTAGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCg                                                                                                                                                   >  1:1720354/1‑150 (MQ=255)   tGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTc                                                                                                                                                 <  1:1225942/150‑1 (MQ=255)   tGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTc                                                                                                                                                 >  2:686225/1‑150 (MQ=255)    gggTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCa                                                                                                                                                <  1:783862/150‑1 (MQ=255)        gCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCTCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAg                                                                                                                                            <  2:274356/150‑1 (MQ=255)           ggTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAgcgc                                                                                                                                         <  1:686225/150‑1 (MQ=255)            gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAccgacaag                                                                                                                                        <  2:1828251/150‑7 (MQ=255)            gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCa                                                                                                                                                >  2:956752/1‑142 (MQ=255)            gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCa                                                                                                                                                <  1:956752/142‑1 (MQ=255)              gtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATc                                                                                                                                      <  1:375975/150‑1 (MQ=255)                 gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGGCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑c                                                                                                                                  <  2:550385/150‑1 (MQ=255)                 gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAg                                                                                                                                            >  2:1638709/1‑141 (MQ=255)                 gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAg                                                                                                                                            <  1:1638709/141‑1 (MQ=255)                   tCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCt                                                                                                                                >  1:638461/1‑150 (MQ=255)                       gTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGtt                                                                                                                            <  2:620510/150‑1 (MQ=255)                        tGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGttt                                                                                                                           >  2:833489/1‑150 (MQ=255)                             caTGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGatac                                                                                                                      >  1:1762100/1‑150 (MQ=255)                             caTGTTGGTGGTTTCGCTAGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGatac                                                                                                                      <  2:1720354/150‑1 (MQ=255)                                            gCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTc                                                                                                       >  1:1726685/1‑150 (MQ=255)                                            gCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTc                                                                                                       >  2:269837/1‑150 (MQ=255)                                            gCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACCTCGATTTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTc                                                                                                       >  2:1617028/1‑150 (MQ=255)                                             cTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCg                                                                                                      >  1:1926032/1‑150 (MQ=255)                                              tGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑CCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGt                                                                                                     >  1:1316214/1‑150 (MQ=255)                                                gtgtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTg                                                                                                   >  2:674554/1‑150 (MQ=255)                                                 tgtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCTTTTGATACATACAGCCAGCGGTCGTTgg                                                                                                  <  2:638461/150‑1 (MQ=255)                                                 tgtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGttt                                                                                                                           <  1:771587/125‑1 (MQ=255)                                                 tgtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGttt                                                                                                                           >  2:771587/1‑125 (MQ=255)                                                  gtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                                 <  2:208280/150‑1 (MQ=255)                                                   ttGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAcc                                                                                                                                              <  2:1499659/105‑1 (MQ=255)                                                   ttGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAcc                                                                                                                                              >  1:1499659/1‑105 (MQ=255)                                                        cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                           >  1:192134/1‑150 (MQ=255)                                                        cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                           >  2:1682039/1‑150 (MQ=255)                                                         gCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAgg                                                                                          >  1:334712/1‑150 (MQ=255)                                                          cccACTTGCCGTCGGGGCTGACCGCCATGCCTTGGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAggg                                                                                         >  2:178312/1‑150 (MQ=255)                                                           ccACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAcc                                                                                                                                              <  2:454010/97‑1 (MQ=39)                                                           ccACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAcc                                                                                                                                              >  1:454010/1‑97 (MQ=39)                                                           ccACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTAAACCACACATCGTAACTCCAGT‑AAACATACAGCCAGCGGTCGTTGGGGTGCAGGGc                                                                                        >  1:1398862/1‑150 (MQ=255)                                                            cACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCg                                                                                       >  2:706996/1‑150 (MQ=255)                                                              cTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGTTTCGT‑CCTCGTTTGATACATACAGCCAGTGGTCGTTGGGGTGCAGGGCGaa                                                                                     <  2:758548/150‑1 (MQ=255)                                                                    gTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTc                                                                               >  1:504836/1‑150 (MQ=255)                                                                      cGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCgg                                                                             <  1:716653/150‑1 (MQ=255)                                                                       ggggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                            <  1:257599/150‑1 (MQ=255)                                                                       ggggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                            <  2:1267202/150‑1 (MQ=255)                                                                       ggggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGAGGAAGTTGGGGTGACGGGCGAACAGCttcgg                                                                            >  2:900677/1‑146 (MQ=255)                                                                        gggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑CCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCgggg                                                                           <  2:1316214/150‑1 (MQ=255)                                                                                cGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                                   >  2:1110371/1‑150 (MQ=255)                                                                                cGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                                   >  2:1620893/1‑150 (MQ=255)                                                                                 gCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACGAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCg                                                                  >  1:1401382/1‑150 (MQ=255)                                                                                         ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCaa                                                          >  2:820853/1‑150 (MQ=255)                                                                                               ttCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtt                                                    <  1:674554/150‑1 (MQ=255)                                                                                                    tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                             >  1:1642780/1‑136 (MQ=255)                                                                                                    tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                             <  2:1642780/136‑1 (MQ=255)                                                                                                       cGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCt                                            >  1:379037/1‑150 (MQ=255)                                                                                                         aCTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGCTGATCttt                                          >  2:862518/1‑149 (MQ=255)                                                                                                              gATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGt                                     <  2:1395478/150‑1 (MQ=255)                                                                                                              gATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGt                                     <  1:820853/150‑1 (MQ=255)                                                                                                               aTCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTg                                    >  2:1497611/1‑150 (MQ=255)                                                                                                                tCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGAAGGGCGAACTGCTCGGGGTACTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTgg                                   >  2:1852323/1‑150 (MQ=255)                                                                                                                      cGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGt                             >  1:1588342/1‑150 (MQ=255)                                                                                                                        aGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTcc                           <  1:706996/150‑1 (MQ=255)                                                                                                                          cACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCAt                         <  2:192134/150‑1 (MQ=255)                                                                                                                            ccTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGGGGCCACGTCCATCa                       >  2:740001/1‑150 (MQ=255)                                                                                                                             tgcgctggTGACGGTGATCGA‑CACCGTCACCGGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCAc                      >  2:749854/7‑150 (MQ=255)                                                                                                                              ttGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGttt                                                                                                                           <  2:1562777/48‑1 (MQ=255)                                                                                                                              ttGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGttt                                                                                                                           >  1:1562777/1‑48 (MQ=255)                                                                                                                              ttGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtga                                                  <  1:1315406/121‑1 (MQ=255)                                                                                                                              ttGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtga                                                  >  2:1315406/1‑121 (MQ=255)                                                                                                                                  cGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTCATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCa                 <  2:1398862/150‑1 (MQ=255)                                                                                                                                   ggTGACGGTGATCGA‑CACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGAGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAc                >  1:1053507/1‑150 (MQ=255)                                                                                                                                    gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAcc               >  2:592865/1‑150 (MQ=255)                                                                                                                                     ggACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAccc              <  2:334712/149‑1 (MQ=255)                                                                                                                                       aCGGCG‑TAGCTCAACGTCACCAACGCATCGG‑TCTCGTTTGATACATACAGCCAGCGGTCGTTGGTGGGCATGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCgg            <  1:1490789/150‑1 (MQ=255)                                                                                                                                        cGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACCTACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGt           <  2:1514203/150‑1 (MQ=255)                                                                                                                                         ggTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGc                                         <  1:505761/119‑1 (MQ=255)                                                                                                                                         ggTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGc                                         >  2:505761/1‑119 (MQ=255)                                                                                                                                          gtg‑tCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                             >  1:1085503/1‑98 (MQ=255)                                                                                                                                          gtg‑tCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                             <  2:1085503/98‑1 (MQ=255)                                                                                                                                          gtg‑tCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc         <  1:1617028/150‑1 (MQ=255)                                                                                                                                          gtg‑tCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc         <  1:1682039/150‑1 (MQ=255)                                                                                                                                              tCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGaa      >  2:453614/1‑150 (MQ=255)                                                                                                                                                gATCACCGTCACCAGCGCATCGT‑GCTCGTTTGATACATACAGCCAGCGGTCGGTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTTCCCCCGGTCCGAATc    >  2:391553/1‑150 (MQ=255)                                                                                                                                                gATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGAGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATc    >  2:1393865/1‑150 (MQ=255)                                                                                                                                                  tCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGc  <  1:1620893/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       TCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGT‑CCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGC  >  NC_002947/3056862‑3057155

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