Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,753,486	C→A	10.5%	intergenic (+47/‑319)	glnA → / → glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,753,486	0	C	A	10.5%	191.7 / 8.2	65	intergenic (+47/‑319)	glnA/glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (3/4);  total (31/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 = 8.58e-01

TTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTC  >  NC_002947/5753344‑5753631                                                                                                                                               |                                                                                                                                                  ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCACGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAgg                                                                                                                                             >  2:316632/1‑149 (MQ=255)   cATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAgg                                                                                                                                             <  2:260598/147‑1 (MQ=255)   cATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGcc                                                                                                                                           <  4:71307/149‑1 (MQ=255)   cATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGcc                                                                                                                                           >  5:88529/1‑149 (MQ=255)   cATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                            <  3:358925/148‑1 (MQ=255)       gATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                       <  5:295974/149‑1 (MQ=255)         ggCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCTTACCTTCGTACACCCGCGGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttttt                                                                                                                                     <  7:195285/148‑1 (MQ=255)            cTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGtt                                                                                                                                  >  1:226515/1‑149 (MQ=255)            cTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGtt                                                                                                                                  >  6:79627/1‑149 (MQ=255)             ttCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTg                                                                                                                                    <  1:300969/146‑1 (MQ=255)             ttCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTg                                                                                                                                    >  2:300969/1‑146 (MQ=255)               cATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGttt                                                                                                                                 <  4:17228/147‑1 (MQ=255)               cATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttttgatcat                                                                                                                                <  6:282395/148‑7 (MQ=255)                 tCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                             >  7:292224/1‑149 (MQ=255)                 tCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                             <  7:31135/149‑1 (MQ=255)                 tCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTgc                                                                                                                             <  6:269527/149‑1 (MQ=255)                   gAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTgcgc                                                                                                                           <  4:363125/149‑1 (MQ=255)                      cTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCAcc                                                                                                                        <  6:168324/149‑1 (MQ=255)                        gAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCtt                                                                                                                      >  6:166178/1‑149 (MQ=255)                                   aagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgcc                                                                                                           >  4:46285/1‑149 (MQ=255)                                    agaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccg                                                                                                          >  6:414142/1‑149 (MQ=255)                                              gTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGcc                                                                                                <  6:88529/149‑1 (MQ=255)                                                 cGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGtt                                                                                                                                        >  2:241935/1‑106 (MQ=255)                                                 cGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGtt                                                                                                                                        <  1:241935/106‑1 (MQ=255)                                                   tACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCa                                                                                           >  6:310705/1‑149 (MQ=255)                                                          gTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATg                                                                                    >  4:119818/1‑149 (MQ=255)                                                           tACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGc                                                                                   <  5:314848/149‑1 (MQ=255)                                                              acCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                 >  6:175243/1‑148 (MQ=255)                                                              acCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                 >  5:281122/1‑148 (MQ=255)                                                               cccGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTa                                                                               <  1:316632/149‑1 (MQ=255)                                                                ccGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTa                                                                               >  2:83117/1‑148 (MQ=255)                                                                      gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                          <  6:281122/147‑1 (MQ=255)                                                                      gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                          <  1:83117/147‑1 (MQ=255)                                                                      gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGt                                                                        <  3:46285/149‑1 (MQ=255)                                                                          agcaGCTGTAGTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCAt                                                                    >  8:377783/4‑149 (MQ=255)                                                                            gAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCAtt                                                                   >  5:274853/1‑148 (MQ=255)                                                                              gCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAt                                                                <  5:79627/149‑1 (MQ=255)                                                                              gCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAt                                                                <  2:226515/149‑1 (MQ=255)                                                                                  tactacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTca                                                             <  5:157316/148‑1 (MQ=255)                                                                                     tacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACg                                                                                         >  4:47535/1‑117 (MQ=255)                                                                                     tacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACg                                                                                         <  3:47535/117‑1 (MQ=255)                                                                                       cAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCa                                                        <  7:377783/148‑1 (MQ=255)                                                                                       cAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCAt                                                       >  8:22043/1‑149 (MQ=255)                                                                                                   tgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGc                                           >  1:349475/1‑149 (MQ=255)                                                                                                   tgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGc                                           >  3:211035/1‑149 (MQ=255)                                                                                                            ccTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCACCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATc                                  <  8:292224/149‑1 (MQ=255)                                                                                                              tCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                >  3:210029/1‑149 (MQ=255)                                                                                                                gcgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTgg                               <  2:349475/148‑1 (MQ=255)                                                                                                                gcgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                <  5:175243/147‑1 (MQ=255)                                                                                                                 cgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGtt                             <  5:414142/149‑1 (MQ=255)                                                                                                                 cgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGtt                             <  7:22043/149‑1 (MQ=255)                                                                                                                          cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTaa                     >  8:322715/1‑148 (MQ=255)                                                                                                                           aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGccc                                                   >  4:113320/1‑117 (MQ=255)                                                                                                                           aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGccc                                                   <  3:113320/117‑1 (MQ=255)                                                                                                                           aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                <  4:38618/136‑1 (MQ=255)                                                                                                                           aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTg                                >  3:38618/1‑136 (MQ=255)                                                                                                                              aTCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCTGCTGCTTCGCCTGTCACGGTATGCTCTATCTTGGTGCATTATTTCACAACATGCCCTCAATTGCTACCCAATCTGGTTAAGGTTAAACGtt                >  7:106867/1‑149 (MQ=255)                                                                                                                                aaaaaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGttt               >  8:334973/2‑148 (MQ=255)                                                                                                                                aaaaaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTc              <  7:334973/148‑1 (MQ=255)                                                                                                                                aaaaaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGATCCCCAATCTGGTTCAGGTTAAACGttt               >  2:378418/2‑148 (MQ=255)                                                                                                                                 aaaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCa             <  3:119818/149‑1 (MQ=255)                                                                                                                                  aaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAc            >  6:339104/1‑149 (MQ=255)                                                                                                                                     cGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTcc                                        <  7:223409/118‑1 (MQ=255)                                                                                                                                     cGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTcc                                        >  8:223409/1‑118 (MQ=255)                                                                                                                                      ggCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAACCGTTTCAca           >  8:322237/1‑145 (MQ=255)                                                                                                                                            cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc  <  5:310705/149‑1 (MQ=255)                                                                                                                                            cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc  <  3:33732/149‑1 (MQ=255)                                                                                                                                            cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc  >  4:33732/1‑149 (MQ=255)                                                                                                                                            cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc  <  4:211035/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  TTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTC  >  NC_002947/5753344‑5753631

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