Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,057,006	Δ1 bp	6.3%	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	.	6.3%	272.9 / 13.4	95	coding (385/990 nt)	PP_2669	outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base . (2/4);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGCT  >  NC_002947/3056860‑3057156                                                                                                                                                    |                                                                                                                                                       ggTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCAc                                                                                                                                                     >  1:841839/1‑150 (MQ=255)  gTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGGCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCAcc                                                                                                                                                    <  1:596562/150‑1 (MQ=255)      gggTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTGGATGTCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                                <  2:833990/150‑1 (MQ=255)       ggTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAc                                                                                                                                               <  2:1758211/150‑1 (MQ=255)         tGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGACCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCa                                                                                                                                             >  1:1900540/1‑150 (MQ=255)            tGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAgcg                                                                                                                                          >  2:1566332/1‑150 (MQ=255)              gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCa                                                                                                                                        >  2:1355509/1‑150 (MQ=255)              gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTGCGATCTGCCCGAGCACCTTGTCGGTGATGGTG‑TCGAGCACCGTCACCAGCGCa                                                                                                                                        >  2:907443/1‑150 (MQ=255)               tgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCAt                                                                                                                                       <  2:269686/150‑1 (MQ=255)                   gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGCCCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTc                                                                                                                                   <  2:137861/150‑1 (MQ=255)                   gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCCACAGCGCATCGTc                                                                                                                                   >  1:1433457/1‑150 (MQ=255)                   gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTc                                                                                                                                   >  1:1921104/1‑150 (MQ=255)                      ccAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTc                                                                                                                                <  1:47313/150‑1 (MQ=255)                      ccAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTc                                                                                                                                >  2:1395556/1‑150 (MQ=255)                       cAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCg                                                                                                                               >  1:427686/1‑150 (MQ=255)                          tGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGttt                                                                                                                            <  1:70150/150‑1 (MQ=255)                             agcaTGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGat                                                                                                                         <  1:549235/150‑1 (MQ=255)                              gcaGGTTGGTGGGTTCGCTGGTGTTGACCGACCCCTTGCCTGGTGTGCGGTCCGCCATGCCTTCGGGTTCCTGGCGGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGata                                                                                                                        <  1:1374263/150‑1 (MQ=255)                                  gTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                                    >  1:1899090/1‑150 (MQ=255)                                    tggtggTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGAtacataca                                                                                                                  >  2:734428/1‑150 (MQ=255)                                         gtTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATAcagccag                                                                                                             >  2:1270782/1‑150 (MQ=255)                                          tttCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATc                                                                                                                                      >  2:46432/1‑124 (MQ=255)                                          tttCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATc                                                                                                                                      <  1:46432/124‑1 (MQ=255)                                                tGGTGTTGACCGCCCACTTGCCGACGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGGGCACCTTGTCGGTGACGGTG‑TCGATCTCCGTAACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGt                                                                                                      >  1:911809/1‑150 (MQ=255)                                                    gttgcgcGCCCACTTGGCGTCGGGGCTGACCGCCAGGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGGCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                                  <  2:293179/144‑1 (MQ=255)                                                    gtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                                  <  1:1355509/150‑1 (MQ=255)                                                    gtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                                  <  2:427686/150‑1 (MQ=255)                                                     ttGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTgggg                                                                                                 >  1:1541833/1‑150 (MQ=255)                                                          cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACGCCCAGCGGTCGTTGGGGTGCAg                                                                                            <  1:537752/150‑1 (MQ=255)                                                          cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                            >  2:1246235/1‑150 (MQ=255)                                                          cGCCCACTGGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                            >  1:643387/1‑150 (MQ=255)                                                           gCCCACTTGCCGTCGGTGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCGCCTTGTCGGAGACGGTG‑TCGATCAACTTCACCAGCGCATCGTCCTCGTTTTATACATACAGCCAGCGGTCGTTGGGGTGCAgg                                                                                           >  1:1455117/1‑150 (MQ=255)                                                           gCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAgg                                                                                           <  1:1566332/150‑1 (MQ=255)                                                             ccACTTGCCGTCGGGGCTGACCGCCATGCCTTAGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATTGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGc                                                                                         >  1:1899984/1‑150 (MQ=255)                                                              cACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCg                                                                                        >  2:1724409/1‑150 (MQ=255)                                                               aCTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGa                                                                                       >  1:998717/1‑150 (MQ=255)                                                                cTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                                <  1:1800358/92‑1 (MQ=255)                                                                cTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                                >  2:1800358/1‑92 (MQ=255)                                                                      gTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTc                                                                                >  1:219141/1‑150 (MQ=255)                                                                          gggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCgggg                                                                            <  2:1285746/150‑1 (MQ=255)                                                                           ggCTGACCGCCACGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGGGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGt                                                                           <  1:568332/150‑1 (MQ=255)                                                                             cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                         >  1:165135/1‑150 (MQ=255)                                                                             cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                         <  1:615362/150‑1 (MQ=255)                                                                             cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                         <  2:1906532/150‑1 (MQ=255)                                                                                 ccGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTCTGATACATACAGCCAGCGGTCGTTGGGCTGCAGGGCGAACTGCTCGGGGTCCTTGc                                                                     <  2:998717/150‑1 (MQ=255)                                                                                  cGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTAGGCTGCAGTGCGAACTGCTCGGGGTCCTTGcc                                                                    <  2:841839/150‑1 (MQ=255)                                                                                      aTGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCCTCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGaa                                                                <  2:219141/150‑1 (MQ=255)                                                                                      aTGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCGTGCCGGaa                                                                <  2:78481/150‑1 (MQ=255)                                                                                        gCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                              >  1:1777836/1‑150 (MQ=255)                                                                                         ccTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                             <  1:1518375/150‑1 (MQ=255)                                                                                         ccTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                             <  2:920024/150‑1 (MQ=255)                                                                                           ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                             >  2:1874877/1‑132 (MQ=255)                                                                                           ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                             <  1:1874877/132‑1 (MQ=255)                                                                                           ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCaa                                                           >  1:1775673/1‑150 (MQ=255)                                                                                            tCGGGTTCGATGCCGACTTCGATCTCCCCGAGTATCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAAc                                                          <  2:1455117/150‑1 (MQ=255)                                                                                                 ttCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCTCCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtt                                                     <  1:881462/150‑1 (MQ=255)                                                                                                 ttCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtt                                                     <  1:1395556/150‑1 (MQ=255)                                                                                                   cGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtggg                                                   <  2:656176/150‑4 (MQ=255)                                                                                                   cGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtga                                                   <  2:1899090/150‑1 (MQ=255)                                                                                                     aTGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                             <  2:1043984/138‑1 (MQ=255)                                                                                                     aTGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                             >  1:1043984/1‑138 (MQ=255)                                                                                                     aTGCCGACTTCGATCTGCCCGAGCACCTTGTCAGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtgatg                                                 >  2:1070978/1‑150 (MQ=255)                                                                                                      tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGaa                                                                                      >  2:956970/1‑112 (MQ=255)                                                                                                      tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGaa                                                                                      <  1:956970/112‑1 (MQ=255)                                                                                                         cGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                                    <  1:1533567/127‑1 (MQ=255)                                                                                                         cGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                                    >  2:1533567/1‑127 (MQ=255)                                                                                                           aCTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTg                                           <  1:180398/150‑1 (MQ=255)                                                                                                                      cccGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCAc                                >  1:1521562/1‑150 (MQ=255)                                                                                                                         gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                            >  2:667010/1‑87 (MQ=255)                                                                                                                         gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                            <  1:667010/87‑1 (MQ=255)                                                                                                                         gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTc                             >  2:837773/1‑150 (MQ=255)                                                                                                                                      gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                                    <  2:1350287/50‑1 (MQ=255)                                                                                                                                      gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                                    >  1:1350287/1‑50 (MQ=255)                                                                                                                                      gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTACCGGAAGGCAACTCCTGGATGATCTGGCGTGTGGCCACGTCCATCACCTGCAcc                >  2:1754254/1‑150 (MQ=255)                                                                                                                                       tGACGGTGATCGA‑CACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTg                                                                                                    <  2:547695/65‑1 (MQ=37)                                                                                                                                       tGACGGTGATCGA‑CACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTg                                                                                                    >  1:547695/1‑65 (MQ=37)                                                                                                                                       tGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCATCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAccc               >  1:1251471/1‑150 (MQ=255)                                                                                                                                       tGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAccc               >  2:1643426/1‑150 (MQ=255)                                                                                                                                        gACGGTG‑TCGATCACCGTCACCAGCGCATCGGCCGCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCg              <  1:679776/150‑1 (MQ=255)                                                                                                                                         aCGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCgg             <  2:1541833/150‑1 (MQ=255)                                                                                                                                           ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTc           <  1:875159/150‑1 (MQ=255)                                                                                                                                           ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTc           >  2:182870/1‑150 (MQ=255)                                                                                                                                           ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACATGCACCCGGTc           >  1:1022499/1‑150 (MQ=255)                                                                                                                                            gtg‑tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTTCTTTACGATATGCTTCTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc          <  2:1699262/150‑1 (MQ=255)                                                                                                                                            gtg‑tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc          <  1:182870/150‑1 (MQ=255)                                                                                                                                            gggtgcGATCACCGTCACCGTCGCGTCGTCCTCGTGTCGTGCATCCAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTc           <  2:585835/145‑1 (MQ=255)                                                                                                                                                tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGaa                                                                >  1:1268560/1‑93 (MQ=255)                                                                                                                                                tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGaa                                                                <  2:1268560/93‑1 (MQ=255)                                                                                                                                                tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGACGGCCACTACTTGTTGCTCTTGCGGCTGGCGACGCCCAGCACCCGCACCCGGTCCGaa       >  2:1153689/1‑150 (MQ=255)                                                                                                                                                tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAGCCGGTCCGaa       >  2:802332/1‑150 (MQ=255)                                                                                                                                                  gATCACCGTCACCAGCGCATCGTCCTCGTTTGATTCATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATc     >  2:38924/1‑150 (MQ=255)                                                                                                                                                  gATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATc     >  1:1014051/1‑150 (MQ=255)                                                                                                                                                  gATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTACTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATATCGGGCACCCGGTCCGAATc     >  1:1532764/1‑150 (MQ=255)                                                                                                                                                   aTCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGaa       >  2:1001842/1‑147 (MQ=255)                                                                                                                                                   aTCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGaa       <  1:1001842/147‑1 (MQ=255)                                                                                                                                                    tCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGc   >  1:949342/1‑150 (MQ=255)                                                                                                                                                    tCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGc   <  1:1357484/150‑1 (MQ=255)                                                                                                                                                    tCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCC‑TGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGCt  <  1:1724409/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCCGAATCGCT  >  NC_002947/3056860‑3057156

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