Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,385,793	Δ1 bp	6.1%	coding (85/576 nt)	PP_2093 ←	two‑component system response regulator NasT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,385,793	0	T	.	6.1%	207.6 / 7.2	66	coding (85/576 nt)	PP_2093	two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base T (34/28);  new base . (2/2);  total (36/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 = 6.82e-01

AACACGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAC  >  NC_002947/2385645‑2385941                                                                                                                                                     |                                                                                                                                                      aacaCGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCgat                                                                                                                                                       >  6:331600/1‑149 (MQ=255) aacaCGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCgat                                                                                                                                                       >  2:409142/1‑149 (MQ=255)     cGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGAtc                                                                                                                                                   <  4:54246/149‑1 (MQ=255)      gATGGGCCGCGGCTGGTCACGGCTGACCCGCACTACCTGTTCCATCACATAGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGAAGCAGGCGGGCAGGTCGATGGTCATGTCTGGGGCTTCGATGAtc‑t                                                                                                                                                 >  1:238515/1‑149 (MQ=255)         gggCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCg                                                                                                                                               >  7:17904/1‑148 (MQ=255)          ggCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGaa                                                                                                                                             <  5:242117/149‑1 (MQ=255)            ccGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGc                                                                                                                                           <  1:409142/149‑1 (MQ=255)             cgcgGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCTCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGA‑GATCATCGAAGcc                                                                                                                                          <  7:264755/149‑1 (MQ=255)             cgcgGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGA‑GATCATCGAAGcc                                                                                                                                          <  5:255649/149‑1 (MQ=255)                      tCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCctcact                                                                                                                                 >  3:331545/1‑149 (MQ=255)                                       tACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGc                                                                                                                >  5:418346/1‑149 (MQ=255)                                        aCCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCg                                                                                                               >  6:189506/1‑149 (MQ=255)                                           tGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGAcc                                                                                                            >  3:342337/1‑149 (MQ=255)                                             ttCCATCACATCGCGGTCT‑GCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAg                                                                                                         <  5:66401/149‑1 (MQ=255)                                                      aTCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGctct                                                                                                 >  5:359575/1‑149 (MQ=255)                                                       tCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTg                                                                                                >  2:6776/1‑149 (MQ=255)                                                          cgGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGc                                                                                                        >  5:40797/1‑138 (MQ=255)                                                          cgGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGc                                                                                                        <  6:40797/138‑1 (MQ=255)                                                                 gCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTgtcgtc                                                                                       <  6:418346/148‑1 (MQ=255)                                                                   gAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCAATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGAt                                                                                    >  1:15157/1‑149 (MQ=255)                                                                        ggTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                               <  5:121563/149‑1 (MQ=255)                                                                        ggTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                               <  3:154080/149‑1 (MQ=255)                                                                            tCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGa                                                                            <  5:331600/148‑1 (MQ=255)                                                                            tCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                                                                           <  2:86064/149‑1 (MQ=255)                                                                            tCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGAt                                                                           >  4:389641/1‑149 (MQ=255)                                                                               aTCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                               <  2:74207/142‑1 (MQ=255)                                                                               aTCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                               >  1:74207/1‑142 (MQ=255)                                                                               aTCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATgcg                                                                        >  3:144524/1‑149 (MQ=255)                                                                                    caccacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACa                                                                   <  6:411004/149‑1 (MQ=255)                                                                                     accacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACAt                                                                  >  8:362981/1‑149 (MQ=255)                                                                                     accacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAAAAt                                                                  >  6:403905/1‑149 (MQ=255)                                                                                      ccacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATg                                                                 >  3:227553/1‑149 (MQ=255)                                                                                       cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGt                                                                >  2:174475/1‑149 (MQ=255)                                                                                       cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGt                                                                >  6:401610/1‑149 (MQ=255)                                                                                         cGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTc                                                              <  8:17904/149‑1 (MQ=255)                                                                                          gTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCAATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCg                                                             <  2:15157/149‑1 (MQ=255)                                                                                              gggCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTccc                                                         >  6:379713/1‑149 (MQ=255)                                                                                                    aCCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCAc                                                   >  4:396067/1‑149 (MQ=255)                                                                                                      cGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCAcgc                                                 >  7:362593/1‑149 (MQ=255)                                                                                                        tttCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGAt                                                                                    <  2:369901/112‑1 (MQ=255)                                                                                                        tttCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGAt                                                                                    >  1:369901/1‑112 (MQ=255)                                                                                                        tttCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGCCCCCCTCACgcgc                                               >  4:357185/1‑149 (MQ=255)                                                                                                             aCGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGaca                                          >  1:386487/1‑149 (MQ=255)                                                                                                              cGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGacac                                         <  6:233499/149‑1 (MQ=255)                                                                                                                     gggCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATc                                  <  4:138922/149‑1 (MQ=255)                                                                                                                         aGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGCCCCCCTCACGCGCCGACACGGGCATCGGGc                              <  3:357185/149‑1 (MQ=255)                                                                                                                           gTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTggg                                                                                              <  1:1894/83‑1 (MQ=255)                                                                                                                           gTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTggg                                                                                              >  2:1894/1‑83 (MQ=255)                                                                                                                           gTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCGCCTCACGCGCCGACACGGGCATCGGGCaa                            >  2:134608/1‑149 (MQ=255)                                                                                                                           gTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCaa                            >  4:161015/1‑149 (MQ=255)                                                                                                                            tCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAAc                           <  4:331545/149‑1 (MQ=255)                                                                                                                             cGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACt                          >  1:241122/1‑149 (MQ=255)                                                                                                                              gATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTc                         >  5:47844/1‑149 (MQ=255)                                                                                                                              gATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCCAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATTCGCAACATGTTCGTCCCCCTCACGCCCCGACACGGGCATCGTGCAACTc                         >  4:174627/1‑149 (MQ=255)                                                                                                                                tGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGc                       <  6:359575/149‑1 (MQ=255)                                                                                                                                   tCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATca                                                                                  >  8:147007/1‑87 (MQ=255)                                                                                                                                   tCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATca                                                                                  <  7:147007/87‑1 (MQ=255)                                                                                                                                     aGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCaa                   <  7:362981/148‑1 (MQ=255)                                                                                                                                       gTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAgg                 >  3:259228/1‑148 (MQ=255)                                                                                                                                        tCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCa               >  3:353821/1‑149 (MQ=255)                                                                                                                                           ccggCTTCGA‑GATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATg             >  5:45923/3‑148 (MQ=255)                                                                                                                                            gtgCTTCGATTATG‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGc            <  3:174627/146‑1 (MQ=255)                                                                                                                                            gggCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCAt              <  1:6776/146‑1 (MQ=255)                                                                                                                                            gggCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGc            >  5:165725/1‑148 (MQ=255)                                                                                                                                             ggCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCaa          <  4:342337/149‑1 (MQ=255)                                                                                                                                             ggCTTCGA‑GATCATCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGc            >  5:262564/1‑147 (MQ=255)                                                                                                                                              gCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGTTGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAg         <  4:144524/149‑1 (MQ=255)                                                                                                                                              gCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAg         <  3:161015/149‑1 (MQ=255)                                                                                                                                                   gatgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGa     <  5:379713/148‑1 (MQ=255)                                                                                                                                                   gatgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGa     >  3:12601/1‑148 (MQ=255)                                                                                                                                                     tgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAc  >  3:49802/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AACACGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAC  >  NC_002947/2385645‑2385941

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