Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,342,601	G→A	11.1%	intergenic (+47/‑160)	panC → / → pgi‑II	pantothenate synthetase/glucose 6‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,342,601	0	G	A	11.1%	165.0 / 9.9	63	intergenic (+47/‑160)	panC/pgi‑II	pantothenate synthetase/glucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base A (4/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

ATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAG  >  NC_002947/5342456‑5342747                                                                                                                                                  |                                                                                                                                                   aTCAATGATCGAGACCTGGTGGTGATTGCGGCGGC‑TACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCAc                                                                                                                                                 <  3:156959/148‑1 (MQ=255)         tCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCt                                                                                                                                         >  5:51961/1‑149 (MQ=255)           gagaCCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTg                                                                                                                                        <  8:136057/148‑1 (MQ=255)            agaCCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTa                                                                                                                                      <  1:50723/149‑1 (MQ=255)                cTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAgggg                                                                                                                                  >  5:275226/1‑149 (MQ=255)                cTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAgggg                                                                                                                                  >  1:111469/1‑149 (MQ=255)                cTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAgggg                                                                                                                                  >  1:199865/1‑149 (MQ=255)                      gtgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgct                                                                                                                            >  3:183720/1‑149 (MQ=255)                       tgATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgctt                                                                                                                           <  8:67244/149‑1 (MQ=255)                            gcggcgGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGgcttgctt                                                                                                                       <  1:178045/148‑1 (MQ=255)                             cggcggCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGcttgcttg                                                                                                                      <  4:183720/148‑1 (MQ=255)                                 ggcTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATg                                                                                                                 >  1:38814/1‑149 (MQ=255)                                    ttACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcggg                                                                                                                              >  4:30268/1‑133 (MQ=38)                                    ttACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcggg                                                                                                                              <  3:30268/133‑1 (MQ=38)                                    ttACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTAGCGGGCAAGCCCGCTCCTGCACCTTTAACTGTAGGGGCGGGCTTGCTTGCAAAGccc                                                                                                              >  7:100014/1‑149 (MQ=255)                                    ttACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTGGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGccc                                                                                                              >  1:56669/1‑149 (MQ=255)                                      aCCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGcccc                                                                                                             <  2:56669/148‑1 (MQ=255)                                      aCCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGccc                                                                                                              <  1:245914/147‑1 (MQ=255)                                      aCCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGccc                                                                                                              <  1:242298/147‑1 (MQ=255)                                            gTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAGGTAAACGTGCAGgagcg                                                                                                                                <  6:134196/123‑1 (MQ=37)                                            gTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTACAGGTAAACGTGCAGgagcg                                                                                                                                >  5:134196/1‑123 (MQ=37)                                             tAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATcc                                                                                                     <  2:38814/149‑1 (MQ=255)                                             tAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcgg                                                                                                                               <  1:9459/123‑1 (MQ=37)                                             tAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCCCCTACAGGTAAACGTGCAGgagcgg                                                                                                                               >  2:9459/1‑123 (MQ=37)                                                 acGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGcccc                                                                                                             >  3:54848/1‑137 (MQ=255)                                                 acGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGcccc                                                                                                             <  4:54848/137‑1 (MQ=255)                                                 acGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATccc                                                                                                    <  2:164190/146‑1 (MQ=255)                                                 acGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATccc                                                                                                    >  1:164190/1‑146 (MQ=255)                                                       ttGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTAc                                                                                                                                           <  6:260020/101‑1 (MQ=255)                                                       ttGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTAc                                                                                                                                           >  5:260020/1‑101 (MQ=255)                                                       ttGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATa                                                                                           <  6:275226/149‑1 (MQ=255)                                                         gATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATAc                                                                                          <  7:251375/148‑1 (MQ=255)                                                              aCAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCGCCCAAAATCCCCCTGCCATACCCATTc                                                                                    <  2:199865/149‑1 (MQ=255)                                                                 aCCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCa                                                                                 <  2:111469/149‑1 (MQ=255)                                                                  ccTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTAACTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAg                                                                                >  5:276271/1‑149 (MQ=255)                                                                      tACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCaaaa                                                                                                        <  6:147091/121‑1 (MQ=255)                                                                      tACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCaaaa                                                                                                        >  5:147091/1‑121 (MQ=255)                                                                                 ggAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGcc                                                                 >  3:213425/1‑149 (MQ=255)                                                                                 ggAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGc                                                                  <  6:51961/148‑1 (MQ=255)                                                                                     gagaagaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGc                                                                  >  6:272873/1‑144 (MQ=255)                                                                                     gagaagaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGc                                                                  <  5:272873/144‑1 (MQ=255)                                                                                          gaCCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCa                                                                                                                    <  1:91822/89‑1 (MQ=255)                                                                                            ccGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGc                                                      <  2:125785/149‑1 (MQ=255)                                                                                                    cGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGa                                              <  8:213359/149‑1 (MQ=255)                                                                                                     ggTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGGTTTCCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGGCTTTGCCTATAATGGTGCCAGCCTGaa                                             >  7:184495/1‑149 (MQ=255)                                                                                                             gCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTATCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCaa                                     <  6:276271/149‑1 (MQ=255)                                                                                                               tGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAAc                                    >  3:159609/1‑148 (MQ=255)                                                                                                               tGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAAc                                    >  6:25048/1‑148 (MQ=255)                                                                                                               tGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATACCCCTGCCCTACCCATTCCCAGCCCGTGGCCTTTGCCTAATATGGTGCCAGCCTGAACCCGGCAACg                                   >  4:198087/1‑149 (MQ=255)                                                                                                                 ccTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCaa                                     >  2:155386/1‑145 (MQ=255)                                                                                                                 ccTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCaa                                     <  1:155386/145‑1 (MQ=255)                                                                                                                     gCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGcc                             >  3:52397/1‑149 (MQ=255)                                                                                                                      cAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGcc                             <  4:213425/148‑1 (MQ=255)                                                                                                                         cccATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCgtgt                         >  1:200736/1‑149 (MQ=255)                                                                                                                            aTTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTc                        <  5:88913/147‑1 (MQ=255)                                                                                                                            aTTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTc                        >  6:88913/1‑147 (MQ=255)                                                                                                                              ttgcaagCAAGCCCGCCCCTACAGGTAAACGTGCAGGAGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAg                    >  4:171040/7‑149 (MQ=37)                                                                                                                                      aaGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTcacc            >  1:57672/1‑149 (MQ=255)                                                                                                                                       aGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGc                                       >  2:113255/1‑121 (MQ=255)                                                                                                                                       aGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGc                                       <  1:113255/121‑1 (MQ=255)                                                                                                                                         cccGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCCCGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACg         >  1:13928/1‑149 (MQ=255)                                                                                                                                                cTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAg  <  4:52397/149‑1 (MQ=255)                                                                                                                                                cTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAg  <  3:788/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ATCAATGATCGAGACCTGGTGGTGATTGCGGCGGCTTACCTGGGTAACACGCGGTTGATCGACAACCTCTACCTGCATTTGGAAGAGAAGACCGCATAACGGTCTTGGGCTGCCTTGCAGCCCATTCGCGGGCAAGCCCGCTCCTGCACGTTTACCTGTAGGGGCGGGCTTGCTTGCAATGCCCCCAAAATCCCCCTGCCATACCCATTCCCAGCACGTGGCCTTTGCCTATAATGGTGCCAGCCTGAACCCGGCAACGACTGCCGTGTCCAAGCCCTCACCACGCACCAAG  >  NC_002947/5342456‑5342747

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