Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	462,592	A→T	10.0%	intergenic (‑128/+13)	pqqF ← / ← davA	coenzyme PQQ synthesis protein F/5‑aminopentanamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	462,592	0	A	T	10.0%	190.9 / 8.5	67	intergenic (‑128/+13)	pqqF/davA	coenzyme PQQ synthesis protein F/5‑aminopentanamidase
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base T (3/4);  total (31/36)
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.80e-01

GCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAG  >  NC_002947/462447‑462733                                                                                                                                                  |                                                                                                                                              gCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGAGTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAg                                                                                                                                            <  8:145303/149‑1 (MQ=255)           tCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCTATGCGGcccattc                                                                                                                                  <  4:6640/149‑5 (MQ=255)           tCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCAATGCGGcccattc                                                                                                                                  <  5:236540/149‑5 (MQ=255)             aGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCTATGCGGcccattcgc                                                                                                                                <  8:149818/149‑7 (MQ=255)               gCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGc                                                                                                                              <  6:288685/149‑1 (MQ=255)               gCATGGGGGCTCCGGTGGCTGCTCTCCAAGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCTATGCGGcccattc                                                                                                                                  <  6:185997/145‑5 (MQ=255)               gCATGGGGGCTCCGGTGGCTGCTCTCCAAGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCTATGCGGcccattc                                                                                                                                  >  5:185997/1‑141 (MQ=255)                         tCCGGTGGCTGCTCTGCATGCTAACCTATTTGTTTGCAGGCCCGGCCTTTTCGGGGGTATGGCGCCGTTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAg                                                                                                                    <  7:203269/149‑1 (MQ=255)                           cGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAg                                                                                                                    <  6:275313/147‑1 (MQ=255)                                  tgctCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                           >  4:243002/1‑149 (MQ=255)                                  tgctCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                           >  4:67184/1‑149 (MQ=255)                                    ctctCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                           <  6:139080/147‑1 (MQ=255)                                         cATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGgc                                                                                                      <  7:41544/147‑1 (MQ=255)                                         cATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTAGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGgcg                                                                                                     >  7:202111/1‑148 (MQ=255)                                          aTGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGt                                                                                                                                  >  3:339007/1‑118 (MQ=255)                                          aTGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGt                                                                                                                                  <  4:339007/118‑1 (MQ=255)                                              tAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATc                                                                                               >  8:143833/1‑149 (MQ=255)                                               aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATACCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                              >  5:25549/1‑149 (MQ=255)                                                aCCGATTTGTTTGCAGGCTCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                              <  3:350918/148‑1 (MQ=255)                                                 ccGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                            <  2:211752/149‑1 (MQ=255)                                                  cGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                            <  1:128083/148‑1 (MQ=255)                                                   gATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCTGTCa                                                                                          <  1:286707/149‑1 (MQ=255)                                                   gATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                            >  7:345564/1‑147 (MQ=255)                                                   gATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                            >  2:128083/1‑147 (MQ=255)                                                     tttgtttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCttt                                                                                                                          >  4:281166/1‑115 (MQ=255)                                                     tttgtttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAgg                                                                                        <  7:143833/149‑1 (MQ=255)                                                         tttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAgg                                                                                    <  1:147074/149‑1 (MQ=255)                                                         tttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAgg                                                                                    <  3:92507/149‑1 (MQ=255)                                                          ttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAg                                                                                     <  4:13214/147‑1 (MQ=255)                                                          ttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAg                                                                                     >  3:13214/1‑147 (MQ=255)                                                           tgCAGGCCCGGCCTTTTCGCGGGTATGGCGTCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAggg                                                                                   >  5:56524/1‑148 (MQ=255)                                                            gCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGa                                                                                  >  3:295831/1‑148 (MQ=255)                                                             cAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGaaa                                                                                <  4:295831/149‑1 (MQ=255)                                                                gCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGc                                                                             >  5:149131/1‑149 (MQ=255)                                                                 cccGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGcc                                                                            >  8:11789/1‑149 (MQ=255)                                                                  ccGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGccc                                                                           >  4:254923/1‑149 (MQ=255)                                                                    ggCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCACGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTg                                                                         >  3:136292/1‑149 (MQ=255)                                                                       cTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTgcgc                                                                      >  4:166763/1‑149 (MQ=255)                                                                        ttttCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAggg                                                                                   >  7:333707/1‑135 (MQ=255)                                                                        ttttCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAggg                                                                                   <  8:333707/135‑1 (MQ=255)                                                                               gggTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGcac                                                              >  3:321369/1‑149 (MQ=255)                                                                                 gTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGcacca                                                            >  8:294015/1‑149 (MQ=255)                                                                                  tATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGTAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGca                                                               >  4:318026/1‑145 (MQ=255)                                                                                  tATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGTAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGca                                                               <  3:318026/145‑1 (MQ=255)                                                                                          cGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCAt                                                   >  6:8088/1‑149 (MQ=255)                                                                                          cGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCAt                                                   <  3:166763/149‑1 (MQ=255)                                                                                               tGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTc                                               <  4:321369/148‑1 (MQ=255)                                                                                                 aaaGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAg                                            <  8:202111/149‑1 (MQ=255)                                                                                                    ggCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGctc                                         >  7:733/1‑149 (MQ=255)                                                                                                     gCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGctct                                        <  7:11789/149‑1 (MQ=255)                                                                                                     gCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGctct                                        <  4:136292/149‑1 (MQ=255)                                                                                                     gCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGctct                                        >  2:304733/1‑149 (MQ=255)                                                                                                      cAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTg                                       <  5:8088/149‑1 (MQ=255)                                                                                                       aGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGc                                      <  6:56524/149‑1 (MQ=255)                                                                                                               tCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACt                              >  2:343368/1‑149 (MQ=255)                                                                                                                 ccGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTgg                            <  6:149131/149‑1 (MQ=255)                                                                                                                 ccGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCAGCTCATGCTCGAGCTCTGCCAACAACTgg                            >  5:110948/1‑149 (MQ=255)                                                                                                                          cGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATcc                   <  3:254923/149‑1 (MQ=255)                                                                                                                                   ccGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTTGCACTCATCCCGCCCGGcc          <  6:25549/149‑1 (MQ=255)                                                                                                                                   ccGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATcc                   <  2:188849/140‑1 (MQ=255)                                                                                                                                   ccGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATcc                   >  1:188849/1‑140 (MQ=255)                                                                                                                                     tgacaggggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCa                          >  7:341078/7‑131 (MQ=255)                                                                                                                                     tgacaggggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCa                          <  8:341078/125‑1 (MQ=255)                                                                                                                                          tGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCa   <  8:345564/149‑1 (MQ=255)                                                                                                                                          tGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCa   <  8:733/149‑1 (MQ=255)                                                                                                                                           gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg  >  1:359816/1‑149 (MQ=255)                                                                                                                                            ggCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg  >  3:290661/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                              GCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTAGGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAG  >  NC_002947/462447‑462733

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