Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	418,296	C→G	6.8%	Q27E (CAG→GAG)	PP_0345 →	heptose kinase WapP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	418,296	0	C	G	6.8%	224.1 / 4.0	72	Q27E (CAG→GAG)	PP_0345	heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base G (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCA  >  NC_002947/418155‑418441                                                                                                                                                |                                                                                                                                                  ctggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACATGGTGAGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCATCCAGGGc                                                                                                                                              >  1:342427/1‑149 (MQ=255)    gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGaa                                                                                                                                           <  5:296889/147‑1 (MQ=255)    gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAAc                                                                                                                                          <  6:97494/148‑1 (MQ=255)     ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAcac                                                                                                                                        <  1:213731/149‑3 (MQ=255)         aGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTga                                                                                                                                    <  2:221015/149‑1 (MQ=255)              aaGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAggg                                                                                                                                               <  6:26506/133‑1 (MQ=255)              aaGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAggg                                                                                                                                               >  5:26506/1‑133 (MQ=255)                   ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgaccc                                                                                                                          >  1:159842/1‑149 (MQ=255)                   ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgaccc                                                                                                                          >  8:11457/1‑149 (MQ=255)                   ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgaccc                                                                                                                          >  2:17309/1‑149 (MQ=255)                    cAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                         >  5:221458/1‑149 (MQ=255)                    cAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                         >  5:82111/1‑149 (MQ=255)                       aaaCTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCt                                                                                                                       >  4:197986/1‑148 (MQ=255)                           tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGCGTTCGCCCTGGAGGGCGAAcac                                                                                                                                        >  7:20488/1‑125 (MQ=255)                           tctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGCGTTCGCCCTGGAGGGCGAAcac                                                                                                                                        <  8:20488/127‑3 (MQ=255)                            ctACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTg                                                                                                                  >  7:97354/1‑148 (MQ=255)                               cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCGCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAg                                                                                                               <  6:221458/148‑1 (MQ=255)                                gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                             >  4:165028/1‑149 (MQ=255)                                     caagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                     >  4:66884/1‑136 (MQ=255)                                     caagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                     <  3:66884/136‑1 (MQ=255)                                           aCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGa                                                                                                  <  7:11457/149‑1 (MQ=255)                                            cGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAg                                                                                                 >  3:397067/1‑149 (MQ=255)                                              cTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCg                                                                                               <  4:125234/149‑1 (MQ=255)                                               tATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGt                                                                                              >  5:271501/1‑149 (MQ=255)                                                   gggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                                           <  3:165028/148‑1 (MQ=255)                                                   gggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                                           <  1:17309/148‑1 (MQ=255)                                                   gggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                                           <  6:271501/148‑1 (MQ=255)                                                    ggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACgg                                                                                         <  5:208557/149‑1 (MQ=255)                                                     gATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCGCAGTGAAGTGGTACGCATCGAGCGTGACgg                                                                                         >  6:41276/1‑148 (MQ=255)                                                          aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACCCATCGAGCGTGACGGGATCaa                                                                                   <  7:81352/149‑1 (MQ=255)                                                              tGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGAGACGGGATCAacta                                                                                >  7:179378/1‑148 (MQ=255)                                                                 tGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                                            <  2:342427/149‑1 (MQ=255)                                                                     gAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGt                                                                                                                   >  5:192575/1‑106 (MQ=255)                                                                     gAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGt                                                                                                                   <  6:192575/106‑1 (MQ=255)                                                                      aGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCaa                                                                       <  6:82111/149‑1 (MQ=255)                                                                      aGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCaa                                                                       <  8:97354/149‑1 (MQ=255)                                                                          ggACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                   >  8:335183/1‑149 (MQ=255)                                                                          ggACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                   >  6:362187/1‑149 (MQ=255)                                                                            acacTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTa                                                                 >  4:37443/1‑149 (MQ=255)                                                                              acTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTAGCCCGCCAGGGCGAAAGGCTGACCCGCGACCCGCACAGTGAAGTGGTAAGCATAGAGCCTGACGGGATCAACTACACAGTCAATCGCTacc                                                               >  6:263371/1‑148 (MQ=255)                                                                               cTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                              <  4:397067/149‑1 (MQ=255)                                                                                tGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACAcc                                                             >  3:145169/1‑149 (MQ=255)                                                                                tGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACAcc                                                             >  8:215360/1‑149 (MQ=255)                                                                                     ccGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGc                                                        >  8:81956/1‑149 (MQ=255)                                                                                           tACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCaa                                                  <  5:41276/149‑1 (MQ=255)                                                                                            aCGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGccttcaa                                                  >  1:346858/1‑143 (MQ=255)                                                                                             cGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGc                                                <  4:145169/149‑1 (MQ=255)                                                                                              gAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                  <  2:47074/114‑1 (MQ=255)                                                                                              gAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGca                                               >  3:102691/1‑149 (MQ=255)                                                                                              gAACA‑CCTG‑GCGGCGGACTTCGGCAGCCACGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGca                                               >  1:231599/1‑149 (MQ=255)                                                                                                     cTG‑GGGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACgg                                                                                         <  1:24812/101‑1 (MQ=255)                                                                                                     cTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACgg                                                                                         >  2:24812/1‑101 (MQ=255)                                                                                                      tG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGa                                                                                       <  7:169995/102‑1 (MQ=255)                                                                                                      tG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGa                                                                                       >  8:169995/1‑102 (MQ=255)                                                                                                           ggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                >  7:154830/1‑41 (MQ=255)                                                                                                           ggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                <  8:154830/41‑1 (MQ=255)                                                                                                                       cAGCCACGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                        <  2:231599/149‑1 (MQ=255)                                                                                                                         gCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCa                      <  3:309279/149‑1 (MQ=255)                                                                                                                          ccTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCAc                     <  4:102691/149‑1 (MQ=255)                                                                                                                           cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                         <  7:279294/144‑1 (MQ=255)                                                                                                                           cTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                         >  8:279294/1‑144 (MQ=255)                                                                                                                            tCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                    >  3:229893/1‑148 (MQ=255)                                                                                                                                   ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                          >  5:240440/1‑103 (MQ=255)                                                                                                                                   ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                          <  6:240440/103‑1 (MQ=255)                                                                                                                                    gtgtTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGc           <  5:362187/149‑1 (MQ=255)                                                                                                                                       ttCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCg         >  1:154456/1‑148 (MQ=255)                                                                                                                                          gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTAc                                  <  4:143349/120‑1 (MQ=255)                                                                                                                                          gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTAc                                  >  3:143349/1‑120 (MQ=255)                                                                                                                                           cccTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATgg    <  7:335183/149‑1 (MQ=255)                                                                                                                                            ccTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATgg    >  2:395840/1‑148 (MQ=255)                                                                                                                                             cTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGc   >  2:276976/1‑148 (MQ=255)                                                                                                                                              tCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCa  >  2:133079/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                  CTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACA‑CCTG‑GCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCA  >  NC_002947/418155‑418441

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