Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	417,375	A→G	7.6%	Q362R (CAG→CGG)	waaG →	lipopolysaccharide glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	417,375	0	A	G	7.6%	192.0 / 5.1	66	Q362R (CAG→CGG)	waaG	lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base G (3/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 = 7.80e-01

GACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417227‑417519                                                                                                                                                     |                                                                                                                                                 gACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgca                                                                                                                                                  <  3:365282/149‑1 (MQ=255)      cGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtg                                                                                                                                             >  6:225886/1‑149 (MQ=255)      cGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtg                                                                                                                                             >  8:166502/1‑149 (MQ=255)       ggCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtgc                                                                                                                                            <  3:261395/149‑1 (MQ=255)         ccTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctg                                                                                                                                          <  3:12959/149‑1 (MQ=255)             gTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTCCGGCAtgctgtaca                                                                                                                                      <  7:319373/149‑5 (MQ=255)             gTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGcc                                                                                                                                        <  7:126004/147‑1 (MQ=255)               gCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACg                                                                                                                                    >  6:35048/1‑149 (MQ=255)                atGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGt                                                                                                                                   <  1:354075/148‑1 (MQ=255)                         ccGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgcagc                                                                                                                                                >  1:164222/1‑127 (MQ=255)                         ccGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgcagc                                                                                                                                                <  2:164222/127‑1 (MQ=255)                                      agcagcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCa                                                                                                             <  7:166502/149‑1 (MQ=255)                                         agcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATg                                                                                                           >  7:333236/1‑148 (MQ=255)                                          gcTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGaa                                                                                                         <  2:226152/149‑1 (MQ=255)                                             cAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAACACCGCATGAAGCt                                                                                                      >  8:87510/1‑149 (MQ=255)                                              aaCGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGc                                                                                                       <  3:39369/147‑1 (MQ=255)                                                    tATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTg                                                                                               >  8:148997/1‑149 (MQ=255)                                                        tGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGcc                                                                                            >  5:148460/1‑148 (MQ=255)                                                        tGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGcc                                                                                            >  8:42393/1‑148 (MQ=255)                                                        tGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGcc                                                                                            >  3:320607/1‑148 (MQ=255)                                                           aGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAc                                                                                        <  3:209146/149‑1 (MQ=255)                                                                aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTc                                                                                   <  7:42393/149‑1 (MQ=255)                                                                aTGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTc                                                                                   <  6:370300/149‑1 (MQ=255)                                                                      gAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGTATGAAGCTGATATTGGCCGAACCCTTCAAAc                                                                               <  4:323856/147‑1 (MQ=255)                                                                         gACCCGCAGGCGCGTGCCAGCAGGTAGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGAAGCATGCAGCCGACGTGATCCTGGCGCAAAAAACCGCATGAAGCTGATATTGTCCGAACCCTACAAACGCCtg                                                                          >  7:164003/1‑149 (MQ=255)                                                                          aCCAGCTGGCACGCGCCTGCGGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCtgt                                                                         >  6:373335/1‑149 (MQ=255)                                                                           cccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCtgtg                                                                        >  4:44381/1‑149 (MQ=255)                                                                            ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTgg                                                                       >  2:276597/1‑149 (MQ=255)                                                                            ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTgg                                                                       >  2:337385/1‑149 (MQ=255)                                                                                  gcgcGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCggg                                                                 >  4:292751/1‑149 (MQ=255)                                                                                         ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGAt                                                           <  3:292751/148‑1 (MQ=255)                                                                                         ccAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATg                                                          <  1:337385/149‑1 (MQ=255)                                                                                          cAGCTGGGCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGc                                                         <  8:16494/149‑1 (MQ=255)                                                                                           aGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGcc                                                        >  6:136437/1‑149 (MQ=255)                                                                                            gCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAACATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCg                                                                                                               >  6:134493/1‑93 (MQ=255)                                                                                            gCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAACATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCg                                                                                                               <  5:134493/93‑1 (MQ=255)                                                                                                gTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCGTTCAAACGCCTGTGGGCCGGGCGTGATGCCTACgc                                                   >  8:241378/1‑148 (MQ=255)                                                                                                gTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCg                                                    <  6:148460/148‑1 (MQ=255)                                                                                                    cgcAAGGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGcc                                               <  7:148997/149‑1 (MQ=255)                                                                                                      cAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGt                                             <  5:136437/149‑1 (MQ=255)                                                                                                           gCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        >  5:327025/1‑149 (MQ=255)                                                                                                           gCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        <  3:298164/149‑1 (MQ=255)                                                                                                           gCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAg                                         <  1:191049/148‑1 (MQ=255)                                                                                                             cTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGgcg                                      >  3:330519/1‑149 (MQ=255)                                                                                                              tGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGgcgc                                     >  4:308284/1‑149 (MQ=255)                                                                                                              tGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGgcgc                                     >  5:26171/1‑149 (MQ=255)                                                                                                                  aTTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGc                                  >  5:357355/1‑148 (MQ=255)                                                                                                                      gCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACGCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTTGAGGCGCTGCAAGGc                             >  6:117557/1‑149 (MQ=255)                                                                                                                      gCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGc                             >  7:130668/1‑149 (MQ=255)                                                                                                                       cTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCg                            <  3:44381/149‑1 (MQ=255)                                                                                                                       cTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTCGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGc                             >  6:334837/1‑148 (MQ=255)                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        >  2:182180/1‑136 (MQ=255)                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAgg                                        <  1:182180/136‑1 (MQ=255)                                                                                                                        tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGc                             >  3:129048/1‑147 (MQ=255)                                                                                                                          aTAGCGACGATCTGTACACAATGCCGCAGCATGCTGCCGACGTTATGCTGGGGCAATAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGACTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                          <  8:164003/148‑1 (MQ=255)                                                                                                                              cgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTa                     >  2:323149/1‑149 (MQ=255)                                                                                                                                   atgtGTACAGCATGCGGCAGCATGCTGCCTACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgcg                <  6:357355/146‑1 (MQ=255)                                                                                                                                       gTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTCGGGCAAAAAAACGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCt            >  1:82820/1‑149 (MQ=255)                                                                                                                                        tccAGGATGCCGCAGCATGGTGCCGACTTGATCGTGGGGCTAAAAACCGCTTGAAGCTGATATTGGCCGAACCGTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGa               <  4:209285/143‑1 (MQ=255)                                                                                                                                        tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgc                 <  4:330519/143‑1 (MQ=255)                                                                                                                                         gcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAACCCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTg           >  4:221883/2‑148 (MQ=255)                                                                                                                                         gcAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCtt                                                      >  7:299949/2‑105 (MQ=255)                                                                                                                                          cctcATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTgg          <  5:225886/145‑1 (MQ=255)                                                                                                                                          cAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTgg          >  2:287512/1‑148 (MQ=255)                                                                                                                                           aGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTgg          >  2:374149/1‑147 (MQ=255)                                                                                                                                           aGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTgg          >  8:227752/1‑147 (MQ=255)                                                                                                                                           aGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGaa        <  6:26171/149‑1 (MQ=255)                                                                                                                                           aGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGaa        <  1:374149/149‑1 (MQ=255)                                                                                                                                                 ccGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGgcgc  >  3:64239/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCGC  >  NC_002947/417227‑417519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑