Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,994,755	Δ1 bp	9.7%	coding (1482/1599 nt)	PP_5251 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,994,755	0	C	.	9.7%	181.2 / 14.1	62	coding (1482/1599 nt)	PP_5251	transporter
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base . (4/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTG  >  NC_002947/5994608‑5994901                                                                                                                                                    |                                                                                                                                                    cTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCaccacc                                                                                                                                                    <  3:270012/149‑1 (MQ=255)      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                                                                              <  4:359013/148‑1 (MQ=255)                                                                       cGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                              <  7:124531/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                                                                                >  4:110702/1‑113 (MQ=255)                                                                                                        ttcttcttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGt                                                                                <  3:110702/113‑1 (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)                                                                                                                                                   accGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGACGGCCATTAAGGCCAATTCGCGGGTg  >  7:88913/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTG  >  NC_002947/5994608‑5994901

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