Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,753,542	C→A	15.6%	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,542	0	C	A	15.6%	165.3 / 17.8	63	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/25);  new base A (3/7);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CGACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTGA  >  NC_002947/5753393‑5753690                                                                                                                                                         |                                                                                                                                                     cgacgaCTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGTACTGTGATACAG‑‑CTGCGGAGCTGGTTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGcctccctc                                                                                                                                                      >  1:625093/1‑150 (MQ=255)  gacgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCg                                                                                                                                                     <  2:948713/150‑1 (MQ=255)   acgacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCgg                                                                                                                                                    <  2:436136/150‑1 (MQ=255)     gacTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGa                                                                                                                                                  <  1:277174/150‑1 (MQ=255)        ttCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGc                                                                                                                                               >  2:729592/1‑150 (MQ=255)          cATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCg                                                                                                                                             >  2:623120/1‑150 (MQ=255)           aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGTGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGt                                                                                                                                            <  1:313707/150‑1 (MQ=255)           aTCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCCAAACCGCCTCCCTCGGGAGGCCGt                                                                                                                                            >  2:551887/1‑150 (MQ=255)             cGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                          <  2:207871/150‑1 (MQ=255)              gATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGtttt                                                                                                                                         <  1:1041820/150‑1 (MQ=255)                   cTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGttt                                                                                                                                    <  1:167583/150‑1 (MQ=255)                                      gaagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAgc                                                                                                                 <  1:835882/150‑1 (MQ=255)                                       aagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCt                                                                                                                                  >  2:534372/1‑132 (MQ=255)                                       aagaagaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCt                                                                                                                                  <  1:534372/132‑1 (MQ=255)                                           agaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGcgccgc                                                                                                            <  1:106152/150‑1 (MQ=255)                                            gaaATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCt                                                                                                           >  2:882967/1‑150 (MQ=255)                                              aaTCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                          <  1:823439/117‑1 (MQ=255)                                              aaTCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                          >  2:823439/1‑117 (MQ=255)                                               aTCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGc                                                                                                                       <  1:178414/135‑1 (MQ=255)                                               aTCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGc                                                                                                                       >  2:178414/1‑135 (MQ=255)                                                   aGGTCCGTGCCTTCCTACACCCGCTGGAATACGAGCTGTAATACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGCCATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGc                                                                                                    <  1:551887/150‑1 (MQ=255)                                                   aGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGc                                                                                                    <  1:469611/150‑1 (MQ=255)                                                   aGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGc                                                                                                    >  1:614633/1‑150 (MQ=255)                                                   aGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGc                                                                                                    <  1:882967/150‑1 (MQ=255)                                                   aGGTCCGTACCTTCGTACACCCGCTCGAATACAAGCTGTACTACAG‑‑CTGCTAATCTGATTGGCGCCTCGCGCCC‑ATGACATGATCAAAACAGCCTCCCTCGGGAGGGCGTTTTAGTCTCTACCCACCTTGCAGAATCGCCACAGATTCGc                                                                                                    >  2:469611/1‑150 (MQ=255)                                                   aGGTCCGCACCTTCGTACGCCCGCTGGAATACGAGCTGTACTACAGGCCTGCTTAGCGG‑‑T‑CATCCTAAATTCGTTTTATATGATCAAAGCCGTCACCCCAGGGAGGCTGTTTTTGGGTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGc                                                                                                    <  1:275534/150‑1 (MQ=255)                                                            ccTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGc                                                                                           >  1:780628/1‑150 (MQ=255)                                                                  tACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGaa                                                                                                                   <  1:554972/120‑1 (MQ=255)                                                                  tACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGaa                                                                                                                   >  2:554972/1‑120 (MQ=255)                                                                   acacCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                          >  1:861449/1‑96 (MQ=255)                                                                   acacCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCGAGGGAGGCCGttt                                                                                                                                          <  2:861449/96‑1 (MQ=255)                                                                     acCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATAACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTa                                                                                  >  1:1180796/1‑150 (MQ=255)                                                                      cccGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGAACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTAt                                                                                 >  1:859724/1‑150 (MQ=255)                                                                       ccGCTGGAATACGAGCTGTACTACAG‑‑CCGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATc                                                                                >  1:293099/1‑150 (MQ=255)                                                                                   gAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCt                                                                               >  1:323068/1‑139 (MQ=255)                                                                                   gAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTaa                                                                    >  1:365545/1‑150 (MQ=255)                                                                                   gAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGATTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCt                                                                               <  2:323068/139‑1 (MQ=255)                                                                                     gCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAtt                                                                  <  1:623120/150‑1 (MQ=255)                                                                                     gCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAAtt                                                                  >  1:943571/1‑150 (MQ=255)                                                                                        gTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTcac                                                               >  1:364956/1‑150 (MQ=255)                                                                                          actacAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCacac                                                             <  1:628953/150‑1 (MQ=255)                                                                                                   ctgctgATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGccc                                                      >  1:148965/1‑150 (MQ=255)                                                                                                         atctgatTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAAtt                                                >  1:181811/1‑150 (MQ=255)                                                                                                          tctgatTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGCCACGCTATGCTCTATCTTGGTGACTTAACACACACCAGGCCCTCAATTg                                               <  1:933715/150‑1 (MQ=255)                                                                                                           ctgatTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGc                                              >  2:184221/1‑150 (MQ=255)                                                                                                             gatTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTAACTGATGAATTCACACCATGCCCTCAATTGCTc                                            <  1:30962/150‑1 (MQ=255)                                                                                                              atTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTcc                                           >  1:29470/1‑150 (MQ=255)                                                                                                                   cgcCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAAt                                      <  2:365545/150‑1 (MQ=255)                                                                                                                    gcCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATc                                     <  2:293099/150‑1 (MQ=255)                                                                                                                         gcgcCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAg                                                                                                                                                 <  1:373789/37‑1 (MQ=255)                                                                                                                         gcgcCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAg                                                                                                                                                 >  2:373789/1‑37 (MQ=255)                                                                                                                         gcgcCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCACTGATTCGCCTGTCACGCTATg                                                                                       <  2:924575/95‑1 (MQ=255)                                                                                                                         gcgcCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCACTGATTCGCCTGTCACGCTATg                                                                                       >  1:924575/1‑95 (MQ=255)                                                                                                                                aTGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGtt                          <  2:1193932/150‑1 (MQ=255)                                                                                                                                aTGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGtt                          >  2:537774/1‑150 (MQ=255)                                                                                                                                  gACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTaa                        >  2:396941/1‑150 (MQ=255)                                                                                                                                      tGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGt                    <  1:396941/150‑1 (MQ=255)                                                                                                                                            aaaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCAcc              >  2:245632/1‑150 (MQ=255)                                                                                                                                             aaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCa             >  1:268742/1‑150 (MQ=255)                                                                                                                                             aaCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCa             >  1:14034/1‑150 (MQ=255)                                                                                                                                              aCGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCaa            <  2:181811/150‑1 (MQ=255)                                                                                                                                              aCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAAc                      <  1:173487/140‑1 (MQ=255)                                                                                                                                              aCGGCCTCCCGAGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAAc                      >  2:173487/1‑140 (MQ=255)                                                                                                                                                ggCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATGCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGa          <  2:780628/150‑1 (MQ=255)                                                                                                                                                     ccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc     <  2:148965/150‑1 (MQ=255)                                                                                                                                                     ccctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTc     <  2:943571/150‑1 (MQ=255)                                                                                                                                                      cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAg                                                                                                                     >  2:298968/1‑37 (MQ=255)                                                                                                                                                      cctcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAg                                                                                                                     <  1:298968/37‑1 (MQ=255)                                                                                                                                                        tcGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTGa  <  1:184221/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                     CGACGACTTCATCGATGCCTTCATCGAGCTGAAGAGCGAAGAAGAAATCAAGGTCCGTACCTTCGTACACCCGCTGGAATACGAGCTGTACTACAG‑‑CTGCTGATCTGATTGGCGCCTCGCGCCG‑ATGACATGATCAAAACGGCCTCCCTCGGGAGGCCGTTTTTGTTTCTGCGCACCTTGCAGAAGCGCCGCTGATTCGCCTGTCACGCTATGCTCTATCTTGGTGCATTAATTCACACCATGCCCTCAATTGCTCCCCAATCTGGTTCAGGTTAAACGTTTCACCAAGACTTTCTGA  >  NC_002947/5753393‑5753690

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