Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,883,595	C→G	6.8%	A60G (GCC→GGC)	PP_5157 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,883,595	0	C	G	6.8%	288.1 / 9.0	88	A60G (GCC→GGC)	PP_5157	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/38);  new base G (3/3);  total (47/41)
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.59e-01

TACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCT  >  NC_002947/5883448‑5883742                                                                                                                                                    |                                                                                                                                                     tACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGATAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGcc                                                                                                                                                     >  8:84344/1‑149 (MQ=255)      aTCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGAGAAGTGCTCATTTCCGTACAGGCCGAGAACTCGGCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGAGATCGAATTTGGCCAGGCCCTgt                                                                                                                                                 <  5:216043/148‑2 (MQ=255)        cAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCAGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCAAGGCCCTGGcc                                                                                                                                               >  8:374078/1‑148 (MQ=255)         aGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGcc                                                                                                                                               >  2:411151/1‑147 (MQ=255)         aGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCa                                                                                                                                              <  8:181916/148‑1 (MQ=255)                    tCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGAtt                                                                                                                                 >  2:107103/1‑149 (MQ=255)                    tCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGAtt                                                                                                                                 >  4:151847/1‑149 (MQ=255)                         cGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTgcgc                                                                                                                            <  4:187952/149‑1 (MQ=255)                         cGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTgcgc                                                                                                                            <  4:181696/149‑1 (MQ=255)                             gCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACCCCTTCAAGCAAGACGACCACATGACAGGTTTCGACATCGAAATTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGAGCGcc                                                                                                                         >  3:229665/1‑148 (MQ=255)                              cTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGa                                                                                                                       <  2:205980/149‑1 (MQ=255)                                 aTCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaaattcga                                                                                                                                       <  3:160825/131‑7 (MQ=255)                                 aTCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaaattcga                                                                                                                                       >  4:160825/1‑125 (MQ=255)                                   cgacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGAttt                                                                                                                                >  5:174727/1‑135 (MQ=255)                                   cgacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGAttt                                                                                                                                <  6:174727/135‑1 (MQ=255)                                     acgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                <  8:225094/149‑1 (MQ=255)                                      cgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                >  2:134314/1‑148 (MQ=255)                                            cgccccccGGGGCGAACTGCCCATTCCCGGACAGGCCGAGAACTCGCCATACCCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGccc                                                                                                         <  6:391668/143‑1 (MQ=255)                                             gCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGcccc                                                                                                        >  5:61938/1‑149 (MQ=255)                                                  ccGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTa                                                                                                   >  3:74417/1‑149 (MQ=255)                                                   cGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc                                                                                                  <  7:152322/149‑1 (MQ=255)                                                     gggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTtcgtcg                                                                                                               >  5:350787/1‑134 (MQ=255)                                                     gggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTtcgtcg                                                                                                               <  6:350787/134‑1 (MQ=255)                                                     gggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCg                                                                                                <  4:74417/149‑1 (MQ=255)                                                       gCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              <  1:411151/149‑1 (MQ=255)                                                          aaCTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGccc                                                                                                         >  7:392819/1‑135 (MQ=255)                                                          aaCTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGccc                                                                                                         <  8:392819/135‑1 (MQ=255)                                                             tGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              >  3:266510/1‑143 (MQ=255)                                                             tGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              <  4:266510/143‑1 (MQ=255)                                                               cgcATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGc                                                                                      <  4:295393/149‑1 (MQ=255)                                                                 cATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACTTCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                    <  3:122269/149‑1 (MQ=255)                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGACACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                    <  1:107103/149‑1 (MQ=255)                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCggg                                                                                   >  4:54862/1‑149 (MQ=255)                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCggg                                                                                   >  4:318957/1‑149 (MQ=255)                                                                      ccGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgtgt                                                                                <  7:374078/148‑1 (MQ=255)                                                                      ccGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTc                                                                               <  3:54862/149‑1 (MQ=255)                                                                            aGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGt                                                                         >  8:119661/1‑149 (MQ=255)                                                                             ggCCGAGAACGCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGATTTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                        <  8:137003/149‑1 (MQ=255)                                                                              gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTgg                                                                       <  7:84344/149‑1 (MQ=255)                                                                                  agaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGCAGTGCTGGCGGGTGTCGAGTGTGGCaa                                                                    >  2:106676/1‑148 (MQ=255)                                                                                    aaCTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTa                                                                 <  1:134314/149‑1 (MQ=255)                                                                                       tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCTAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                              >  5:417514/1‑149 (MQ=255)                                                                                       tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              <  8:252508/117‑1 (MQ=255)                                                                                       tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              >  7:252508/1‑117 (MQ=255)                                                                                       tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                              >  7:378262/1‑149 (MQ=255)                                                                                       tCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCACCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                              >  2:344869/1‑149 (MQ=255)                                                                                         gCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                              <  3:151847/147‑1 (MQ=255)                                                                                                gCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt                                                     <  5:101835/149‑1 (MQ=255)                                                                                                  cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt                                                     <  4:383303/147‑1 (MQ=255)                                                                                                  cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt                                                     >  3:383303/1‑147 (MQ=255)                                                                                                  cTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCa                                                   >  3:166903/1‑149 (MQ=255)                                                                                                     caagcaagACGACCACCTGACAGGTTTCGACATCTAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc                                                  <  6:417514/147‑1 (MQ=255)                                                                                                       agcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCg                                              >  2:39582/1‑149 (MQ=255)                                                                                                         caagACGACAACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCACGTc                                            >  2:274437/1‑149 (MQ=255)                                                                                                          aagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                            >  4:361860/1‑148 (MQ=255)                                                                                                          aagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                            >  3:229756/1‑148 (MQ=255)                                                                                                           agACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                            >  8:333080/1‑147 (MQ=255)                                                                                                            gacgacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAg                                         >  3:394927/1‑149 (MQ=255)                                                                                                                  caccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATc                                   >  1:169414/1‑149 (MQ=255)                                                                                                                    ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAt                                                      >  3:74539/1‑128 (MQ=255)                                                                                                                    ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAt                                                      <  4:74539/128‑1 (MQ=255)                                                                                                                    ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTAGTCGAAGCCCCGGATACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACTACATCTCATTCACCCAGTAGAGTTAATCGc                                 >  5:7679/1‑149 (MQ=255)                                                                                                                       gACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              <  3:402235/85‑1 (MQ=255)                                                                                                                       gACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                              >  4:402235/1‑85 (MQ=255)                                                                                                                        aCAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                              <  2:169414/148‑1 (MQ=255)                                                                                                                        aCAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTca                              >  7:12516/1‑148 (MQ=255)                                                                                                                           ggTTTCGACATCGAATTTGGCCAGGCCCTGGGCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAgg                          <  4:6323/149‑1 (MQ=255)                                                                                                                            gTTTCTACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTcaca                            <  8:216494/146‑1 (MQ=255)                                                                                                                            gTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTcaca                            >  7:216494/1‑146 (MQ=255)                                                                                                                            gTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGc                         >  1:171210/1‑149 (MQ=255)                                                                                                                             tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc                                                  <  6:203765/123‑1 (MQ=255)                                                                                                                             tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc                                                  >  5:203765/1‑123 (MQ=255)                                                                                                                                  aCATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGc                   >  1:169870/1‑149 (MQ=255)                                                                                                                                     tCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCg                >  5:64461/1‑149 (MQ=255)                                                                                                                                      gttccttGGCCAGGGCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACg                     <  3:278982/138‑1 (MQ=255)                                                                                                                                      gttccttGGCCAGGGCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACg                     >  4:278982/6‑143 (MQ=255)                                                                                                                                       gAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCTTTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc               <  4:229756/148‑1 (MQ=255)                                                                                                                                       gAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc                                                  >  5:145505/1‑113 (MQ=255)                                                                                                                                       gAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCAc                                                  <  6:145505/113‑1 (MQ=255)                                                                                                                                        aaTTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTc             <  7:333080/149‑1 (MQ=255)                                                                                                                                         aTTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCg            <  1:39582/149‑1 (MQ=255)                                                                                                                                          cttGGCCAGGGCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATc                                                          <  1:372909/101‑1 (MQ=255)                                                                                                                                          cttGGCCAGGGCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATc                                                          >  2:372909/2‑102 (MQ=255)                                                                                                                                             ggCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATAGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg         >  6:333711/1‑148 (MQ=255)                                                                                                                                               ccAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGt        >  6:190565/1‑147 (MQ=255)                                                                                                                                                cAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAg     >  1:362289/1‑149 (MQ=255)                                                                                                                                                 aGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGt                                        >  8:136230/1‑113 (MQ=255)                                                                                                                                                 aGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGt                                        <  7:136230/113‑1 (MQ=255)                                                                                                                                                   gCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCt  >  6:86251/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGA‑ACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCT  >  NC_002947/5883448‑5883742

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