Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,846,303:1	+A	5.8%	coding (594/1377 nt)	pfeS‑II →	histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,846,303	1	.	A	5.8%	279.1 / 11.3	69	coding (594/1377 nt)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base . (34/31);  new base A (2/2);  total (36/33)
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.92e-01

AGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAAT  >  NC_002947/1846158‑1846443                                                                                                                                                   |                                                                                                                                             aGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCg                                                                                                                                           >  3:317113/1‑149 (MQ=255)     ccGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGccc                                                                                                                                       >  4:333027/1‑149 (MQ=255)             tgtgGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCGG‑TCGGCCCGCCTGGa                                                                                                                                >  5:245581/1‑148 (MQ=255)                        gTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCa                                                                                                                     <  3:333027/148‑1 (MQ=255)                         tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATc                                                                                                                   <  4:333156/149‑1 (MQ=255)                          ctccAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCg                                                                                                                  <  5:231630/147‑1 (MQ=255)                           aCCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCcagttcgc                                                                                                                 >  5:51656‑M1/1‑141 (MQ=255)                             cAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCg                                                                                                                  >  6:231630/1‑146 (MQ=255)                                 ggCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCa                                                                                                           >  8:175937/1‑149 (MQ=255)                                  gCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCAt                                                                                                          <  5:25878/149‑1 (MQ=255)                                          ccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCAGCCTGGACTCGCGCACCATCGCCCGGCAGGACGAAtt                                                                                                  >  4:8720/1‑149 (MQ=255)                                                 tGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGaa                                                                                           <  7:175937/149‑1 (MQ=255)                                                   ttCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGTGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                                         >  4:63601/1‑149 (MQ=255)                                                   ttCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACt                                                                                         >  3:27630/1‑149 (MQ=255)                                                      aCCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgc                                                                                      <  5:335054‑M2/149‑3 (MQ=255)                                                        cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAACTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGc                                                                                              >  2:77329/1‑139 (MQ=255)                                                        cTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAACTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGc                                                                                              <  1:77329/139‑1 (MQ=255)                                                           gCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgcgcg                                                                                   >  6:245581‑M2/1‑142 (MQ=255)                                                           gCTGTGGTGTGTTGGCTTGTACCGCCGGTGTATCGTGCCGCTGAACCAACTTCGCGAGCTGGCCATTGCCTGGGGGGCCGATCTGGTG‑TCGGCCCGCCTGGACTCGCGCcagttcgcccaattcgtcatgccgggcgatggtgcgcga                                                                                  <  6:51656‑M1/148‑40 (MQ=255)                                                            cTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGc                                                                                <  4:27630/149‑1 (MQ=255)                                                              gttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                               >  2:298015/1‑132 (MQ=255)                                                              gttgtGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                               <  1:298015/132‑1 (MQ=255)                                                                 gtgtgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCttt                                                                            <  4:317113/148‑1 (MQ=255)                                                                  tgtgtAGGCCTGTACCGCATGCTGATAGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGa                                                                          >  1:5224/1‑149 (MQ=255)                                                                    tgtAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                         >  2:247480/1‑148 (MQ=255)                                                                         gCCTGTACCGCATGCTGATCGTGCCGCTGACCCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCcagttcgcccaattcgtcatgccgggcgatggtgcgcgagtccag                                                                            <  7:294792‑M1/140‑46 (MQ=255)                                                                         gCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCcagttcgcccaattcgtcatgccgggcgatggtgcgcgagtccag                                                                            >  8:294792‑M1/1‑95 (MQ=255)                                                                            tGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGcc                                                                >  5:362736/1‑149 (MQ=255)                                                                              tACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                         <  2:35110/138‑1 (MQ=255)                                                                              tACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGAc                                                                         >  1:35110/1‑138 (MQ=255)                                                                                 cGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCg                                                           >  4:170471/1‑149 (MQ=255)                                                                                      gCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCttt                                                                            <  2:50603/127‑1 (MQ=255)                                                                                      gCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCttt                                                                            >  1:50603/1‑127 (MQ=255)                                                                                          aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAgg                                                  <  3:170471/149‑1 (MQ=255)                                                                                          aTCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAgg                                                  >  8:131285/1‑149 (MQ=255)                                                                                               gCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACgg                                             >  2:129256/1‑149 (MQ=255)                                                                                                      aaCCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGt                                                                                  >  4:286033/1‑105 (MQ=255)                                                                                                       aCCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATg                                     <  2:5224/149‑1 (MQ=255)                                                                                                        ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAg                                                   >  6:340223/1‑134 (MQ=255)                                                                                                        ccAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAg                                                   <  5:340223/134‑1 (MQ=255)                                                                                                         cAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTgg                                          >  1:240338/1‑142 (MQ=255)                                                                                                         cAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTgg                                          <  2:240338/142‑1 (MQ=255)                                                                                                          aaCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcag                                  <  1:247480/149‑1 (MQ=255)                                                                                                           aCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTg                                           <  1:314104/139‑1 (MQ=255)                                                                                                           aCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTg                                           >  2:314104/1‑139 (MQ=255)                                                                                                           aCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGAGGCTGCAAGGCACGGTGGCTATgcag                                  >  5:175354/1‑148 (MQ=255)                                                                                                                 cgAGGAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                           <  3:8720/149‑1 (MQ=255)                                                                                                                  gAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGtt                          >  2:265744/1‑149 (MQ=255)                                                                                                                     cacgCCCATGCCTGGCGGGCCGATCAGCTG‑GCGGCCCGCCTGGACTCGCGCACCATCGCCCCGCATGACGAATTCGCCGAACTCGCCCGTGCCATTCACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgct                       <  4:251330/146‑1 (MQ=255)                                                                                                                      aGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctg                      >  1:56502/1‑149 (MQ=255)                                                                                                                        gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCTCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCaa                                                    <  8:143825/117‑1 (MQ=255)                                                                                                                        gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCACGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCaa                                                    >  7:143825/1‑117 (MQ=255)                                                                                                                           aaTGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGa                 <  6:362736/149‑1 (MQ=255)                                                                                                                                cTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                           <  1:210858/134‑1 (MQ=255)                                                                                                                                cTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                           >  5:251868/1‑134 (MQ=255)                                                                                                                                cTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                           <  6:251868/134‑1 (MQ=255)                                                                                                                                cTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                           >  2:210858/1‑134 (MQ=255)                                                                                                                                cTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGt            <  3:63601/149‑1 (MQ=255)                                                                                                                                 aggCGGGCCGA‑CAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCATGGCACGGTg                                           >  5:305789/2‑117 (MQ=255)                                                                                                                                 aggCGGGCCGA‑CAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTg                                           <  6:305789/116‑1 (MQ=255)                                                                                                                                 aggCGGGCCGA‑CAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTc           <  4:84028/148‑1 (MQ=255)                                                                                                                                 aggCGGGCCGA‑CAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTc           >  3:84028/2‑149 (MQ=255)                                                                                                                                   gCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACgg                                             <  2:270013/113‑1 (MQ=255)                                                                                                                                   gCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACgg                                             >  1:270013/1‑113 (MQ=255)                                                                                                                                   gCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTc           <  6:335054/147‑1 (MQ=255)                                                                                                                                      ggCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTg                                                       <  8:110405/100‑1 (MQ=255)                                                                                                                                      ggCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTg                                                       >  7:110405/1‑100 (MQ=255)                                                                                                                                         cGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGAGACCTGTCCCACGaaa   >  7:308296/1‑149 (MQ=255)                                                                                                                                          gATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAAt  <  2:56502/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAAT  >  NC_002947/1846158‑1846443

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