Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,378,651	T→A	5.2%	V69D (GTT→GAT) ‡	PP_2086 →	CrfX family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,378,651	0	T	A	5.2%	241.2 / 2.2	76	V69D (GTT→GAT) ‡	PP_2086	CrfX family protein
Reads supporting (aligned to +/- strand):  ref base T (30/42);  new base A (2/2);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

ATGACCGGGACGACCGGGACGACGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGAACTT  >  NC_002947/2378504‑2378791                                                                                                                                                       |                                                                                                                                             aTGACCGGGACGACCGGGACGACGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCACTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGtt                                                                                                                                             >  1:384316/1‑149 (MQ=255)     ccgggacgaccgggacgaCGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTgcgc                                                                                                                                         >  6:326516/1‑149 (MQ=255)     ccgggacgaccgggacgaCGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTgcgc                                                                                                                                         <  3:260047/149‑1 (MQ=255)      cgggacgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGcc                                                                                                                                        <  5:326516/149‑1 (MQ=255)      cgggacgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGcc                                                                                                                                        <  7:351735/149‑1 (MQ=255)        ggacgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCgg                                                                                                                                      >  5:384124/1‑149 (MQ=255)         gacgaccgggacgacGCCGCTTGCCTCGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCAACATGCGCCgag                                                                                                                                     <  2:305724/149‑3 (MQ=255)          acgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGCGCCAGGCGCTGATGAGTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTg                                                                                                                                    >  7:220080/1‑149 (MQ=255)          acgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTg                                                                                                                                    <  6:384124/149‑1 (MQ=255)           cgaccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGc                                                                                                                                   <  1:208373/149‑1 (MQ=255)              ccgggtcgacTCCGCTTGCCGGGGTGGCGTGCTTAAAACCCCCCAACGCCAGGTTGGGGGGGG‑TGATATGTTTATCAAGTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCAACATGCGCCgag                                                                                                                                     <  1:354501/146‑3 (MQ=255)              ccgggacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTc                                                                                                                                <  4:294830/149‑1 (MQ=255)                  gacgacGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGaca                                                                                                                            >  6:256322/1‑149 (MQ=255)                     gacgCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGGATGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACATTGGCTCGGCGCATGTTGCGCCGGTGCGTCGACAct                                                                                                                          >  1:120207/1‑149 (MQ=255)                     gacgCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACActc                                                                                                                         >  6:298013/1‑149 (MQ=255)                        gccgcTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTc                                                                                                                      <  5:269868/149‑1 (MQ=255)                         ccgcTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCCGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGGCGGTGCGTCGACACTCTTcc                                                                                                                     <  8:400670/149‑1 (MQ=255)                          cgcTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTccc                                                                                                                    <  6:88051/149‑1 (MQ=255)                                    gggcgcgTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTccc                                                                                                                    <  7:258183/136‑1 (MQ=255)                                    ggTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTccc                                                                                                                    >  8:258183/1‑139 (MQ=255)                                      tCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGcc                                                                                                         >  5:411285/1‑148 (MQ=255)                                         cgTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATg                                                                                                     <  1:376483/149‑1 (MQ=255)                                            gCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgca                                                                                                  <  7:62288/149‑1 (MQ=255)                                            gCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCACTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgca                                                                                                  <  2:384316/149‑1 (MQ=255)                                             cttaaaaCCGCCACCCGCCAGTTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgca                                                                                                  <  5:278495/145‑1 (MQ=255)                                             cTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgca                                                                                                  <  3:98569/148‑1 (MQ=255)                                               gAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgcagc                                                                                                >  5:275283/1‑148 (MQ=255)                                                aaaaCCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGgcagca                                                                                               >  6:395581/1‑148 (MQ=255)                                                aaaaCCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCaa                                                                                              <  4:336746/149‑1 (MQ=255)                                                  aaccgccaaccgccAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAg                                                                                            <  7:183591/149‑1 (MQ=255)                                                  aaccgccaaccgccAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAg                                                                                            <  5:276703/149‑1 (MQ=255)                                                   accgccaaccgccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGc                                                                                           <  2:84989/149‑1 (MQ=255)                                                      gccaaccgccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAg                                                                                            >  5:380637/1‑145 (MQ=255)                                                       ccaaccgccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGAt                                                                                       <  5:395581/149‑1 (MQ=255)                                                        caaccgccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATa                                                                                      <  8:137444/149‑1 (MQ=255)                                                            cgccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGc                                                                                  >  5:312610/1‑149 (MQ=255)                                                             gccaGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGcc                                                                                 >  3:78619/1‑149 (MQ=255)                                                                            g‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAAc                                                                  >  4:13796/1‑149 (MQ=255)                                                                                aTCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTc                                                                <  3:13796/148‑1 (MQ=255)                                                                                   tGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGAtc                                                            >  4:44552/1‑149 (MQ=255)                                                                                     ttCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCg                                                               <  1:417399/144‑1 (MQ=255)                                                                                     ttCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGAtct                                                           <  5:256322/148‑1 (MQ=255)                                                                                     ttCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGAt                                                             <  2:277596/146‑1 (MQ=255)                                                                                       cAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATgg                                                                     <  1:155706/136‑1 (MQ=255)                                                                                       cAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATgg                                                                     >  2:155706/1‑136 (MQ=255)                                                                                       cAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATctct                                                         <  6:380637/148‑1 (MQ=255)                                                                                             cTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGAcc                                                    <  4:78619/149‑1 (MQ=255)                                                                                                 taggCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGc                                                                                           <  1:335491/104‑1 (MQ=255)                                                                                                 taggCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGc                                                                                           >  2:335491/3‑106 (MQ=255)                                                                                                  tGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCagag                                               >  5:321671/1‑149 (MQ=255)                                                                                                  tGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCagag                                               >  3:420091/1‑149 (MQ=255)                                                                                                      cGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTg                                           <  6:175368/149‑1 (MQ=255)                                                                                                      cGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTg                                           <  1:325923/149‑1 (MQ=255)                                                                                                       gCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCt                                            >  4:147679/1‑147 (MQ=255)                                                                                                       gCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCt                                            <  3:147679/147‑1 (MQ=255)                                                                                                           gAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTc                                        <  3:314085/147‑1 (MQ=255)                                                                                                           gAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCAc                                      >  8:77345/1‑149 (MQ=255)                                                                                                              ccAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTcgggag                                   >  1:156191/1‑144 (MQ=255)                                                                                                               cAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACcgcg                                  <  6:312610/149‑1 (MQ=255)                                                                                                                aGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGa                                 >  7:41598/1‑149 (MQ=255)                                                                                                                      tgctgATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTg                           >  4:210045/1‑149 (MQ=255)                                                                                                                            ttGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCttt                     <  4:420091/149‑1 (MQ=255)                                                                                                                              gCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTg                    <  8:881/148‑1 (MQ=255)                                                                                                                               cggTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGGAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCGAGAGCCTGGTCACCGCGATGACTGCCCTTTGGa                  <  1:403564/147‑1 (MQ=255)                                                                                                                               ccGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGa                  <  6:411285/149‑1 (MQ=255)                                                                                                                               ccGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGa                  >  6:129440/1‑149 (MQ=255)                                                                                                                               ccGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGa                  >  5:83642/1‑149 (MQ=255)                                                                                                                                  tGAACAATGGCTCGTCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCa               <  6:275283/149‑1 (MQ=255)                                                                                                                                     aCAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAg             <  5:298013/148‑1 (MQ=255)                                                                                                                                         tGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGc        >  7:274479/1‑149 (MQ=255)                                                                                                                                          ggCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGc        <  2:120207/148‑1 (MQ=255)                                                                                                                                          gcaccGGCGCAACATGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGc        >  5:65431/5‑148 (MQ=255)                                                                                                                                            cTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGCCCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGaa     >  5:172371/1‑149 (MQ=255)                                                                                                                                            cTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGaa     >  5:165639/1‑149 (MQ=255)                                                                                                                                            accGGCGCAACATGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGCGCCTGGTCCCCGCGATGCCTGCCCTTTGGACCCAGGTAGCGac     >  4:114025/3‑148 (MQ=255)                                                                                                                                               ggCGCATGTTGCTCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTTAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGAAc    <  5:130304/147‑1 (MQ=255)                                                                                                                                                gcgcATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGAACtt  >  4:394366/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                             ATGACCGGGACGACCGGGACGACGCCGCTTGCCTGGGTCGCGTGCTGAAAACCGCCAACCGCCAGGTTGGCGCGGG‑TGATCTGTTCAGC‑‑CTGCTTGGCCGCTGGAGCCAGGCGCTGCTGATTGCCGTGAACAATGGCTCGGCGCATGTTGCGCCGGTGCGTCGACACTCTTCCCGCAACGCTGCCAATGGCAGCAAAGCAGATAAGGCCGATTGAATATGGAACTCGATCTCTGGACCCAGAGCCTGGTCACCGCGATGACTGCCCTTTGGACCAAGGTAGCGAACTT  >  NC_002947/2378504‑2378791

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