Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,655,636	C→A	6.6%	intergenic (‑91/+61)	epd ← / ← tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,655,636	0	C	A	6.6%	164.3 / 4.0	62	intergenic (‑91/+61)	epd/tktA	D‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base A (2/2);  total (32/30)
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.65e-02

CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCC  >  NC_002947/5655488‑5655781                                                                                                                                                          |                                                                                                                                                  cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCAGCTCCCACCG‑GGGTTCCTGc                                                                                                                                                   >  6:240330/1‑149 (MQ=255)    aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCTGGGCTCACTAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGAGGTGTGT‑‑‑‑CATTAGCGGGCAAGCCAGCTCCCACAT‑GGTTTCATTCAAc                                                                                                                                                >  7:124167/1‑149 (MQ=255)      cGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACct                                                                                                                                              <  5:66307/149‑1 (MQ=255)        gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACctct                                                                                                                                            >  3:253668/1‑149 (MQ=255)          agATACGGACATATACGTTGAGTGGAACTTTTTTGGGACGGTTGGGGATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAATCCCGCTCCCTCAG‑TGTTTCATGCAACATCtgt                                                                                                                                          <  4:34225/148‑1 (MQ=255)                ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTGATGCAACCTCTGTGGGAGc                                                                                                                                    >  8:148835/1‑149 (MQ=255)                 gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCg                                                                                                                                   <  1:74538/149‑1 (MQ=255)                  ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGa                                                                                                                                      <  4:34647/149‑1 (MQ=21)                    aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAAc                                                                                                                                                <  2:5325/133‑1 (MQ=255)                    aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAAc                                                                                                                                                >  1:5325/1‑133 (MQ=255)                       tcagtgtaGTGCAACTTTGTAGGGACGGTGGGGCTTTCGTGGTTCACGAGACAGTTCAATTGTCGTCAGCATGAGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCAGAG‑GTTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                              <  8:124167/141‑1 (MQ=38)                       aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                              <  7:277326/148‑1 (MQ=255)                        ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt                                                                                                                                          >  7:109758/1‑135 (MQ=255)                        ccGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt                                                                                                                                          <  8:109758/135‑1 (MQ=255)                                tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCTGGCTTGCCCGCGAAt                                                                                                                    >  5:334431/1‑149 (MQ=255)                                   aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                   >  7:40200/1‑147 (MQ=255)                                        tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCaa                                                                                                                                                 <  3:83292/112‑1 (MQ=255)                                        tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCaa                                                                                                                                                 >  4:83292/1‑112 (MQ=255)                                        tGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCCGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGGTTCCACAG‑GGTTTCATGCAACCTCTGTGGGATCGGGCTTGCCCGCGAATGAAGGc                                                                                                              <  5:241277/147‑1 (MQ=255)                                          tAGGGACGGGGGGGCATTCGGGGCGCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGac                                                                                                           <  7:204274/148‑1 (MQ=255)                                                cGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGCTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgccc                                                                                                    <  5:240330/149‑1 (MQ=255)                                                       gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCGCCCCCAC‑GGGTTGATGGAACCACTGTGGGAGCCGGGGTGGCCGCGAATTAAgg                                                                                                               >  1:339691/1‑131 (MQ=255)                                                             ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCCCACCCACAG‑GGTTTAATGCAACCTCTGTGGGAGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATGGTTaa                                                                                        >  8:41029/1‑148 (MQ=255)                                                             ggggCTAACGAGACATTTCAATTGTCGTCTGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGGGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGGCACCGCCCGCAGATCTTTaaa                                                                                       <  8:40200/149‑1 (MQ=38)                                                              gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                      <  7:41029/149‑1 (MQ=255)                                                                    aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCACAGTGGGAGCGGGCTTTCCCGCGAAAGAATGCGCCACCGCCCGCAGATCGTTAAATAAGTCt                                                                                >  8:160897/1‑149 (MQ=38)                                                                    aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAACCCCCCACCCACAG‑GGGTTCATGGAACCTCTGTGGGACCCGGGATTCCCGCCAATGTAGGCGACACCGCCCGCAGATCTTTAAATCattgt                                                                                >  2:214148/1‑145 (MQ=255)                                                                      gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                              <  1:214148/149‑1 (MQ=255)                                                                       agaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCCCAG‑GGTTTCATACCCCCGCGGTGGGAGCGGGCTTGCCCGCGGATGACGGCGCCGCCGCCCGCGGAGCGGTAAAGGAGTCtggg                                                                             >  4:102550/1‑146 (MQ=255)                                                                             ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                                       >  5:333718/1‑149 (MQ=255)                                                                             ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                                       >  6:331782/1‑149 (MQ=255)                                                                             ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAAACCCCACCCACAG‑GGGTTTAAGGAAACGCTGTGGGAGCGGGGGTGGCCCCGGAATAAGGCGAACCCGCCCGCAGATCTTTTAATAATTCTTGcagcag                                                                        >  6:24945/1‑148 (MQ=255)                                                                                     cgttctgCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCTAGCCCGCTCCCACAG‑GGTCTCATGGAAGCTCTGTGTGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgt                                                                <  7:217104/142‑2 (MQ=38)                                                                                       gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                                <  2:172512/82‑1 (MQ=255)                                                                                       gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                              <  7:160897/148‑1 (MQ=255)                                                                                       gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                             >  6:342809/1‑149 (MQ=255)                                                                                       gtcAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAC‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGGCCGCGAATGAAGGCGACACCGCCCGCCGATCGTTAAATAAGTCTTCCAGCAGCTCGTCGGCAg                                                             >  6:1302/1‑149 (MQ=38)                                                                                          aGCCTGCGGAATTTGCGGTGGGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTa                                                           <  3:6127/148‑1 (MQ=255)                                                                                              tGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCCGCTCTTCGGCAGTACCGAg                                                      >  5:79002/1‑149 (MQ=255)                                                                                                cGGAATTTGCGGGGTGG‑‑‑‑CATTCGCGGTCAAGCCCGCTCCCACAG‑GGGTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGat                                                      <  6:334431/147‑2 (MQ=255)                                                                                                     tttGCGGGGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtc                                               <  5:24945/149‑1 (MQ=255)                                                                                                           gtgtGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         <  3:169023/149‑1 (MQ=255)                                                                                                              tgG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                                       <  5:331782/148‑1 (MQ=255)                                                                                                                      aTTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         >  6:346575/1‑142 (MQ=255)                                                                                                                      aTTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         <  5:346575/142‑1 (MQ=255)                                                                                                                          gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                              <  7:148835/149‑1 (MQ=255)                                                                                                                          gcgGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaaa                              <  4:253668/149‑1 (MQ=255)                                                                                                                                  gCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg                       >  3:104466/1‑148 (MQ=255)                                                                                                                                   cccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGa                      >  4:52096/1‑148 (MQ=255)                                                                                                                                    ccGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                     >  5:295871/1‑148 (MQ=255)                                                                                                                                       cTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt                 >  3:219671/1‑149 (MQ=255)                                                                                                                                        tCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAAT‑‑‑‑GCCACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg             >  2:226121/1‑148 (MQ=21)                                                                                                                                        tCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTATTCGAACAATg                >  5:317204/1‑149 (MQ=255)                                                                                                                                         cccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGc               >  3:138222/1‑149 (MQ=38)                                                                                                                                          ccACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                         <  2:122932/138‑1 (MQ=37)                                                                                                                                          ccACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc              <  5:342809/149‑1 (MQ=255)                                                                                                                                            acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTcc                                            >  1:294963/1‑117 (MQ=255)                                                                                                                                            acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTcc                                            <  2:294963/117‑1 (MQ=255)                                                                                                                                             caG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAACAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa           >  4:333608/1‑149 (MQ=255)                                                                                                                                              aG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         >  7:31472/1‑118 (MQ=255)                                                                                                                                              aG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                         <  8:31472/118‑1 (MQ=255)                                                                                                                                                       aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc  >  5:160013/1‑149 (MQ=255)                                                                                                                                                       aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc  >  7:173032/1‑149 (MQ=255)                                                                                                                                                          |                                                                                                                                                  CACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCC  >  NC_002947/5655488‑5655781

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