Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,551,117	C→T	5.7%	Q76* (CAG→TAG) ‡	PP_3135 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,551,117	0	C	T	5.7%	257.7 / 4.4	88	Q76* (CAG→TAG) ‡	PP_3135	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base T (3/2);  total (44/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 = 4.54e-01

CGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGG  >  NC_002947/3550970‑3551265                                                                                                                                                    |                                                                                                                                                     cgcCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCGCTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAg                                                                                                                                                    <  1:1044502/150‑1 (MQ=255)        ccGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             <  1:1182932/150‑1 (MQ=255)         cGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgc                                                                                                                                            >  1:144870/1‑150 (MQ=255)         cGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgc                                                                                                                                            >  2:124473/1‑150 (MQ=255)           aCAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgcgc                                                                                                                                          <  2:867045/150‑1 (MQ=255)             agcgagcgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgcgcgc                                                                                                                                        >  1:1696227/1‑148 (MQ=255)             agcgagcgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCAc                                                                                                                                        <  2:571553/150‑1 (MQ=255)                 agcgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCa                                                                                                                                    <  1:792879/150‑1 (MQ=255)                         gCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAg                                                                                                                            <  1:1595754/150‑1 (MQ=255)                         gCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCAATGCTGCTGAGCAACAAGCTGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAg                                                                                                                            <  1:1672954/150‑1 (MQ=255)                          ccAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGt                                                                                                                           <  1:612721/150‑1 (MQ=255)                          ccAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGt                                                                                                                           >  1:1144809/1‑150 (MQ=255)                           cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCTGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                          >  2:568649/1‑150 (MQ=255)                           cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCCCTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                          >  2:679454/1‑150 (MQ=255)                            aGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTAt                                                                                                                         >  2:1516327/1‑150 (MQ=255)                              ttACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGc                                                                                                                       <  1:334429/150‑1 (MQ=255)                                aCCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCt                                                                                                                     >  2:1489634/1‑150 (MQ=255)                                aCCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCt                                                                                                                     >  2:1489605/1‑150 (MQ=255)                                 ccTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCtt                                                                                                                    <  1:765914/150‑1 (MQ=255)                                  cTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCGGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTtg                                                                                                                   >  2:911592/1‑150 (MQ=255)                                  cTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTtg                                                                                                                   >  1:1782689/1‑150 (MQ=255)                                   tggtggTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTtgc                                                                                                                  >  1:1114894/1‑150 (MQ=255)                                      tggtCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTtgctgc                                                                                                               <  2:1782689/150‑1 (MQ=255)                                         tCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             <  2:542358/117‑1 (MQ=255)                                         tCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                             >  1:542358/1‑117 (MQ=255)                                              aaTGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATAcc                                                                                                       >  2:52633/1‑150 (MQ=255)                                                 gCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGAt                                                                                                    <  1:1047898/150‑1 (MQ=255)                                                  ccTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATg                                                                                                   <  2:1770874/150‑1 (MQ=255)                                                   cTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGc                                                                                                  >  1:341192/1‑150 (MQ=255)                                                       cgcgATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTc                                                                                              >  1:1302627/1‑150 (MQ=255)                                                          gATGACTCGTTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCCTGCTGCTCAATACCGATGCCTTCGtt                                                                                           <  2:173186/150‑1 (MQ=255)                                                          gATGACTCGGTGCAGCGAATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGtt                                                                                           <  2:374230/150‑1 (MQ=255)                                                           aTGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTg                                                                                          >  2:850845/1‑150 (MQ=255)                                                                tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGt                                                                                            <  1:1919553/143‑1 (MQ=255)                                                                tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGt                                                                                            >  2:1919553/1‑143 (MQ=255)                                                                tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                     >  2:722734/1‑150 (MQ=255)                                                                      cAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCg                                                                               <  1:1300548/150‑1 (MQ=255)                                                                         cgcATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTg                                                                            <  1:439907/150‑1 (MQ=255)                                                                           cATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCa                                                                          >  2:1439690/1‑150 (MQ=255)                                                                             tGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCaaa                                                                        <  1:744312/150‑1 (MQ=255)                                                                                cccAGGGGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATCCCGATCCCTGCCTTGCCGCCGACTCGCTGCaaaaaa                                                                     >  2:16563/1‑150 (MQ=255)                                                                                 ccAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAAc                                                                    <  1:911592/150‑1 (MQ=255)                                                                                   aGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGc                                                                  <  1:184911/150‑1 (MQ=255)                                                                                    ggTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 <  1:1004283/150‑1 (MQ=255)                                                                                          ccccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTa                                                           >  2:1059198/1‑150 (MQ=255)                                                                                           cccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTAc                                                          >  2:1674677/1‑150 (MQ=255)                                                                                           cccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTAc                                                          >  1:1243199/1‑150 (MQ=255)                                                                                            ccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACa                                                         <  1:903209/150‑1 (MQ=255)                                                                                             cGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCAGACTCGCTGCAAAAAACGCTCGACTACAt                                                        >  1:1815949/1‑150 (MQ=255)                                                                                                 gCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGCCTACATGGAc                                                    >  1:428230/1‑150 (MQ=255)                                                                                                 gCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCACTTCAAAAAACGCTCGACTACATGGAc                                                    <  1:1697486/150‑1 (MQ=255)                                                                                                   aCTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGCCTCGCTGCCAAAAACGCTCGACTACATGGACGc                                                  >  2:800314/1‑150 (MQ=255)                                                                                                    ctgctgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGcc                                                 >  1:1878792/1‑150 (MQ=255)                                                                                                     tgctgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGccc                                                >  1:614407/1‑150 (MQ=255)                                                                                                     tgctgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGccc                                                >  1:1605920/1‑150 (MQ=255)                                                                                                     tgctgctgAGCAACAAGCAGAACGAGGGCTCCGGCCGTTCCAACAAACAGTTGTTGGCGCCGCTCACGGGCGAGTAGGCCGCTCTTCTCAATACCGACGCCTTCGTTGCCGCCGACTCTCTGCAACAAACGCCCGACTACATGTCCGccc                                                >  2:365015/1‑150 (MQ=255)                                                                                                        tgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                             >  1:846420/1‑150 (MQ=255)                                                                                                         gctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCCACAACCAGTTGTTGGCGTACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGGCGCGCTGCAAAAAACGCTCGACTACATGGACGCCCAccc                                            >  1:782108/1‑150 (MQ=255)                                                                                                          ctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTAATCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCg                                           <  1:124473/150‑1 (MQ=255)                                                                                                            gAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGa                                         <  2:1488439/150‑1 (MQ=255)                                                                                                            gAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGa                                         <  2:614407/150‑1 (MQ=255)                                                                                                                aacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGTAGTgc                                     <  1:52633/150‑1 (MQ=255)                                                                                                                aacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgc                                     <  1:679454/150‑1 (MQ=255)                                                                                                                 acaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcg                                    <  2:846420/150‑1 (MQ=255)                                                                                                                  caaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 >  1:1594761/1‑120 (MQ=255)                                                                                                                  caaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                 <  2:1594761/120‑1 (MQ=255)                                                                                                                   aaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcggc                                  <  2:1114894/150‑1 (MQ=255)                                                                                                                   aaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcggc                                  <  1:1308430/150‑1 (MQ=255)                                                                                                                     gCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGt                                >  1:95345/1‑150 (MQ=255)                                                                                                                      cAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACCAGGACGAGGTCAATGGCGACCACATCGACGCCCACCCGGAGTGCGGCGTa                               >  1:1374681/1‑150 (MQ=255)                                                                                                                          aCGTGGGTTTCGGCCGGGCCAACAACCCGTTGTTGGTGCAGCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTgg                           <  2:150928/150‑1 (MQ=255)                                                                                                                            gTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTgggg                         <  1:1516327/150‑1 (MQ=255)                                                                                                                                gTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTgc                     >  2:1439256/1‑150 (MQ=255)                                                                                                                                        tgcgCCAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             >  2:1072450/4‑150 (MQ=255)                                                                                                                                        cGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTNGAg             <  2:1200827/150‑1 (MQ=255)                                                                                                                                        cGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             <  2:1605920/150‑1 (MQ=255)                                                                                                                                        cGGGCCAACAACCAGTCGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             <  2:972190/150‑1 (MQ=255)                                                                                                                                        cGGGCCAACAACCAGCTGTTGGCGCACCTCAATGGCAAGTATGCCTTGCTGCTCAATACCGATGCGAGTTCAGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg             <  1:1059198/150‑1 (MQ=255)                                                                                                                                         gggCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAgg            >  2:615391/1‑150 (MQ=255)                                                                                                                                         gggCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAgg            <  1:1457247/150‑1 (MQ=255)                                                                                                                                            ccAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTc         >  1:238114/1‑150 (MQ=255)                                                                                                                                            ccAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTc         <  2:1144809/150‑1 (MQ=255)                                                                                                                                              aacaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcgc       >  1:1499182/1‑150 (MQ=255)                                                                                                                                              aacaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcgc       >  2:1808002/1‑150 (MQ=255)                                                                                                                                               acaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcgcg      >  1:1340348/1‑150 (MQ=255)                                                                                                                                                caacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGa     <  1:1808002/150‑1 (MQ=255)                                                                                                                                                 aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACACCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt    >  2:51736/1‑150 (MQ=255)                                                                                                                                                  acCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATg   <  2:1878792/150‑1 (MQ=255)                                                                                                                                                   ccAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATgg  <  2:1243199/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGG  >  NC_002947/3550970‑3551265

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