Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,994,761	G→C	11.0%	V496V (GTG→GTC)	PP_5251 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,994,761	0	G	C	11.0%	171.6 / 12.9	64	V496V (GTG→GTC)	PP_5251	transporter
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (5/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

AGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGC  >  NC_002947/5994613‑5994908                                                                                                                                                     |                                                                                                                                                     agcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTg                                                                                                                                                      <  2:37060/149‑1 (MQ=255)   cagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                                    >  8:345274/1‑149 (MQ=255)     ggcggcGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTTGCTGGAGCATCTGCGACATCACCACCGGCTGAccc                                                                                                                                                  >  5:276776/1‑149 (MQ=255)      gcggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑g                                                                                                                                                <  8:186604/149‑1 (MQ=255)      gcggcgCGGGAGCTGGAAGCGATGCGCAAGAGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATGCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTgccgcg                                                                                                                                                 >  2:175608/1‑144 (MQ=255)       cggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gt                                                                                                                                               <  1:212921/149‑1 (MQ=255)        ggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                                    >  8:169169/1‑144 (MQ=255)        ggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                                    <  7:169169/144‑1 (MQ=255)           gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtggt                                                                                                                                            >  3:59591/1‑148 (MQ=255)           gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtggt                                                                                                                                            >  3:359013/1‑148 (MQ=255)             gggAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGGTGGTGat                                                                                                                                         <  7:31809/149‑2 (MQ=255)                aGCTGGAAGCGATTCGACAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATGCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGGTGGTGatgtc                                                                                                                                      <  8:171519/149‑5 (MQ=255)                               gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCAGCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCgg                                                                                                                       <  1:175608/149‑1 (MQ=255)                                 aagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGt                                                                                                                     <  2:177687/149‑1 (MQ=255)                                    cgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAg                                                                                                                  >  1:360435/1‑149 (MQ=255)                                       aaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                                    >  6:271733/1‑113 (MQ=255)                                       aaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGAc                                                                                                                                                    <  5:271733/113‑1 (MQ=255)                                       aaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGAGACAC‑GTGGTGAAATCAAGATCGAGGCGCGGGTCACGAt                                                                                                               >  6:380808/1‑149 (MQ=255)                                        aTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtg                                                                                                              <  2:375151/149‑1 (MQ=255)                                           gcgcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                           <  6:276987/149‑1 (MQ=255)                                              cATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                           >  4:357229/1‑146 (MQ=255)                                              cATTGTCGTGGGGGGTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGATCATCTGCGACATCACCAACTGGTGACCC‑GTGGTGAACTGAAGATCGAGTCGCGGGTCAGGAtgctg                                                                                                           >  5:345876/1‑146 (MQ=255)                                               aTTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTATTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGAGGGTCAGGATGCTGCAGa                                                                                                       >  1:304560/1‑149 (MQ=255)                                                ttGTCGTGGGGGGTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGGCAGGATGCTGCAGAc                                                                                                      <  7:345274/149‑1 (MQ=255)                                                    cGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCa                                                                                                                   <  1:36962/132‑1 (MQ=255)                                                    cGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCa                                                                                                                   >  2:36962/1‑132 (MQ=255)                                                     gTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATg                                                                                                 <  3:379571/149‑1 (MQ=255)                                                            ttGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGa                                                                                          <  6:317613/149‑1 (MQ=255)                                                             tGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                         <  8:244227/149‑1 (MQ=255)                                                               aTGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCGGACCGATGGGGTTGAGGt                                                                                       >  8:287449/1‑149 (MQ=255)                                                               aTGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGt                                                                                       <  2:360435/149‑1 (MQ=255)                                                                  cGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                                     <  7:124531/148‑1 (MQ=255)                                                                  cGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                                     <  4:359013/148‑1 (MQ=255)                                                                  cGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCg                                                                                    <  5:380808/149‑1 (MQ=255)                                                                     gACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGAt                                                                                 <  3:357229/149‑1 (MQ=255)                                                                        tgTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGCGCGGATTCg                                                                              <  3:268239/148‑1 (MQ=255)                                                                          tGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGCGCGGATTCg                                                                              >  4:268239/1‑147 (MQ=255)                                                                           gCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTAGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCa                                                                           >  2:358655/1‑149 (MQ=255)                                                                            cTGATCCGTCGTTGTGGCTGGAGTTCTTCTTCGATGAATCGTGGCTTGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCa                                                                           <  6:276776/148‑1 (MQ=255)                                                                             tGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGAc                                                                                                      <  1:388738/120‑1 (MQ=255)                                                                             tGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGAc                                                                                                      >  2:388738/1‑120 (MQ=255)                                                                                   ggCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTgg                                                                   >  8:328049/1‑149 (MQ=255)                                                                                    gCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTg                                                                                           >  5:82621/1‑124 (MQ=255)                                                                                    gCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTg                                                                                           <  6:82621/124‑1 (MQ=255)                                                                                             cTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACTTCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGTGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTggg                                                           >  2:315846/1‑147 (MQ=255)                                                                                                 aGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCa                                                      >  7:345880/1‑148 (MQ=255)                                                                                                   ttcttcttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGt                                                                                       <  3:110702/113‑1 (MQ=255)                                                                                                   ttcttcttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGt                                                                                       >  4:110702/1‑113 (MQ=255)                                                                                                         ttcGAGCAATCGTGGCTGGAGCATCTGCGACATCCCCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCTTAAGGcaca                                             >  6:319603/1‑149 (MQ=255)                                                                                                          tcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCacac                                            >  5:355241/1‑149 (MQ=255)                                                                                                             aGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCAt                                         >  7:175363/1‑149 (MQ=255)                                                                                                             aGGAATCGTGGCTGGAGCATCAGCGACATCACCACCGGGTGACCC‑GAGGTGAACTGAAGATCGAGGCGCGGGTCACGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCAt                                         >  1:326663/1‑149 (MQ=255)                                                                                                                aaTCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGa                                      <  5:347490/149‑1 (MQ=255)                                                                                                                    gTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGACGATCGTGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGGGCGGATTCGGCATCTATTGGCTGGTGGGCAGCATAAGGCACACCATTGATGc                                   <  5:319603/148‑1 (MQ=255)                                                                                                                    gTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGcc                                  <  6:355241/149‑1 (MQ=255)                                                                                                                       gCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCgg                                <  3:137374/148‑1 (MQ=255)                                                                                                                           gatcATCTGCGACATCACCACCGGTTGACCC‑GTGGTGAACTGAAGATCGAGGAGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGTCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGccc                             <  6:345876/144‑1 (MQ=255)                                                                                                                                      agttCACCA‑CGGGTCACCC‑‑‑GGTGAACTGAAGATCCAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgcg             >  2:258288/4‑149 (MQ=255)                                                                                                                                      aCATCACCACCGGGTGACGC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTc                <  2:304560/149‑1 (MQ=255)                                                                                                                                       gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGTGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcg               >  8:342514/3‑149 (MQ=255)                                                                                                                                        ttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgc              >  2:189559/2‑149 (MQ=255)                                                                                                                                        ttCACCA‑CGGGTCACCCGGTGGTGAAATGAAGATCGAGGCGAGGGACAGGATGATGCAGACCGATGGGGGTGCGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACAATTTATGCCGGCCCTTACGGCCCATTagc              >  2:361352/2‑149 (MQ=255)                                                                                                                                          caccacCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgc              >  5:39757/1‑147 (MQ=255)                                                                                                                                            accacGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgcgcg           >  4:230427/5‑147 (MQ=255)                                                                                                                                              accGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGACGGCCATTAAGGCCAATTCGCGGGTg         >  7:88913/1‑148 (MQ=255)                                                                                                                                                    tGAGCC‑GTGGTTTACTGATGATCGAGGCGCGGGTCAGGATGCTGCAGTCCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCTCACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGc  <  2:326663/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCGC  >  NC_002947/5994613‑5994908

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