Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	298,117	C→G	6.3%	intergenic (‑114/+40)	gshA ← / ← PP_0244	glutamate‑‑cysteine ligase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	298,117	0	C	G	6.3%	203.5 / 4.4	63	intergenic (‑114/+40)	gshA/PP_0244	glutamate‑‑cysteine ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

TCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297976‑298265                                                                                                                                              |                                                                                                                                                     tCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCAc                                                                                                                                               >  4:202477/1‑149 (MQ=255)              gaggagTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                   <  2:88800/148‑5 (MQ=255)                ggagTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCg                                                                                                                                 >  5:278256/1‑148 (MQ=255)                   gTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             <  3:202477/149‑1 (MQ=255)                    ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  7:4689/1‑148 (MQ=255)                    ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  7:39392/1‑148 (MQ=255)                    ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  3:41414/1‑148 (MQ=255)                       gCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGTTACt                                                                                                                         >  8:302144/1‑149 (MQ=255)                         tCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGc                                                                                                                       >  8:218358/1‑149 (MQ=255)                           aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                   <  3:149304/135‑5 (MQ=255)                           aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                   >  4:149304/1‑131 (MQ=255)                           aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCtt                                                                                                                     <  7:18770/149‑1 (MQ=255)                            aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  5:132230/149‑1 (MQ=255)                            aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  2:204152/149‑1 (MQ=255)                             gATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTgg                                                                                                                   >  1:239018/1‑149 (MQ=255)                                gTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGGGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                  >  5:311568/1‑147 (MQ=255)                                  aTTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGc                                                                                                              >  2:184264/1‑149 (MQ=255)                                     ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            <  5:270069/148‑1 (MQ=255)                                     ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            <  4:51383/148‑1 (MQ=255)                                     ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            <  6:278256/148‑1 (MQ=255)                                      cTTGACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTagt                                                                                                          <  4:241091/149‑2 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  3:326528/1‑129 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             <  4:326528/129‑1 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  6:189718/138‑1 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    >  5:189718/1‑138 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                  >  5:336612/1‑140 (MQ=255)                                       ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                  <  6:336612/140‑1 (MQ=255)                                                 aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGaa                                                                                                <  7:218358/148‑1 (MQ=255)                                                  gTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAgg                                                                                              >  5:230180/1‑149 (MQ=255)                                                      ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTcc                                                                                          >  2:175664/1‑149 (MQ=255)                                                       cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTg                                                                                                                 <  2:127839/125‑1 (MQ=255)                                                       cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTg                                                                                                                 >  1:127839/1‑125 (MQ=255)                                                       cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCt                                                                                         >  7:99034/1‑149 (MQ=255)                                                        ccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        <  8:99034/149‑1 (MQ=255)                                                          aaTATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTgc                                                                                      <  4:41414/149‑1 (MQ=255)                                                               gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        >  6:7193/1‑142 (MQ=255)                                                               gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        <  5:7193/142‑1 (MQ=255)                                                                ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                <  8:4689/149‑1 (MQ=255)                                                                ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                <  7:144400/149‑1 (MQ=255)                                                                ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                <  4:85079/149‑1 (MQ=255)                                                                      gagaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccac                                                                           <  1:184264/148‑1 (MQ=255)                                                                         aTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGc                                                                       >  1:102035/1‑149 (MQ=255)                                                                            aCGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTg                                                                    <  6:311568/149‑1 (MQ=255)                                                                                    cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             <  8:39392/148‑1 (MQ=255)                                                                                    cAAGCGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             >  3:50905/1‑148 (MQ=255)                                                                                       gggTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAAc                                                         <  3:123324/149‑1 (MQ=255)                                                                                               agaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCgg                                                 >  3:212751/1‑149 (MQ=255)                                                                                                   aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                             >  3:61468/1‑149 (MQ=255)                                                                                                   aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCACTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                             >  1:302521/1‑149 (MQ=255)                                                                                                       aCTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGac                                          >  5:60415/1‑148 (MQ=255)                                                                                                           gCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                       >  5:260007/1‑147 (MQ=255)                                                                                                             gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   <  2:67905/133‑1 (MQ=255)                                                                                                             gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   >  1:67905/1‑133 (MQ=255)                                                                                                              ttGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGTTACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAggg                                  <  7:302144/149‑1 (MQ=255)                                                                                                                gTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGt                                 <  3:246544/148‑1 (MQ=255)                                                                                                                   gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTgcgc                             <  2:239018/149‑1 (MQ=255)                                                                                                                    ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTgg                                                                               >  2:131676/1‑98 (MQ=255)                                                                                                                    ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTgg                                                                               <  1:131676/98‑1 (MQ=255)                                                                                                                     cTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccg                           >  2:237228/1‑149 (MQ=255)                                                                                                                      tGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCGGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccgc                          >  2:50951/1‑149 (MQ=255)                                                                                                                               cgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgc                 <  4:212751/149‑1 (MQ=255)                                                                                                                                ggTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTg                  >  4:90566/1‑147 (MQ=255)                                                                                                                                 ttttcagcagTGCGCTGGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcag               >  8:338472/5‑149 (MQ=255)                                                                                                                                              cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg  <  2:102035/149‑1 (MQ=255)                                                                                                                                              cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg  >  5:366900/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     TCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297976‑298265

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