Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,255,974	G→T	9.8%	intergenic (+34/‑670)	asd → / → PP_1990	aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,255,974	0	G	T	9.8%	241.1 / 13.9	82	intergenic (+34/‑670)	asd/PP_1990	aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base G (41/33);  new base T (4/4);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-01

CGCAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGC  >  NC_002947/2255827‑2256113                                                                                                                                                    |                                                                                                                                                    cgcAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGt                                                                                                                                                    <  4:221601/149‑1 (MQ=255)   catgCTGAAGATGGGGTCGCAGTAGCTGGGCGCGTTCACGGTGGGTGACGAGCTGCTGTGGGTCGCCGCCGATCCGCTGCGCCGCATGTTGCGGTTTCTGCTGGAGCTTTGATCGTTGAGGCTTCACCTGAAAAATACCGCACTGGTGa                                                                                                                                                  <  3:394465/146‑1 (MQ=255)     aGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGa                                                                                                                                                  >  1:34286/1‑147 (MQ=255)      gCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑g                                                                                                                                              >  3:129497/1‑149 (MQ=255)        tGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTg                                                                                                                                            >  2:166268/1‑149 (MQ=255)          aaCATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCg                                                                                                                                          >  5:118151/1‑149 (MQ=255)          aaCATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCg                                                                                                                                          <  2:155601/149‑1 (MQ=255)           aCATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCgg                                                                                                                                         >  1:77505/1‑149 (MQ=255)             aTGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTa                                                                                                                                       <  1:166268/149‑1 (MQ=255)               ggggTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAtt                                                                                                                                     >  4:168271/1‑149 (MQ=255)               ggggTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAtt                                                                                                                                     >  1:180339/1‑149 (MQ=255)                  gTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAttttt                                                                                                                                  >  2:57253/1‑149 (MQ=255)                   tCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAttttt‑‑‑‑t                                                                                                                             <  5:380127/149‑1 (MQ=255)                        gTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGc                                                                                                                        <  2:180339/149‑1 (MQ=255)                          aCCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATgat                                                                                                                       >  8:163332/1‑146 (MQ=255)                           ccTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCt                                                                                                                      >  2:342452/1‑148 (MQ=255)                           ccTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCCTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCt                                                                                                                      <  7:340665/148‑1 (MQ=255)                            cTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAttttt‑‑‑‑tt                                                                                                                            >  7:210192/1‑141 (MQ=255)                            cTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTAttttt‑‑‑‑tt                                                                                                                            <  8:210192/141‑1 (MQ=255)                              gggCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGaaa                                                                                                                  <  3:168271/149‑1 (MQ=255)                                gcgcGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAg                                                                                                                 >  3:101984/1‑148 (MQ=255)                                    gTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGAACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATg                                                                                                             >  5:362523/1‑148 (MQ=255)                                        aCCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCAttt                                                                                                        >  7:7967/1‑149 (MQ=255)                                          cGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑gg                                                                                                     <  6:118151/149‑1 (MQ=255)                                                    cAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑cc                                                                                         >  2:309134/1‑149 (MQ=255)                                                    cAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑cc                                                                                         >  4:19281/1‑149 (MQ=255)                                                     aGCTTCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCt                                                                                        <  2:268160/149‑1 (MQ=255)                                                        tgctgTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAAt                                                                                     >  4:96979/1‑149 (MQ=255)                                                          ctgTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCg                                                                                   >  2:49496/1‑149 (MQ=255)                                                          ctgTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCg                                                                                   <  8:7967/149‑1 (MQ=255)                                                             tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGTCACCA‑GTGCGGTAtttttca                                                                                                                                <  4:218561/108‑3 (MQ=38)                                                             tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGTCACCA‑GTGCGGTAtttttca                                                                                                                                >  3:218561/1‑106 (MQ=38)                                                             tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCg                                                                                <  1:57253/149‑1 (MQ=255)                                                                ggCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAAt                                                                                     <  5:417201/141‑1 (MQ=255)                                                                ggCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAAt                                                                                     >  6:417201/1‑141 (MQ=255)                                                                 gcgcCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGGTACATAATCGCCgta                                                                              >  1:343356/1‑147 (MQ=255)                                                                  cgccgcCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATaa                                                                           >  8:63743/1‑149 (MQ=255)                                                                   gccgccGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATTCCTGAAAGTATGCATTT‑GGTTCTCGTG‑‑CCTAATCTCCGGATAAt                                                                          >  8:130187/1‑149 (MQ=255)                                                                    ccgccgAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACACCTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAAt                                                                          <  1:115707/149‑1 (MQ=255)                                                                       ccgAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCAt                                                                      >  2:145785/1‑149 (MQ=255)                                                                       ccgAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCAt                                                                      >  3:128122/1‑149 (MQ=255)                                                                        cgAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATg                                                                     >  2:286955/1‑149 (MQ=255)                                                                         gAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGc                                                                    >  6:85676/1‑149 (MQ=255)                                                                             cGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCt                                                                >  7:203367/1‑149 (MQ=255)                                                                               cTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGt                                                              >  8:165474/1‑149 (MQ=255)                                                                                 gcgcCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGttt                                                            <  3:218650/149‑1 (MQ=255)                                                                                    ccgcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAAt                                                         >  8:306477/1‑149 (MQ=255)                                                                                    ccgcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAAt                                                         >  7:358817/1‑149 (MQ=255)                                                                                    ccgcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGAGGG‑‑CCTAATCTCCGGATAATTAATTCCGCTGTTTAAt                                                         >  8:312146/1‑149 (MQ=255)                                                                                          gTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGa                                                    >  7:1212/1‑148 (MQ=255)                                                                                           ttGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGaaa                                                  >  5:372790/1‑149 (MQ=255)                                                                                             gCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAAt                                                 <  6:372790/148‑1 (MQ=255)                                                                                             gCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAAt                                                 <  4:129497/148‑1 (MQ=255)                                                                                              cGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGCG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGt                                               <  1:286955/149‑1 (MQ=255)                                                                                               ggATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTTAGGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGt                                               <  1:309134/149‑1 (MQ=255)                                                                                               ggATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTc                                              <  2:34286/149‑1 (MQ=255)                                                                                                     ggctgGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAAtt                                                                         >  8:155570/2‑116 (MQ=255)                                                                                                     ggctgGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAAtt                                                                         <  7:155570/115‑1 (MQ=255)                                                                                                       ctgGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGctct                                        <  1:342452/147‑1 (MQ=255)                                                                                                          gAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGAc                                     <  2:77505/147‑1 (MQ=255)                                                                                                                     gTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgcc                         <  3:19281/148‑1 (MQ=255)                                                                                                                      ttCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgccgt                       >  6:360891/1‑149 (MQ=255)                                                                                                                            cTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgccgttgccgt                 <  8:358817/149‑1 (MQ=255)                                                                                                                             ttCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCAGGATAATTAATGCCGCTGTTTAATGGCGAAATGTCTGCTCTCACAGAgccgttgccgttgccgtt                >  4:216651/1‑149 (MQ=255)                                                                                                                                 ccTGAAAAATACCGCACGGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCGGCTGTTTAATGGCGAAATGTCTGCTCTGACAGCGCCGTTGCCGTTGCCGttgctt            >  5:221257/1‑149 (MQ=255)                                                                                                                                   tGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctttt          >  6:150506/1‑149 (MQ=255)                                                                                                                                    gAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctttt          <  1:145785/148‑1 (MQ=255)                                                                                                                                     aaaaaTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgc        <  5:360891/149‑1 (MQ=255)                                                                                                                                     aaaaaTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgc        >  5:34588/1‑149 (MQ=255)                                                                                                                                     aaaaaTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgc        >  4:349359/1‑149 (MQ=255)                                                                                                                                        aaTACCGCACTGTCACCA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATg                                                        >  7:175848/1‑98 (MQ=38)                                                                                                                                        aaTACCGCACTGTCACCA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATg                                                        <  8:175848/98‑1 (MQ=38)                                                                                                                                         aTACCGCACTGTCACCA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCt                                             >  7:177963/1‑108 (MQ=39)                                                                                                                                         aTACCGCACTGTCACCA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCt                                             <  8:177963/108‑1 (MQ=39)                                                                                                                                         aTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgctttt    <  8:1212/149‑1 (MQ=255)                                                                                                                                         aTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgctttt    <  8:203367/149‑1 (MQ=255)                                                                                                                                         aTACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgctttt    <  1:49496/149‑1 (MQ=255)                                                                                                                                          tACCGCACTGTCACCA‑GTGCGGTATTTTTCAGGTGAAGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctttt          >  7:96728/1‑146 (MQ=25)                                                                                                                                          tACCGCACTGTCACCA‑GTGCGGTATTTTTCAGGTGAAGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctttt          <  8:96728/146‑1 (MQ=25)                                                                                                                                           aCCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAAt                                                         >  2:189393/1‑94 (MQ=255)                                                                                                                                           aCCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAAt                                                         <  1:189393/94‑1 (MQ=255)                                                                                                                                            ccGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgcttttgc  <  7:63743/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGCAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACA‑GTGCGGTATTTTT‑‑‑‑TTATGCCTGAAAGTATGCATTT‑GGGTCGCGGG‑‑CCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGC  >  NC_002947/2255827‑2256113

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