Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,385,797:1	+A	5.8%	coding (81/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,797	1	.	A	5.8%	281.0 / 11.2	69	coding (81/576 nt)	PP_2093	two‑component system response regulator NasT
Reads supporting (aligned to +/- strand):  ref base . (36/29);  new base A (2/2);  total (38/31)
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.33e-01

CGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAC  >  NC_002947/2385649‑2385941                                                                                                                                                      |                                                                                                                                                 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                                                                               <  3:154080/149‑1 (MQ=255)                                                                    ggTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATcagca                                                                               <  5:121563/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                                                                >  6:401610/1‑149 (MQ=255)                                                                                   cacGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGt                                                                >  2:174475/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     >  3:12601/1‑148 (MQ=255)                                                                                                                                               gatgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGa     <  5:379713/148‑1 (MQ=255)                                                                                                                                                 tgatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAc  >  3:49802/1‑149 (MQ=255)                                                                                                                                                  gatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCAt              >  3:8629/1‑136 (MQ=255)                                                                                                                                                  gatC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCAt              <  4:8629/136‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 CGATGGGCCGCGGCTGGTCACGGCTGACCAGCACTACCTGTTCCATCACATCGCGGTCTGGCGAGTCGGTATCGATCAGCACCACGTCCGGGCGCACCGTTTCGACGCAGGCGGGCAGGTCGATGGTCAGGTCTGGGGCTTCGATGATC‑TCGAAGCCGGCCTCACTCAACGCAGCCCTCAGGCGACCGAGCTTGCTCTGGGTGTCGTCGATCAGCAGGATGCGCAACATGTTCGTCCCCCTCACGCGCCGACACGGGCATCGGGCAACTCGCCCAAGGCATGCAAGCTGAAAC  >  NC_002947/2385649‑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: ‑