Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,753,486	C→A	7.2%	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	7.2%	206.8 / 6.0	70	intergenic (+47/‑319)	glnA/glnL	glutamine synthetase/two‑component system sensor histidine kinase/phosphatase GlnL/GlnG
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base A (3/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 = 9.13e-02

CGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCT  >  NC_002947/5753340‑5753632                                                                                                                                                   |                                                                                                                                                   cgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                  <  3:38160/149‑1 (MQ=255)     ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAgg                                                                                                                                              >  3:172102/1‑149 (MQ=255)      tCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                             >  2:312495/1‑149 (MQ=255)           gATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttt                                                                                                                                        >  8:346330/1‑149 (MQ=255)           gATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGtt                                                                                                                                         <  2:139776/148‑1 (MQ=255)           gATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtt                                                                                                                                         <  4:318314/148‑1 (MQ=255)             tGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGttttt                                                                                                                                      <  2:157634/149‑1 (MQ=255)                cTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGtt                                                                                                                                   >  7:321040/1‑149 (MQ=255)                    aTCGAGCTGAAGAGCGATGAAGAAATCACGGTCCGTCCCTTCGTACACCGGCTTGATTTCGAGCTGTTCTTGAGCTGCTGATCTGATTGGCTCCTCGCGCCGATGACATGATCATAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCt                                                                                                                                <  4:82109/148‑1 (MQ=255)                                  cGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTCGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAg                                                                                                                   >  2:244135/1‑147 (MQ=255)                                    aagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAgc                                                                                                               <  6:374425/149‑1 (MQ=255)                                    aagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtttt                                                                                                                                       >  5:99842/1‑125 (MQ=255)                                    aagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtttt                                                                                                                                       <  6:99842/125‑1 (MQ=255)                                       aagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAgcgc                                                                                                             >  3:331639/1‑148 (MQ=255)                                       aagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAgcgc                                                                                                             <  4:331639/148‑1 (MQ=255)                                        agaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgcc                                                                                                            <  3:46441/148‑1 (MQ=255)                                         gaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccgc                                                                                                          >  7:17219/1‑149 (MQ=255)                                              cAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGCCATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGAtt                                                                                                     >  7:238002/1‑149 (MQ=255)                                                      gTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTc                                                                                             <  4:172102/149‑1 (MQ=255)                                                       tACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTc                                                                                             <  8:212033/148‑1 (MQ=255)                                                       tACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTc                                                                                             >  7:212033/1‑148 (MQ=255)                                                              gTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTTCAGAATCGCCACTGATTCGCCTGTCACGCTATg                                                                                     <  7:57622/149‑1 (MQ=255)                                                                  acCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                  >  7:246961/1‑148 (MQ=255)                                                                  acCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                  >  7:99209/1‑148 (MQ=255)                                                                   cccGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTa                                                                                >  1:364848/1‑149 (MQ=255)                                                                    ccGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccgc                                                                                                          >  8:207939/1‑122 (MQ=255)                                                                    ccGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccgc                                                                                                          <  7:207939/122‑1 (MQ=255)                                                                      gCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCt                                                                             >  4:250384/1‑149 (MQ=255)                                                                        tGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                           >  4:260182/1‑149 (MQ=255)                                                                          gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCTGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTg                                                                           <  7:317275/147‑1 (MQ=255)                                                                          gAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGt                                                                         <  8:99209/149‑1 (MQ=255)                                                                            aTACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGATCTATCTTGGTg                                                                        >  7:149785/1‑148 (MQ=255)                                                                              aCGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCAt                                                                     >  2:344279/1‑149 (MQ=255)                                                                                 aGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACTCTATGATCTATCGTCGTGCATTaa                                                                  >  6:262524/1‑149 (MQ=255)                                                                                 aGCTGGACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTaa                                                                  <  6:87225/149‑1 (MQ=255)                                                                                   cTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAtt                                                                <  8:246961/149‑1 (MQ=255)                                                                                     gTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTca                                                              >  7:334152/1‑149 (MQ=255)                                                                                        ctacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCacac                                                           <  8:321040/149‑1 (MQ=255)                                                                                        ctacAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCacac                                                           <  2:347236/149‑1 (MQ=255)                                                                                              ctgctgATCTGATTGGCGGCTCGCGGCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGcc                                                     <  7:250844/149‑1 (MQ=255)                                                                                                  tgatctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTc                                                  >  1:177128/1‑148 (MQ=255)                                                                                                   gatctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCGCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTg                                                                                                                      >  1:262903/1‑79 (MQ=255)                                                                                                   gatctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTg                                                                                                                      <  2:262903/79‑1 (MQ=255)                                                                                                   gatctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCaa                                                >  7:179685/1‑149 (MQ=255)                                                                                                    atctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCaa                                                <  8:238002/148‑1 (MQ=255)                                                                                                     tctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAtt                                                                <  3:325282/131‑1 (MQ=255)                                                                                                     tctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAtt                                                                >  4:325282/1‑131 (MQ=255)                                                                                                      ctgatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAAtt                                              <  2:177128/148‑1 (MQ=255)                                                                                                        gatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                  <  3:329129/110‑1 (MQ=255)                                                                                                        gatTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGctc                                                                                  >  4:329129/1‑110 (MQ=255)                                                                                                         atTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCt                                           <  8:334152/148‑1 (MQ=255)                                                                                                            ggCGCCTCGCGCCGATTACATGATCACACCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTTCATTAAATCACACCATGCCCTCAATTGcccc                                         <  8:17219/147‑1 (MQ=255)                                                                                                               gcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAAt                                    >  6:374454/1‑149 (MQ=255)                                                                                                               gcCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAAt                                    >  2:177935/1‑149 (MQ=255)                                                                                                                 cTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCt                                  >  3:346538/1‑149 (MQ=255)                                                                                                                    gcgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGt                               >  5:324313/1‑149 (MQ=255)                                                                                                                     cgcCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGtt                              <  1:344279/149‑1 (MQ=255)                                                                                                                              cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAAtt                                              <  1:351479/124‑1 (MQ=255)                                                                                                                              cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAAtt                                              >  2:351479/1‑124 (MQ=255)                                                                                                                              cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTaaa                     <  3:260182/149‑1 (MQ=255)                                                                                                                              cATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTaa                      <  3:250384/148‑1 (MQ=255)                                                                                                                               aTGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAAc                    <  3:198823/149‑1 (MQ=255)                                                                                                                                     aaaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCa              >  1:257277/1‑149 (MQ=255)                                                                                                                                     aaaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCa              >  6:195284/1‑149 (MQ=255)                                                                                                                                      aaaCGGCCTCCCTCGGGATGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAc             <  7:346330/149‑1 (MQ=255)                                                                                                                                      aaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAc             >  3:61086/1‑149 (MQ=255)                                                                                                                                      aaaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCACAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCAAGCCCTCAATTGCTCCCCAATCTGGTTCAGTTTAAACGTTTCAc             >  8:167589/1‑149 (MQ=255)                                                                                                                                       aaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTTATTCACCCCATGCCGTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAcc            >  8:290146/1‑149 (MQ=255)                                                                                                                                       aaCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCag             >  8:248910/1‑147 (MQ=255)                                                                                                                                         cgacctcccgaGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAcc            >  4:177333/4‑147 (MQ=255)                                                                                                                                         cGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCaa          <  5:121665/149‑1 (MQ=255)                                                                                                                                               ccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACtt     >  3:309631/1‑148 (MQ=255)                                                                                                                                                cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACttt    <  4:346538/148‑1 (MQ=255)                                                                                                                                                cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc   <  8:179685/149‑1 (MQ=255)                                                                                                                                                cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc   <  4:61086/149‑1 (MQ=255)                                                                                                                                                 ctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCt  >  2:270583/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAGCTGCTGATCTGATTGGCGCCTCGCGCCGATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCT  >  NC_002947/5753340‑5753632

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