Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,512,979	G→A	5.5%	I37I (ATC→ATT)	yegQ ←	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,512,979	0	G	A	5.5%	224.6 / 2.0	73	I37I (ATC→ATT)	yegQ	peptidase
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base A (2/2);  total (34/39)
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.75e-01

GATCAGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCCAGCAGTTCGGGC  >  NC_002947/2512832‑2513126                                                                                                                                                    |                                                                                                                                                    gATCAGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGa                                                                                                                                                    >  7:181876/1‑149 (MQ=255)  aTCAGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGa                                                                                                                                                    <  2:15572/148‑1 (MQ=255)    cAGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATc                                                                                                                                                  <  5:358110/148‑1 (MQ=255)     aGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATccc                                                                                                                                                <  1:87657/149‑1 (MQ=255)      ggCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATcccc                                                                                                                                               <  4:113035/149‑1 (MQ=255)       gCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATcccc                                                                                                                                               >  6:243378/1‑148 (MQ=255)       gCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGAGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGCGCGCCTCCTCGATcccc                                                                                                                                               >  8:289100/1‑148 (MQ=255)             gTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCa                                                                                                                                        <  3:169550/149‑1 (MQ=255)                ggACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGt                                                                                                                                     >  7:305442/1‑149 (MQ=255)                       aTCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTGTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTg                                                                                                                              >  6:40928/1‑149 (MQ=255)                       aTCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGt                                                                                                                               <  8:305442/148‑1 (MQ=255)                         cAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCAATCCCCAATGCCAGGTTGGCGTGAt                                                                                                                            >  4:90544/1‑149 (MQ=255)                         cAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGATTGGCATTGTGAGGGGCGATGTTGACCACCAAGTAGAAGCGCTTGCCCACGGCGTGCGCCTCCTCAATACCCAATGCCAGGTTGGCGTGAt                                                                                                                            >  4:85357/1‑149 (MQ=255)                          aGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGAt                                                                                                                            <  4:271035/148‑1 (MQ=255)                          aGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGAt                                                                                                                            <  4:266124/148‑1 (MQ=255)                          aGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGAt                                                                                                                            <  7:165512/148‑1 (MQ=255)                           gcgcGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCg                                                                                                                          <  2:108793/149‑1 (MQ=255)                            cgcgTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCg                                                                                                                          >  4:29093/1‑148 (MQ=255)                                  cGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCgt                                                                                                                   >  2:314752/1‑149 (MQ=255)                                   ggCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCgtt                                                                                                                  >  2:198265/1‑149 (MQ=255)                                           ctcGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCa                                                                                                          >  7:157626/1‑149 (MQ=255)                                            tcGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCAc                                                                                                         >  3:108151/1‑149 (MQ=255)                                              gATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGGTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCACTGCCAGGTTGGCGTGATCGAAGTCGTTGTTGCGCAccc                                                                                                       <  6:65188/149‑1 (MQ=255)                                               aTCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCg                                                                                                      >  5:90864/1‑149 (MQ=255)                                                   cGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTg                                                                                                                                   >  2:265734/1‑116 (MQ=255)                                                   cGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAACCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTg                                                                                                                                   <  1:265734/116‑1 (MQ=255)                                                    gggTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTTGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGc                                                                                                 <  5:120849/149‑1 (MQ=255)                                                      gTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTg                                                                                               <  1:198265/149‑1 (MQ=255)                                                              tCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGCGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAATCTGTAACGCgg                                                                                       <  7:289100/149‑1 (MQ=255)                                                                cTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCt                                                                                     >  7:95233/1‑149 (MQ=255)                                                                    aGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGc                                                                                  <  8:181876/148‑1 (MQ=255)                                                                      gAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCgg                                                                               >  3:338440/1‑149 (MQ=255)                                                                      gAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCgg                                                                               <  8:157626/149‑1 (MQ=255)                                                                      gAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCgg                                                                               <  8:176949/149‑1 (MQ=255)                                                                       aaGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCACTGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGc                                                                              <  7:244968/149‑1 (MQ=255)                                                                        aGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCg                                                                             >  7:93576/1‑149 (MQ=255)                                                                          gTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTa                                                                           <  8:93576/149‑1 (MQ=255)                                                                           tCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTgg                                                                                   <  7:335088/140‑1 (MQ=255)                                                                           tCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTgg                                                                                   >  8:335088/1‑140 (MQ=255)                                                                               gAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATc                                                                                                                                                  >  8:280846/1‑73 (MQ=255)                                                                               gAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATc                                                                                                                                                  <  7:280846/73‑1 (MQ=255)                                                                               gAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGAcc                                                                       <  5:243378/148‑1 (MQ=255)                                                                               gAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGAcc                                                                       <  1:178600/148‑1 (MQ=255)                                                                                 gCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCg                                                                             >  6:191304/1‑140 (MQ=255)                                                                                 gCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCg                                                                             <  5:191304/140‑1 (MQ=255)                                                                                     ggCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATc                                                                  <  5:58868/147‑1 (MQ=255)                                                                                     ggCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCgg                                                                >  8:210215/1‑149 (MQ=255)                                                                                        aTTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGAc                                                                        >  6:373391/1‑138 (MQ=255)                                                                                        aTTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGAc                                                                        <  5:373391/138‑1 (MQ=255)                                                                                               ggggCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCgg                                                                                       >  5:291631/1‑116 (MQ=255)                                                                                               ggggCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCgg                                                                                       <  6:291631/116‑1 (MQ=255)                                                                                                      tGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAAcgcg                                               >  5:120948/1‑149 (MQ=255)                                                                                                        ttGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAAcgcg                                               >  8:345561/1‑147 (MQ=255)                                                                                                        ttGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTa                                             <  8:95233/149‑1 (MQ=255)                                                                                                        ttGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTa                                             <  5:40928/149‑1 (MQ=255)                                                                                                           accaccACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAgcg                                          <  6:120948/149‑1 (MQ=255)                                                                                                               ccacGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATg                                      <  1:314752/149‑1 (MQ=255)                                                                                                               ccacGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATg                                      >  5:120481/1‑149 (MQ=255)                                                                                                                cacGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATgg                                     >  5:329387/1‑149 (MQ=255)                                                                                                                   gTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCt                                  >  5:164559/1‑149 (MQ=255)                                                                                                                      gAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCGAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGa                               <  3:11974/149‑1 (MQ=255)                                                                                                                      gAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCtt                                 >  7:70552/1‑147 (MQ=255)                                                                                                                      gAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGa                               <  4:338440/149‑1 (MQ=255)                                                                                                                       aaGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAg                              <  5:350051/149‑1 (MQ=255)                                                                                                                         gcgcTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAggg                            >  2:136236/1‑149 (MQ=255)                                                                                                                                  cAAGGCGTGCGCCTCCTCGCTCCCCAATGCCAGGTTGGCGTGATCGACCTCGTTGTTGCGCACCCGGAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCcggcgg                   >  5:112687/1‑149 (MQ=255)                                                                                                                                  cAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCcggcgg                   >  3:268848/1‑149 (MQ=255)                                                                                                                                  cAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCcggcgg                   >  1:249974/1‑149 (MQ=255)                                                                                                                                           cgcCTCCTCAATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTTGCCGGCGTAGACCGCAGCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCcagc           <  3:85357/148‑1 (MQ=255)                                                                                                                                           cgcCTCCTCAATCCCCAATGCCAGGGTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCcagc           <  3:90544/148‑1 (MQ=255)                                                                                                                                                cctcGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCCAGCAGTTCg     <  1:136236/149‑1 (MQ=255)                                                                                                                                                 ctcGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGAGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCCAGCAGTTCgg    >  1:72948/1‑149 (MQ=255)                                                                                                                                                   cGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCCAGCAGTTCGGGc  <  7:210215/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GATCAGGCCGGGGTCGGACATGATCAGCGCGTCCGGCGCCATCTCGATCACGGGTGCCAGGTCCTTGAGGAAGGTCTTGAGCTTGGCATTGTGCGGGGCGATGTTGACCACCACGTAGAAGCGCTTGCCCAAGGCGTGCGCCTCCTCGATCCCCAATGCCAGGTTGGCGTGATCGAACTCGTTGTTGCGCACCCGCAAGCTGTAACGCGGCTGGCCGGCGTAGACCGCATCGGCGCCGTAGGCAAACGCGTAGCGCATGGTCTTGAGGGTGCCGGCGGGGGCCAGCAGTTCGGGC  >  NC_002947/2512832‑2513126

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