Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,170,225	G→A	6.5%	R24C (CGT→TGT)	PP_5408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,170,225	0	G	A	6.5%	180.4 / 2.4	62	R24C (CGT→TGT)	PP_5408	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GTCGGTATGCAAATGCACCAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATGGCAGGAGCCGA  >  NC_002947/6170085‑6170373                                                                                                                                             |                                                                                                                                                     gTCGGTATGCAAATGCACCAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGt                                                                                                                                               <  4:56546/148‑1 (MQ=255)  tCGGTATGCAAATGCACCAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAg                                                                                                                                             <  4:33013/149‑1 (MQ=255)     gTATGCAAATGCACCAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGt                                                                                                                                           >  1:54383/1‑148 (MQ=255)                   cAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCaa                                                                                                                            <  5:105208/149‑1 (MQ=255)                    aaCTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAAc                                                                                                                           <  4:258180/149‑1 (MQ=255)                    aaCTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAAc                                                                                                                           <  2:104630/149‑1 (MQ=255)                       ttGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACAt                                                                                                                         >  4:252527/1‑148 (MQ=255)                          ttAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAg                                                                                                                      <  3:252527/148‑1 (MQ=255)                            aaaGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGaa                                                                                                                    >  4:15602/1‑148 (MQ=255)                            aaaGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGaa                                                                                                                    <  3:15602/148‑1 (MQ=255)                             aaGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAAc                                                                                                                   >  5:185442/1‑148 (MQ=255)                             aaGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAAc                                                                                                                   <  6:185442/148‑1 (MQ=255)                              aGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACg                                                                                                                  >  1:218318/1‑148 (MQ=255)                               gCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGa                                                                                                                <  2:28464/149‑1 (MQ=255)                                ccGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAg                                                                                                               >  2:261011/1‑149 (MQ=255)                                 cGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGc                                                                                                              >  7:144387/1‑149 (MQ=255)                                   tAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGa                                                                                                            >  1:203696/1‑149 (MQ=255)                                      cAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATtc                                                                                                         <  3:245317/149‑1 (MQ=255)                                        gCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATtcat                                                                                                       >  1:72371/1‑149 (MQ=255)                                          cATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGt                                                                                                                                               <  8:3968/107‑1 (MQ=255)                                          cATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGt                                                                                                                                               >  7:3968/1‑107 (MQ=255)                                             gCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGGTTCGTTCAACATAAGAACGGCGCGATTCATCAGGt                                                                                                  >  4:42246/1‑149 (MQ=255)                                                  cccATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGAt                                                                                             >  1:220147/1‑149 (MQ=255)                                                   ccATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGAtt                                                                                            <  1:261011/149‑1 (MQ=255)                                                        aTACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGCTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCaa                                                                                       <  4:62266/149‑1 (MQ=255)                                                            tGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTc                                                                                                                                  >  2:100378/1‑102 (MQ=255)                                                            tGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTc                                                                                                                                  <  1:100378/102‑1 (MQ=255)                                                               cTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATg                                                                                >  6:140149/1‑149 (MQ=255)                                                                         tACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTg                                                                        <  2:72371/147‑1 (MQ=255)                                                                           cctgcctgTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGTATTGTAGGTAAAGGTTTCGTTCAACATAAGAACCGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTc                                                                     <  1:191054/148‑1 (MQ=255)                                                                           cctgcctgTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACAAAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTcc                                                                    <  2:54383/149‑1 (MQ=255)                                                                              gcctgTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTccc                                                                   >  3:179827/1‑147 (MQ=255)                                                                                  gTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACAGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGa                                                                              <  1:94763/132‑1 (MQ=255)                                                                                  gTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACAGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGa                                                                              >  2:94763/1‑132 (MQ=255)                                                                                       tCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAg                                                         >  4:192827/1‑148 (MQ=255)                                                                                         cAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAgg                                                        <  8:144387/147‑1 (MQ=255)                                                                                           gCTGAAGTACCGGAGCTCCGCCGACTTCAGAATATCTCAGTCAGAGTTACTGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAg                                                    <  4:187077/149‑1 (MQ=255)                                                                                            cTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAgg                                                   >  2:114889/1‑149 (MQ=255)                                                                                               aaGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGt                                                  <  5:140149/147‑1 (MQ=255)                                                                                                 gTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACCGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTgagc                                              >  8:145810/1‑149 (MQ=255)                                                                                                      ggAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGaaa                                                     <  7:77502/137‑1 (MQ=255)                                                                                                      ggAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGaaa                                                     >  8:77502/1‑137 (MQ=255)                                                                                                       gAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACAGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGagcgagcg                                         >  4:219689/1‑148 (MQ=255)                                                                                                        aGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGga                                       <  1:114889/149‑1 (MQ=255)                                                                                                          cTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGgaga                                     <  3:192827/149‑1 (MQ=255)                                                                                                          cTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACAGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGgaga                                     <  3:219689/149‑1 (MQ=255)                                                                                                            ccgccgACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACCTAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAAt                                   <  8:55317/149‑1 (MQ=255)                                                                                                                cgACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAgg                               <  4:179827/149‑1 (MQ=255)                                                                                                                    ttCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATATGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGa                            >  6:206249/1‑148 (MQ=255)                                                                                                                    ttCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGa                            >  7:213599/1‑148 (MQ=255)                                                                                                                    ttCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAACAACGGTGAGATTCATCAGGTAGGATTCGCAACAGCATGGATTGCTGGTCCCTAGCGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGa                            >  7:172582/1‑148 (MQ=255)                                                                                                                      cAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAggg                                                              <  7:171688/112‑1 (MQ=255)                                                                                                                      cAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAggg                                                              >  8:171688/1‑112 (MQ=255)                                                                                                                        gAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGaaa                                                     <  7:133871/119‑1 (MQ=255)                                                                                                                        gAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGaaa                                                     <  5:111865/119‑1 (MQ=255)                                                                                                                        gAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGaaa                                                     >  6:111865/1‑119 (MQ=255)                                                                                                                        gAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGc                       >  8:154842/1‑149 (MQ=255)                                                                                                                            atCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTgg                    >  8:26957/1‑148 (MQ=255)                                                                                                                              cTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGAATGGATTCCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGAGGAGAATCACGCGAAAGGCTGGGAc                 >  3:192755/1‑149 (MQ=255)                                                                                                                                     agagTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATgg            >  4:37808/1‑147 (MQ=255)                                                                                                                                     agagTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATgg            >  3:34281/1‑147 (MQ=255)                                                                                                                                        gTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATGGCAGGa       <  2:220147/149‑1 (MQ=255)                                                                                                                                             ggAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATGGCAGGAGCCGa  <  8:213599/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     GTCGGTATGCAAATGCACCAACTTGTTAAAGCCGTAACAGCCATGCCAGCCCATGATACTGCCTACGAGCTTTACCTGCCTGTTGCTCCAGCTGAAGTAACGGAGCTCCGCCGACTTCAGAATATCTAAGTCAGAGTTACGGAATTGTAGGTAAAGGTTTCGTTCAACATAAGAACGGAGCGATTCATCAGGTTGGATTCGCAACAGCATGGATTGCTGGTCCCTAGGGAGCAGGAAAGGTGAGCGAGCGGAGAATCAGGCGAAAGGCTGGGACCATGGCAGGAGCCGA  >  NC_002947/6170085‑6170373

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