Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,573,593	T→G	6.9%	C117G (TGC→GGC) ‡	PP_2257 →	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,573,593	0	T	G	6.9%	226.7 / 6.6	73	C117G (TGC→GGC) ‡	PP_2257	aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base G (3/2);  total (40/33)
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.67e-01

GTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTGGGGG  >  NC_002947/2573445‑2573740                                                                                                                                                     |                                                                                                                                                    gTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTAt                                                                                                                                                     <  2:212941/149‑1 (MQ=255)  ttGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATg                                                                                                                                                    >  1:351137/1‑149 (MQ=255)  ttGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATg                                                                                                                                                    >  1:360388/1‑149 (MQ=255)          gcgcTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcg                                                                                                                                             <  1:3712/148‑1 (MQ=255)            gcTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTAcgcgc                                                                                                                                          >  7:226017/1‑149 (MQ=255)                  tGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGTGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGGGGGTGACGAGTCGGTGCGCGTCTGCCCTACGCGCGAGc                                                                                                                                      >  4:276230/1‑147 (MQ=255)                   gAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAg                                                                                                                                   >  6:64129/1‑149 (MQ=255)                    aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                                                                  <  8:226017/149‑1 (MQ=255)                    aaaGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGt                                                                                                                                  <  5:64129/149‑1 (MQ=255)                         cGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGccc                                                                                                                             >  8:365196/1‑149 (MQ=255)                          ggCAAGAGCTGGATGGGTATCGTCAAGAACCCCTCCCGTAACGGCGACAACTACTGGGTCAATGCCAACGTCACGCCCATCCTCGAAGGAGGCCTGGTAGTGGGTTACGAGTCGCTGAGCGTAAGCCCTACCCGCGAGCAAGTAGcact                                                                                                                            >  6:307434/1‑146 (MQ=255)                           gCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGt                                                                                                                           >  2:230588/1‑149 (MQ=255)                           gCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGt                                                                                                                           <  5:71524/149‑1 (MQ=255)                             aaGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTAGTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGGGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                         >  6:114725/1‑149 (MQ=255)                                 gCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGaa                                                                                                                     <  6:367279/149‑1 (MQ=255)                                 gCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGaa                                                                                                                     <  5:243087/149‑1 (MQ=255)                                      aTGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGTGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCt                                                                                                                <  4:149616/149‑1 (MQ=255)                                       tGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTg                                                                                                               >  2:144147/1‑149 (MQ=255)                                        gggTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTg                                                                                                               >  3:56900/1‑148 (MQ=255)                                              cGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATgcgc                                                                                                        >  6:333822/1‑149 (MQ=255)                                               gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCTCGAGCAAGTAGCCCGTGcc                                                                                                                        >  1:309834/1‑132 (MQ=255)                                               gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCTCGAGCAAGTAGCCCGTGcc                                                                                                                        <  2:309834/132‑1 (MQ=255)                                               gTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATgcgc                                                                                                        <  8:168028/148‑1 (MQ=255)                                                  aaGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCt                                                                                                    >  2:129919/1‑149 (MQ=255)                                                       ccGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTgcgc                                                                                                >  7:62760/1‑148 (MQ=255)                                                         gCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAAc                                                                                             >  4:178900/1‑149 (MQ=255)                                                             ccGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCGAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGc                                                                                          >  4:138046/1‑148 (MQ=255)                                                             ccGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGcc                                                                                         <  8:90369/149‑1 (MQ=255)                                                              cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                         >  5:258552/1‑116 (MQ=255)                                                              cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGc                                                                                                                         <  6:258552/116‑1 (MQ=255)                                                              cGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGc                                                                                          >  5:5226/1‑147 (MQ=255)                                                                    ggCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAggcggc                                                                                  <  8:62760/149‑1 (MQ=255)                                                                     gCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCa                                                                                 <  3:346897/149‑1 (MQ=255)                                                                      cGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCAt                                                                                <  7:365196/149‑1 (MQ=255)                                                                             tactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACgg                                                                                           <  3:265665/131‑1 (MQ=255)                                                                             tactGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACgg                                                                                           >  4:265665/1‑131 (MQ=255)                                                                               ctGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcg                                                                        <  1:230588/148‑1 (MQ=255)                                                                                tGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgcc                                                                      >  2:233586/1‑149 (MQ=255)                                                                                  ggTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGgcgccgc                                                                    >  1:177225/1‑149 (MQ=255)                                                                                           cTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCt                                                             >  1:371911/1‑147 (MQ=255)                                                                                           cTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGa                                                           <  1:233586/149‑1 (MQ=255)                                                                                              cGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGc                                                         <  2:360388/148‑1 (MQ=255)                                                                                              cGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGc                                                         <  2:351137/148‑1 (MQ=255)                                                                                                  ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCt                                                    >  7:30186/2‑149 (MQ=255)                                                                                                  ccGCCCACCTTCGAGGATGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCt                                                    <  8:30186/148‑1 (MQ=255)                                                                                                      ccATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcg                                                >  3:328500/1‑149 (MQ=255)                                                                                                      ccATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcg                                                >  4:157482/1‑149 (MQ=255)                                                                                                       cATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGgcgc                                               <  5:114725/149‑1 (MQ=255)                                                                                                         tCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCGAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAg                                             <  3:138046/149‑1 (MQ=255)                                                                                                            tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                          >  7:342457/1‑149 (MQ=255)                                                                                                            tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                          >  6:354746/1‑149 (MQ=255)                                                                                                            tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggcc                                          >  3:360982/1‑149 (MQ=255)                                                                                                            tCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAggc                                           >  4:342624/1‑148 (MQ=255)                                                                                                               atggTGGGCGGGTTGTGTGTTACGAGTCGGTGCGCGTTTGCCCTTCGCTCGATCAAGTAGCCGGTGCCGAAGCCCTGTATGCGCGTCTGAGATACGGCGAGGCGGCATTGTCGGGGCCGCGCCTGTTGAGCTTGCTGGCGCAggccttg                                       <  5:307434/147‑1 (MQ=255)                                                                                                                     ggCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGgccttggccttg                                 >  7:114940/1‑149 (MQ=255)                                                                                                                      gCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCa                                              >  7:200085/1‑135 (MQ=255)                                                                                                                      gCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCa                                              <  8:200085/135‑1 (MQ=255)                                                                                                                        gggTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCa                              <  3:276230/149‑1 (MQ=255)                                                                                                                             gTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTgg                          <  6:5226/148‑1 (MQ=255)                                                                                                                                ggTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCg                       <  1:129919/148‑1 (MQ=255)                                                                                                                                  ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    >  2:137916/1‑149 (MQ=255)                                                                                                                                  ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    >  5:57168/1‑149 (MQ=255)                                                                                                                                  ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    >  6:43590/1‑149 (MQ=255)                                                                                                                                  ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    <  2:371911/149‑1 (MQ=255)                                                                                                                                  ttACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    <  3:157482/149‑1 (MQ=255)                                                                                                                                        gTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGgctgc              >  3:74016/1‑149 (MQ=255)                                                                                                                                            tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    <  6:109326/137‑1 (MQ=255)                                                                                                                                            tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    >  5:109326/3‑139 (MQ=255)                                                                                                                                            tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    >  5:114678/3‑139 (MQ=255)                                                                                                                                            tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTg                    <  6:114678/137‑1 (MQ=255)                                                                                                                                            tcgcgcGTAGGGCATACGCGCGAGCAAGTAGCCCATGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGc          >  2:101697/3‑149 (MQ=255)                                                                                                                                                 cgTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTgg     <  1:144147/149‑1 (MQ=255)                                                                                                                                                  gTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTggg    <  6:57168/149‑1 (MQ=255)                                                                                                                                                    aTGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTggggg  >  3:339639/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTTGATGTGGCGCTACCTGAAAGCCGGCAAGAGCTGGATGGGTATCGTCAAGAACCGCTGCCGTAACGGCGACTACTACTGGGTCAATGCCTACGTCACGCCCATCCTCGAAGGTGGGCGGGTAGTGGGTTACGAGTCGGTGCGCGTATGCCCTACGCGCGAGCAAGTAGCCCGTGCCGAAGCCCTGTATGCGCGGCTGCGCAACGGCCAGGCGGCATTGTCGGCGCCGCGCCGGCTGAGCTTGCTGGCGCAGGCCTTGGCCTTGCCACTGGCCGGTGGGCTGCTGGCGGTGGGGG  >  NC_002947/2573445‑2573740

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