Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,928	G→T	6.2%	R87R (CGC→CGA)	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,928	0	G	T	6.2%	204.5 / 3.3	65	R87R (CGC→CGA)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base T (2/2);  total (36/29)
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.28e-01

CCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTCCCAAAGCCCGAAGCTG  >  NC_002947/1235780‑1236069                                                                                                                                                     |                                                                                                                                              ccGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTg                                                                                                                                               >  6:251962/1‑149 (MQ=255)   gAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTgcg                                                                                                                                             <  1:356893/149‑1 (MQ=255)       ccAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTgc                                                                                                                                              >  8:42388/1‑144 (MQ=255)       ccAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTgc                                                                                                                                              <  7:42388/144‑1 (MQ=255)       ccAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTTCGCTTg                                                                                                                                         >  3:377765/1‑149 (MQ=255)       ccAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTTCGCTTg                                                                                                                                         >  3:371767/1‑149 (MQ=255)        cAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGt                                                                                                                                        <  5:251962/149‑1 (MQ=255)           cTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCATGTCg                                                                                                                                      >  3:71796/1‑148 (MQ=255)               aGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAAccc                                                                                                                                 >  3:250993/1‑149 (MQ=255)                ggCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCt                                                                                                                                <  1:151126/149‑1 (MQ=255)                  cTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGTCAGCTTGCGCTTGTCGAACCCTGc                                                                                                                              <  7:57749/149‑1 (MQ=255)                      cTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGaa                                                                                                                           >  4:256809/1‑148 (MQ=255)                          gCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAgt                                                                                                                       <  3:256809/148‑1 (MQ=255)                              ccGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGa                                                                                                                   <  6:88029/148‑1 (MQ=255)                              ccGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAg                                                                                                                  <  7:314207/149‑1 (MQ=255)                               cgcgCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTACGCTTGTCGAACCCTGCGAACAGTGTGAgg                                                                                                                 <  8:122218/149‑1 (MQ=255)                                 cgcAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGc                                                                                                               >  7:10640/1‑149 (MQ=255)                                  gcagcaTGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCg                                                                                                              <  8:10640/149‑1 (MQ=255)                                       aTGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACt                                                                                                         >  5:120907/1‑149 (MQ=255)                                         ggCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGACAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCa                                                                                                       >  2:229059/1‑149 (MQ=255)                                                       cGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACAAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGTGCGATACTCACGGCTGGGTTCAgg                                                                                         >  6:337880/1‑149 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTTCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        <  4:377765/149‑1 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTTCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        <  4:371767/149‑1 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCt                                                                                                                                >  7:17619/1‑109 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCt                                                                                                                                <  8:17619/109‑1 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        >  7:319189/1‑149 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        >  4:396186/1‑149 (MQ=255)                                                        gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        >  3:295297/1‑149 (MQ=255)                                                            tCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTgcgc                                                                                    <  6:4041/149‑1 (MQ=255)                                                             cGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCa                                                                                   >  6:419274/1‑149 (MQ=255)                                                              gaagaaCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTCGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCAc                                                                                  >  7:52630/1‑149 (MQ=255)                                                                 gaaCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACTGCTGGGTTCAGGTGCGCACCg                                                                                >  8:334925/1‑148 (MQ=255)                                                                 gaaCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        <  1:155691/140‑1 (MQ=255)                                                                 gaaCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                                        >  2:155691/1‑140 (MQ=255)                                                                   aCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGa                                                                              <  4:295297/148‑1 (MQ=255)                                                                    cAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGaaa                                                                            >  4:122831/1‑149 (MQ=255)                                                                              tACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGaa                                                                             <  6:339625/138‑1 (MQ=255)                                                                              tACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGaa                                                                             >  5:339625/1‑138 (MQ=255)                                                                              tACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACcggcg                                                                    <  4:250993/147‑1 (MQ=255)                                                                                 aGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAg                                                               >  1:155038/1‑149 (MQ=255)                                                                                   cGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAgg                                                              >  2:148607/1‑148 (MQ=255)                                                                                       tAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGa                                                          <  1:148607/148‑1 (MQ=255)                                                                                         gACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGACAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCTCCGGAAATACCGGCGGTGAGGTAGATgg                                                       <  1:229059/149‑1 (MQ=255)                                                                                           ccAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGc                                                      >  7:186665/1‑148 (MQ=255)                                                                                            cAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCa                                                    >  2:184759/1‑149 (MQ=255)                                                                                               tgctgcGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGc                                                 <  4:258575/149‑1 (MQ=255)                                                                                                 ctgcGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTg                                               <  8:186665/149‑1 (MQ=255)                                                                                                  tgcGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGa                                              >  5:89986/1‑149 (MQ=255)                                                                                                  tgcGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGa                                              >  4:189201/1‑149 (MQ=255)                                                                                                    cgcCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacg                                            >  4:26543/1‑149 (MQ=255)                                                                                                    cgcCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacg                                            >  2:329432/1‑149 (MQ=255)                                                                                                    cgcCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacg                                            >  4:189076/1‑149 (MQ=255)                                                                                                      ccACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcc                                          >  7:139465/1‑149 (MQ=255)                                                                                                       cacaGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgccc                                         <  6:120907/149‑1 (MQ=255)                                                                                                           gAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccacgc                                     >  4:267872/1‑149 (MQ=255)                                                                                                            aaTGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccacgcc                                    >  4:382176/1‑149 (MQ=255)                                                                                                            aaTGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccacgcc                                    >  8:199808/1‑149 (MQ=255)                                                                                                              tGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCATGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGcccc                                  <  4:71796/149‑1 (MQ=255)                                                                                                                cgccgcATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCa                                 <  3:26543/148‑1 (MQ=255)                                                                                                                   cgcATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAgat                              >  3:296001/1‑148 (MQ=255)                                                                                                                    gcATACCTGGGCCAGCATGTAGAATGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAgatga                            <  5:419274/149‑1 (MQ=255)                                                                                                                     cATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAgatgat                           <  5:337880/149‑1 (MQ=255)                                                                                                                         ccTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAga                               <  6:33857/141‑1 (MQ=255)                                                                                                                         ccTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAga                               >  5:33857/1‑141 (MQ=255)                                                                                                                          cTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGa                      >  5:40979/1‑149 (MQ=255)                                                                                                                                cAGCATGTAGAAGGGCAGCTCGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTccc                <  8:52630/149‑1 (MQ=255)                                                                                                                                    aTGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTCCCAAAg            <  3:122831/149‑1 (MQ=255)                                                                                                                                     tGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCAATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTCCCAAAGc           >  1:93363/1‑149 (MQ=255)                                                                                                                                              gCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTCCCAAAGCCCGAAGCTg  <  2:155038/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTGATTTCCCAAAGCCCGAAGCTG  >  NC_002947/1235780‑1236069

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