Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,096,649	T→G	5.5%	S86A (TCG→GCG) ‡	PP_5608 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,096,649	0	T	G	5.5%	265.1 / 4.5	86	S86A (TCG→GCG) ‡	PP_5608	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/38);  new base G (2/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATAGCCGCCGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAGCGATCTTC  >  NC_002947/4096508‑4096793                                                                                                                                              |                                                                                                                                                 aTAGCCGCCGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCt                                                                                                                                           >  3:194471/1‑149 (MQ=255)  tAGCCGCCGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCGAGCGGCtg                                                                                                                                          <  3:214584/149‑2 (MQ=255)   aGCCGCCGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTc                                                                                                                                          >  7:322380/1‑148 (MQ=255)    gccgccGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAACAAACCTACAACAGCCGCTCGCGGCTCaa                                                                                                                                        <  1:220552/149‑1 (MQ=255)     ccgccgAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAg                                                                                                                                       >  2:139951/1‑149 (MQ=255)       gccgAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGc                                                                                                                                      >  5:337731/1‑148 (MQ=255)         cgAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCcag                                                                                                                                   <  5:143300/149‑1 (MQ=255)          gAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCcagc                                                                                                                                  <  4:61401/149‑1 (MQ=255)              cccAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCagcagc                                                                                                                               >  6:3947/1‑148 (MQ=255)               ccAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGa                                                                                                                             <  3:160761/149‑1 (MQ=255)                 aaCAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGAt                                                                                                                            >  5:252690/1‑148 (MQ=255)                   cAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCaa                                                                                                                         <  7:343768/149‑1 (MQ=255)                    aGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGACTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAg                                                                                                                        <  3:330764/149‑1 (MQ=255)                      gCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAg                                                                                                                        >  8:179444/1‑147 (MQ=255)                      gCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAg                                                                                                                        >  4:10215/1‑147 (MQ=255)                      gCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAg                                                                                                                        >  1:315005/1‑147 (MQ=255)                       ccAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAg                                                                                                                        >  7:330389/1‑146 (MQ=255)                        cAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGc                                                                                                                                      <  2:264807/131‑1 (MQ=255)                        cAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGc                                                                                                                                      >  1:264807/1‑131 (MQ=255)                        cAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGa                                                                                                                    >  5:20477/1‑149 (MQ=255)                               tCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTAc                                                                                                               <  1:371937/147‑1 (MQ=255)                                ccGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACtt                                                                                                             <  2:176761/148‑1 (MQ=255)                                ccGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACtt                                                                                                             <  5:29921/148‑1 (MQ=255)                                 cGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGTGCAAACCTACAACAGCCGCTCGCGGCTCAAGc                                                                                                                                      <  1:237748/122‑1 (MQ=255)                                 cGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGc                                                                                                                                      >  2:237748/1‑122 (MQ=255)                                 cGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTc                                                                                                            <  2:315005/148‑1 (MQ=255)                                   gggTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAg                                                                                                         <  5:160497/149‑1 (MQ=255)                                          ccTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCg                                                                                                  <  8:322380/149‑1 (MQ=255)                                           cTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGt                                                                                                 <  3:10215/149‑1 (MQ=255)                                           cTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGt                                                                                                 >  8:323655/1‑149 (MQ=255)                                            tGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGt                                                                                                 >  6:18897/1‑148 (MQ=255)                                                cGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGc                                                                                            <  6:20477/149‑1 (MQ=255)                                                  ccACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCAt                                                                                          >  4:146490/1‑149 (MQ=255)                                                   cACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATg                                                                                         >  7:311150/1‑149 (MQ=255)                                                   cACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATg                                                                                         <  8:311150/149‑1 (MQ=255)                                                      gcgcAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCg                                                                                      <  4:250291/149‑1 (MQ=255)                                                      gcgcAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCg                                                                                      >  4:253926/1‑149 (MQ=255)                                                          aagaagCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTg                                                                                   >  6:68950/1‑148 (MQ=255)                                                          aagaagCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCGAGCGGCtgtt                                                                                                                                        <  3:349817/95‑4 (MQ=255)                                                          aagaagCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCGAGCGGCtgtt                                                                                                                                        >  4:349817/1‑92 (MQ=255)                                                             aagCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCCTGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGc                                                                                  >  6:284391/1‑146 (MQ=255)                                                               gCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGAc                                                                             <  7:179444/149‑1 (MQ=255)                                                                   cAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGa                                                                           <  6:337731/147‑1 (MQ=255)                                                                      gTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAAc                                                                       >  2:264454/1‑148 (MQ=255)                                                                         cGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAActct                                                                    <  5:68950/148‑1 (MQ=255)                                                                            aaaTTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTAc                                                                >  5:150003/1‑149 (MQ=255)                                                                                  cccAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATc                                                           <  5:18897/148‑1 (MQ=255)                                                                                  cccAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCa                                                          >  5:330259/1‑149 (MQ=255)                                                                                  cccAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGAATGCCGTTGCAAGACGATAACTCATTACTTATCa                                                          >  1:277149/1‑149 (MQ=255)                                                                                   ccAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCAc                                                         >  8:75657/1‑149 (MQ=255)                                                                                   ccAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCAc                                                         >  1:126919/1‑149 (MQ=255)                                                                                     aaaCCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGc                                                       <  6:33690/149‑1 (MQ=255)                                                                                       aCCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCt                                                     >  7:199844/1‑149 (MQ=255)                                                                                       aCCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCt                                                     <  3:7679/149‑1 (MQ=255)                                                                                        cccTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTg                                                    >  6:134660/1‑149 (MQ=255)                                                                                             gATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCaaa                                                <  6:252690/148‑1 (MQ=255)                                                                                             gATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAg                                               >  7:42983/1‑149 (MQ=255)                                                                                               tGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGc                                             <  5:3947/149‑1 (MQ=255)                                                                                               tGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGc                                             <  1:139951/149‑1 (MQ=255)                                                                                                  gcATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCg                                           >  7:3878/1‑148 (MQ=255)                                                                                                    aTGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCAt                                        >  3:42399/1‑149 (MQ=255)                                                                                                     tGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATc                                       <  2:126919/149‑1 (MQ=255)                                                                                                       aCTGCCTTCGGCATTACCGAGCTAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCa                                      <  1:264454/148‑1 (MQ=255)                                                                                                           ccTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCTATCAAGCAGACGTACTTCGAGCAGTGCGTTCACCATGCCGTTGCAAGACTATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGt                                  <  5:377213/148‑1 (MQ=255)                                                                                                                   aTTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGtc                                 <  4:164322/141‑1 (MQ=255)                                                                                                                   aTTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGtc                                 >  3:164322/1‑141 (MQ=255)                                                                                                                      aCCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAgg                       <  7:323655/148‑1 (MQ=255)                                                                                                                        cGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGGAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgcc                    >  7:198760/1‑149 (MQ=255)                                                                                                                        cGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgcc                    >  7:119767/1‑149 (MQ=255)                                                                                                                           gCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgcctgc                 <  8:119767/149‑1 (MQ=255)                                                                                                                            cAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgcctgcc                <  5:134660/149‑1 (MQ=255)                                                                                                                            cAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgcctgcc                <  4:42399/149‑1 (MQ=255)                                                                                                                              aaCCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCATCGTCTCCATGGAGGGCCTGCCtt              >  3:5310/1‑149 (MQ=255)                                                                                                                              aaCCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgc                     <  4:237750/142‑1 (MQ=255)                                                                                                                              aaCCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCtt              >  4:357547/1‑149 (MQ=255)                                                                                                                              aaCCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGATGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGgccagatcc              >  3:237750/1‑143 (MQ=255)                                                                                                                                 cTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTg             >  4:134849/1‑147 (MQ=255)                                                                                                                                 cTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGaa           <  6:330259/149‑1 (MQ=255)                                                                                                                                  tACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAg          >  4:294804/1‑149 (MQ=255)                                                                                                                                  tACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACCCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAg          >  7:68607/1‑149 (MQ=255)                                                                                                                                    caacaGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAg          <  6:150003/147‑1 (MQ=255)                                                                                                                                     ttgaGCCGCGAGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGaa           <  7:50600/142‑1 (MQ=255)                                                                                                                                     ttgaGCCGCGAGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGaa           >  8:50600/4‑145 (MQ=255)                                                                                                                                       caGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAGCGATc     >  6:134090/1‑149 (MQ=255)                                                                                                                                       caGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAGCGATc     >  1:161269/1‑149 (MQ=255)                                                                                                                                          ccGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAGCGATCTTc  >  1:228202/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                 ATAGCCGCCGAGGCCCAACAGGCCAAGCTTTCCGGGGTAACCCTGTACGCCACGCGCAAGAAGCTGCAGGTACGCAAATTCCCCAAACCCTGGATGCGCATGACTGCCTTCGGCATTACCGAGCAAACCTACAACAGCCGCTCGCGGCTCAAGCCAGCAGCGATCAAGCAGACGTACTTCGAGCAGTGCGTTCAGCATGCCGTTGCAAGACGATAACTCTTTACTTATCACGCCTGCAAAGGCCGCATCAGCGTCTCCATGGAGGGCCTGCCTTGAAGCGATCTTC  >  NC_002947/4096508‑4096793

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