Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,846,297	Δ1 bp	5.7%	coding (588/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,297	0	T	.	5.7%	214.6 / 7.2	70	coding (588/1377 nt)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base . (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAAT  >  NC_002947/1846152‑1846443                                                                                                                                                  |                                                                                                                                                    cGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAg                                                                                                                                                  <  4:270503/149‑1 (MQ=255)       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)                                                                                                                                                  |                                                                                                                                                    CGGGCCAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCGCATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAAT  >  NC_002947/1846152‑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: ‑