Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,286,678	G→C	8.9%	A353G (GCC→GGC)	PP_3758 ←	two component system response regulator/sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,286,678	0	G	C	8.9%	207.2 / 5.9	67	A353G (GCC→GGC)	PP_3758	two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base C (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAAC  >  NC_002947/4286536‑4286825                                                                                                                                                  |                                                                                                                                                    gTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCctgct                                                                                                                                               <  3:292637/149‑1 (MQ=255)  ttttGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCctgctg                                                                                                                                              <  7:131503/149‑1 (MQ=255)   tttGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACGAATCGGGTGACGGAGATGAACAGCCCCAAGCCCACACCACCGCTGGCctgc                                                                                                                                                >  1:342128/1‑149 (MQ=255)        tCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGc                                                                                                                                        >  2:70721/1‑149 (MQ=255)         cTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGc                                                                                                                                        <  7:101058/148‑1 (MQ=255)          tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGCCCTCCGCTGGCCTGCTGGGCAGCGa                                                                                                                                      >  8:59360/1‑149 (MQ=255)          tttGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGa                                                                                                                                      >  3:313179/1‑149 (MQ=255)            tGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                    >  2:226411/1‑146 (MQ=255)            tGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCgg                                                                                                                                      <  1:70721/147‑2 (MQ=255)             gAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                    <  1:226411/148‑4 (MQ=255)             gAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGCCCAGCAGGCCAGCggag                                                                                                                                    <  5:139839/148‑4 (MQ=255)                 gTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                               >  3:257629/1‑149 (MQ=255)                     ccAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                     <  1:319608/139‑1 (MQ=255)                     ccAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGAc                                                                                                                                     >  2:319608/1‑139 (MQ=255)                      cAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTggg                                                                                                                                            <  8:71152/131‑1 (MQ=255)                       aGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCAACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActg                                                                                                                         >  4:374093/1‑149 (MQ=255)                        gCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTCCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgc                                                                                                                        >  6:47228/1‑149 (MQ=255)                         cgtattttGCTGGGTTATTGTCCCACTCTTGCCCTCCTCTCTATCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgct                                                                                                                       <  1:136502/143‑1 (MQ=255)                          ggCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAActgctg                                                                                                                      <  4:81970/149‑1 (MQ=255)                                  cAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                              >  5:108114/1‑149 (MQ=255)                                  cAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                              >  6:363389/1‑149 (MQ=255)                                                cACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTggcg                                                                                                <  3:374093/149‑1 (MQ=255)                                                    cTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGgcgc                                                                                            <  2:282314/149‑1 (MQ=255)                                                       gCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCt                                                                                           <  3:154401/147‑1 (MQ=255)                                                        cccTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgat                                                                                        <  7:119838/149‑1 (MQ=255)                                                         cctcctCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgata                                                                                       >  7:320342/1‑149 (MQ=255)                                                          ctcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgatac                                                                                      <  4:257629/149‑1 (MQ=255)                                                           tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGAAAGATCCGCTGCTGGTTGGCGGCGCTgatag                                                                                      >  5:267058/1‑147 (MQ=255)                                                           tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGAAAGATCCGCTGCTGGTTGGCGGCGCTga                                                                                         <  6:267058/145‑1 (MQ=255)                                                           tcctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCAAAGCCCAGACCTCCGCTGGCCTCCTGGGCACCGACCCGCACGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCCATgatacc                                                                                     >  6:35229/1‑149 (MQ=255)                                                             ctcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccga                                                                                   <  6:108114/149‑1 (MQ=255)                                                              tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGa                                                                                                                    <  6:289874/115‑1 (MQ=255)                                                              tcGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGa                                                                                                                    >  5:289874/1‑115 (MQ=255)                                                                gCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTgataccgatac                                                                                >  7:128945/1‑149 (MQ=255)                                                                  tttCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                               <  6:210508/100‑1 (MQ=255)                                                                  tttCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTc                                                                                                                               >  5:210508/1‑100 (MQ=255)                                                                  tttCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGataccgataccg                                                                              >  3:123921/1‑149 (MQ=255)                                                                    tCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGataccgataccg                                                                              <  5:363389/147‑1 (MQ=255)                                                                     cAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTg                                                                            <  7:214021/148‑1 (MQ=255)                                                                      aaCCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGAt                                                                          >  3:357906/1‑149 (MQ=255)                                                                         cTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcct                                                                       >  2:370921/1‑149 (MQ=255)                                                                            tATCCGGCCACCGTGCGCCTGGACG‑ATGTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctgca                                                                    <  6:40439/149‑1 (MQ=255)                                                                                       cGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                         <  1:370921/149‑1 (MQ=255)                                                                                       cGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                         <  7:59360/149‑1 (MQ=255)                                                                                       cGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAg                                                         <  4:123921/149‑1 (MQ=255)                                                                                              cTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgc                                                  <  4:313179/149‑1 (MQ=255)                                                                                                  aCG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                              >  3:34917/1‑85 (MQ=255)                                                                                                   actaTCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATcc                                                                                                              <  4:34917/82‑1 (MQ=255)                                                                                                    g‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCCTGCCGCCCt                                            >  8:218831/1‑149 (MQ=255)                                                                                                    g‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCt                                            >  7:16126/1‑149 (MQ=255)                                                                                                           cT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCAt                                                    >  5:155535/1‑135 (MQ=255)                                                                                                           cT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCAt                                                    <  6:155535/135‑1 (MQ=255)                                                                                                           cT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTg                                           <  8:128945/144‑1 (MQ=255)                                                                                                           cT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCt                                         <  2:342128/146‑1 (MQ=255)                                                                                                            t‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGaa                                     >  5:157038/1‑149 (MQ=255)                                                                                                                    aTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                              >  2:336700/1‑134 (MQ=255)                                                                                                                    aTGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATgccgcc                                              <  1:336700/134‑1 (MQ=255)                                                                                                                        aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                         >  1:180139/1‑103 (MQ=255)                                                                                                                        aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATc                                                                         <  2:180139/103‑1 (MQ=255)                                                                                                                        aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                        >  3:241093/1‑136 (MQ=255)                                                                                                                        aCAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTc                                        <  4:241093/136‑1 (MQ=255)                                                                                                                         cAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCt                          >  4:359228/1‑149 (MQ=255)                                                                                                                            ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                      >  2:139829/1‑102 (MQ=255)                                                                                                                            ccccAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATcctg                                                                      <  1:139829/102‑1 (MQ=255)                                                                                                                                   cccAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGAtt                <  5:47228/149‑1 (MQ=255)                                                                                                                                        aCCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCt           >  5:358401/1‑149 (MQ=255)                                                                                                                                         ccTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAATATGCTGTGATGTTTGGCGGCGCTGATACAGATACAGTGATGCTGGTACTGCTACCAGGCCTTGCAGCCCTTCTCGAACTCCTGCACGTGCACCGGATTGctctt          <  8:320342/149‑2 (MQ=255)                                                                                                                                                  gCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAAc  <  8:16126/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTTTTGTTCTTTGAATGTTGCCAGCGGCAATTGCAGGGTGAATGTGGCACCCTTGCCCTCCTCGCTTTCAACCTGTATCCGGCCACCGTGCGCCTGGACG‑ATCTGCT‑‑CGGAGATGTACAGCCCCAAGCCCAGACCTCCGCTGGCCTGCTGGGCAGCGACCCGCTCGAACTGCTGGAAGATCCGCTGCTGGTTGGCGGCGCTGATACCGATACCGTGATCCTGCACCTGCACCCAGGCCATGCCGCCCTGCTCGAACACCCGCACCTGCACCGGATTGCGCTCGCCGTAAC  >  NC_002947/4286536‑4286825

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