Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,720,490	G→A	6.0%	V199V (GTG→GTA)	PP_4177 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,720,490	0	G	A	6.0%	200.6 / 2.4	66	V199V (GTG→GTA)	PP_4177	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTGGTACTGATTGGCGGCGCAC  >  NC_002947/4720343‑4720635                                                                                                                                                    |                                                                                                                                                  cGTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATATCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCAtgtg                                                                                                                                                   >  8:6070/1‑148 (MQ=255)  gTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGTGACCACCTTGCTGGTACCGACCGAGCAtgtg                                                                                                                                                   >  2:161677/1‑147 (MQ=255)  gTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCa                                                                                                                                                 <  3:306676/149‑1 (MQ=255)  gTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATATCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCa                                                                                                                                                 <  7:6070/149‑1 (MQ=255)     aGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCagc                                                                                                                                               >  1:233956/1‑148 (MQ=255)     aGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCagc                                                                                                                                               <  1:152771/148‑1 (MQ=255)           aGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCTCCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCt                                                                                                                                        <  1:301606/149‑1 (MQ=255)           aGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCt                                                                                                                                        <  5:38139/149‑1 (MQ=255)             ccggggTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCGTGGCGTCACCGGTGCCGGTGCTGTCGCGGTCCCTTTATGCGCGTGGCCGGAGGCCACCGGCGCCCCCCTTGCTGGTGCCGCCCGGGCGTGGGCAGCGCGCGCTgg                                                                                                                                      >  8:283286/4‑148 (MQ=255)               tGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTgtac                                                                                                                                    >  7:336818/1‑149 (MQ=255)                   tgtgCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTgtacgtac                                                                                                                                <  5:77131/149‑1 (MQ=255)                         tATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATg                                                                                                                          <  4:64393/149‑1 (MQ=255)                          aTCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATgg                                                                                                                         >  5:99180/1‑149 (MQ=255)                            cTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGt                                                                                                                        >  4:92477/1‑148 (MQ=255)                             tGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTACAGCGAGCGCTGTACGTACTCGATGGt                                                                                                                        >  6:247067/1‑147 (MQ=255)                               cATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACTGGCGACCACCTTGCTGGTACCGACCGAGCATGTACAGCGCTCGCTgcacatgct                                                                                                                               <  6:138885/138‑6 (MQ=255)                               cATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACTGGCGACCACCTTGCTGGTACCGACCGAGCATGTACAGCGCTCGCTgcacatgct                                                                                                                               >  5:138885/1‑133 (MQ=255)                               cATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTg                                                                                                                       >  8:42087/1‑146 (MQ=255)                               cATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTg                                                                                                                       <  7:42087/146‑1 (MQ=255)                               cATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGact                                                                                                                    <  6:99180/149‑3 (MQ=255)                                aTGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGTGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGtt                                                                                                                   <  1:161677/149‑1 (MQ=255)                                  gATTTCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCACTACCCTTTATGCGCATTTCAGGATGCAACCGGCGACCACCTTTCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTAGATGGTCAGTgga                                                                                                                 >  4:186913/1‑149 (MQ=255)                                   aTTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTtgat                                                                                                                >  6:378252/1‑149 (MQ=255)                                       cTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTtgatgat                                                                                                             <  2:233956/148‑1 (MQ=255)                                           cAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACg                                                                                                        <  4:63787/149‑1 (MQ=255)                                            aGCGGGTTCTGCCTGGTGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGa                                                                                                       >  7:284202/1‑149 (MQ=255)                                            aGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGa                                                                                                       >  8:315530/1‑149 (MQ=255)                                                        cTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTgtac                                                                                                                                    <  3:355130/108‑1 (MQ=255)                                                        cTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTgtac                                                                                                                                    >  4:355130/1‑108 (MQ=255)                                                         tGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCCTGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAg                                                                                          >  2:319592/1‑149 (MQ=255)                                                             gTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGt                                                                                      >  2:38149/1‑149 (MQ=255)                                                              tCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGAt                                                                                                      >  7:251233/1‑132 (MQ=255)                                                              tCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGAt                                                                                                      <  8:251233/132‑1 (MQ=255)                                                                  cggtgccggtgcTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGc                                                                                 <  8:284202/149‑1 (MQ=255)                                                                    gtgccggtgcTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCt                                                                               >  7:36614/1‑149 (MQ=255)                                                                      gccggtgcTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCtgg                                                                             <  2:165486/149‑1 (MQ=255)                                                                      gccggtgcTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGACtgg                                                                             >  5:288577/1‑149 (MQ=255)                                                                       ccggtgcTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTTATGATGAACGATGAAGACGTGGAGCCGTGCAGCCtggt                                                                            <  3:92477/149‑1 (MQ=255)                                                                                  cccTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACATGGAGCCGTGCAGCCTGGTGGTACTGCCg                                                                  >  5:14168/1‑148 (MQ=255)                                                                                   ccTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCgg                                                                 >  4:53683/1‑148 (MQ=255)                                                                                    cTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGaa                                                               <  3:53683/149‑1 (MQ=255)                                                                                     tACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAg                                                              >  2:357875/1‑149 (MQ=255)                                                                                      aCCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGTATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAgg                                                             >  4:22806/1‑149 (MQ=255)                                                                                      aCCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTACAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAgg                                                             <  5:247067/149‑1 (MQ=255)                                                                                       cccTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGt                                                                                                                        >  1:151412/1‑89 (MQ=255)                                                                                       cccTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGt                                                                                                                        <  2:151412/89‑1 (MQ=255)                                                                                          tttATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGa                                                                                                                      <  7:219673/88‑1 (MQ=255)                                                                                          tttATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGa                                                                                                                      >  8:219673/1‑88 (MQ=255)                                                                                             aTGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAaga                                                       <  5:378252/148‑1 (MQ=255)                                                                                              tGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAagag                                                      <  8:36614/148‑1 (MQ=255)                                                                                                cgcATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACATGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGt                                                    <  6:14168/148‑1 (MQ=255)                                                                                                  cATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGAc                                                 <  7:315530/149‑1 (MQ=255)                                                                                                       ccGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTg                                             >  2:173813/1‑148 (MQ=255)                                                                                                            tGCAACCGGCGAGCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCTATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAaga                                                       <  3:211135/133‑1 (MQ=255)                                                                                                            tGCAACCGGCGACCACCTTGCTGGTACCTACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTCAGTTGATGATTAACGATGAAGACGTGCAGCCGTGCAGCAAGGTGGTACTGCCGGAAGGGGAaga                                                       >  4:211135/1‑133 (MQ=255)                                                                                                               aaCCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGATCGCTGTACGTAATCGATGGTGAGTTGATTATGAACGATGAAGACGTGTAGCAGTGCAGCCTGGTGGTACTTCATGAAGGGGaa                                                         <  6:37751/128‑1 (MQ=255)                                                                                                               aaCCGGCGACCACCTTGATGGTACCGACCGATCATGTGCATAGAGCGATGTACGTAATCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGACTGGTGGTACTGCCGGAAGTGGaa                                                         >  5:37751/1‑128 (MQ=255)                                                                                                               aaCCAGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAg                                    >  3:176029/1‑149 (MQ=255)                                                                                                                  cGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGt                                 <  1:38149/149‑1 (MQ=255)                                                                                                                     cGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGaa                              <  8:336818/149‑1 (MQ=255)                                                                                                                      gACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAAt                             >  5:348956/1‑149 (MQ=255)                                                                                                                        ccaccTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGc                           >  1:125334/1‑149 (MQ=255)                                                                                                                         caccTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTATGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGTGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGc                           >  4:328118/1‑148 (MQ=255)                                                                                                                                cTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTgg                    <  1:357875/148‑1 (MQ=255)                                                                                                                                cTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTgg                    <  1:173813/148‑1 (MQ=255)                                                                                                                                 tGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTGGt                   >  8:104013/1‑148 (MQ=255)                                                                                                                                   gTACCGACCGAGCATGTGCAGCGCGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGACGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGCGGGGCCGTTGTATGCGGAGGGTGAATGCCAATTGGTAcg                >  5:289709/1‑148 (MQ=255)                                                                                                                                               cATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATCAGTACCAATTggcattca    <  4:176029/149‑6 (MQ=255)                                                                                                                                                 tgtgCAGCGAGCGCCGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTGGTACTGATTGGCGGCGCAc  >  8:146121/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGTCAGTGACAGCCTGGGTGTGCGTATCTGCATGATTGCTGGCAGCGGGTTCTGCCTGGCGTCACCGGTGCCGGTGCTGTCCCCTACCCTTTATGCGCATGTCCGGATGCAACCGGCGACCACCTTGCTGGTACCGACCGAGCATGTGCAGCGAGCGCTGTACGTACTCGATGGTGAGTTGATGATGAACGATGAAGACGTGGAGCCGTGCAGCCTGGTGGTACTGCCGGAAGGGGAAGAGGTGACGTTGTATGCGGAGGGTGAATGCCAATTGGTACTGATTGGCGGCGCAC  >  NC_002947/4720343‑4720635

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