Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,837,686	A→C	5.8%	K102T (AAG→ACG)	PP_1642 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,837,686	0	A	C	5.8%	261.2 / 4.7	81	K102T (AAG→ACG)	PP_1642	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base C (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGGCTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCGGG  >  NC_002947/1837538‑1837834                                                                                                                                                     |                                                                                                                                                     gcgGGGCTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCaa                                                                                                                                                      <  3:32337/149‑1 (MQ=255)   ggggCTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGc                                                                                                                                                    <  7:75409/149‑1 (MQ=255)   ggggCTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGc                                                                                                                                                    >  2:211594/1‑149 (MQ=255)       cTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCTAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCACATTATCCGCCACCTGGGCAAGCTCaa                                                                                                                                                <  7:199106/149‑1 (MQ=255)       cTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCaa                                                                                                                                                <  5:315717/149‑1 (MQ=255)          aaGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAaga                                                                                                                                              <  6:144656/148‑1 (MQ=255)          aaGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAag                                                                                                                                               <  6:201631/147‑1 (MQ=255)          aaGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTAGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAaga                                                                                                                                              <  1:180056/148‑1 (MQ=255)            gCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCg                                                                                                                                           >  2:103046/1‑149 (MQ=255)             cacaGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGt                                                                                                                                          <  5:869/149‑1 (MQ=255)              acaGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTg                                                                                                                                         <  1:246006/149‑1 (MQ=255)                  ccTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCACGCTCTTGAGCTTGccca                                                                                                                                     <  3:252925/149‑3 (MQ=255)                      gtggATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCaa                                                                                                                                 >  5:95057/1‑149 (MQ=255)                         gATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATgc                                                                                                                              >  2:262243/1‑149 (MQ=255)                             ccAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCaa                                                                                                                          <  1:51705/149‑1 (MQ=255)                             ccAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCaa                                                                                                                          <  7:17886/149‑1 (MQ=255)                              cAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCaaa                                                                                                                         >  3:335382/1‑149 (MQ=255)                              cAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCaaa                                                                                                                         >  6:200597/1‑149 (MQ=255)                                aGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATg                                                                                                                       <  2:35779/149‑1 (MQ=255)                                 gTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATgg                                                                                                                      >  3:205623/1‑149 (MQ=255)                                     ccGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCTGAGAAAGTGGTGCTGATGTGCGCCGAGCATCTGGAGCTGTTGATGCAGTATTAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCt                                                                                                                  <  5:95783/149‑1 (MQ=255)                                      cGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCACGCTCAAGCGCGTGCCTCGCAATGCGCAAATGGTGCTg                                                                                                                 >  8:348335/1‑149 (MQ=255)                                                   ggTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATg                                                                                                                       <  4:204183/130‑1 (MQ=255)                                                   ggTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATg                                                                                                                       >  3:204183/1‑130 (MQ=255)                                                   ggTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCAGGCCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTCCTGGATATCGCCAAgg                                                                                                    >  3:26668/1‑149 (MQ=255)                                                    gtgtTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGgc                                                                                                   >  3:186768/1‑149 (MQ=255)                                                      gtTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCGAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGgcg                                                                                                  >  1:256897/1‑148 (MQ=255)                                                      gtTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATgc                                                                                                                              >  8:277673/1‑120 (MQ=255)                                                      gtTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATgc                                                                                                                              <  7:277673/120‑1 (MQ=255)                                                        tCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTgg                                                                                                                >  6:68156/1‑132 (MQ=255)                                                        tCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTgg                                                                                                                <  5:68156/132‑1 (MQ=255)                                                         cTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTgg                                                                                                                <  2:327427/131‑1 (MQ=255)                                                         cTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTgg                                                                                                                >  1:327427/1‑131 (MQ=255)                                                               cTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGtt                                                                                        >  2:369721/1‑149 (MQ=255)                                                                ttCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGttt                                                                                       <  1:251514/149‑1 (MQ=255)                                                                ttCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGttt                                                                                       >  7:72543/1‑149 (MQ=255)                                                                 tCGACCCGGAGATAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGtttt                                                                                      >  5:71139/1‑149 (MQ=255)                                                                    aCCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGt                                                                                   <  1:103046/149‑1 (MQ=255)                                                                       cGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCGAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCa                                                                                <  2:256897/149‑1 (MQ=255)                                                                       cGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCa                                                                                >  8:253828/1‑149 (MQ=255)                                                                              aGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCg                                                                         >  8:220716/1‑149 (MQ=255)                                                                                   tgcggctGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCCAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGAt                                                                    <  6:93843/143‑1 (MQ=255)                                                                                    gCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGAtt                                                                   >  1:229581/1‑149 (MQ=255)                                                                                      tGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGc                                                                                                                                   <  8:52109/83‑1 (MQ=255)                                                                                      tGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGc                                                                                                                                   >  7:52109/1‑83 (MQ=255)                                                                                        aTGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGGAGTTTTGGTGCATTCATCGCCGATTGGc                                                                <  1:211594/148‑1 (MQ=255)                                                                                        aTGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTgg                                                                 <  6:71139/147‑1 (MQ=255)                                                                                         tGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCa                                                              <  4:186768/149‑1 (MQ=255)                                                                                             cgcCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGa                                                          >  3:235000/1‑149 (MQ=255)                                                                                                     aTCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACAtt                                                  >  4:23745/1‑149 (MQ=255)                                                                                                          gAGCGTTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTc                                               <  7:220716/147‑1 (MQ=255)                                                                                                          gAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCgg                                             <  5:200597/149‑1 (MQ=255)                                                                                                               gTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTtgt                                        <  2:229581/149‑1 (MQ=255)                                                                                                                ttGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTtgtg                                       >  2:49142/1‑149 (MQ=255)                                                                                                                ttGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTtgtg                                       >  4:271327/1‑149 (MQ=255)                                                                                                                  gATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTgg                                      >  1:349911/1‑148 (MQ=255)                                                                                                                   aTGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGaa                                    <  3:271327/149‑1 (MQ=255)                                                                                                                        ggATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCTTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGt                                                <  6:225653/132‑1 (MQ=255)                                                                                                                        ggATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGt                                                >  5:225653/1‑132 (MQ=255)                                                                                                                        ggATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTAcc                               <  8:72543/149‑1 (MQ=255)                                                                                                                        ggATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTAcc                               >  7:249575/1‑149 (MQ=255)                                                                                                                            gAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGt                                                                                   >  3:15634/1‑93 (MQ=255)                                                                                                                            gAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGt                                                                                   <  4:15634/93‑1 (MQ=255)                                                                                                                               cgcATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCaa                        <  4:335382/149‑1 (MQ=255)                                                                                                                                gcATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGACGTACCTGGCCAAg                       >  2:70673/1‑149 (MQ=255)                                                                                                                                  aTTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCg                                                       >  6:15014/1‑115 (MQ=255)                                                                                                                                  aTTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCg                                                       <  5:15014/115‑1 (MQ=255)                                                                                                                                     aTCCGCCACCTGGGCAAGCTCAAGAGCGTGCGACGCAATGCGCAATTGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACg                   <  4:26668/148‑1 (MQ=255)                                                                                                                                     aTCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACg                   <  1:262243/148‑1 (MQ=255)                                                                                                                                      tCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACgg                  <  7:253828/148‑1 (MQ=255)                                                                                                                                        cGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCaa               >  7:172513/1‑149 (MQ=255)                                                                                                                                        cGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCaa               <  6:85121/149‑1 (MQ=255)                                                                                                                                        cGCCACCTGGGCAAGCGCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGGCGGCTTGTGGAAGTACCTGGCCAAGCACGGCaa               <  7:348335/149‑1 (MQ=255)                                                                                                                                          ccACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGTAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAAc              >  8:23387/1‑148 (MQ=255)                                                                                                                                             ccTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTTATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGt          <  1:369721/149‑1 (MQ=255)                                                                                                                                              cTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGtt         <  6:95057/149‑1 (MQ=255)                                                                                                                                               gtggCACGCTCTTGAGCTTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGaa                                    <  2:326183/119‑1 (MQ=255)                                                                                                                                               gtggCACGCTCTTGAGCTTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGaa                                    >  1:326183/3‑121 (MQ=255)                                                                                                                                                tggCACGCTCTTGAGCTTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTgg       >  6:137858/2‑149 (MQ=255)                                                                                                                                                   cAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGTAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCg    <  7:23387/149‑1 (MQ=255)                                                                                                                                                   cAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCg    >  5:374794/1‑149 (MQ=255)                                                                                                                                                   cAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCg    >  4:207955/1‑149 (MQ=255)                                                                                                                                                    aaGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCgg   <  6:374794/149‑1 (MQ=255)                                                                                                                                                     aGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCggg  <  4:235000/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGGGGCTCAAGCACAGCCTGGTGGATGCCAAGTGGCCGGCGTTCGAGCAGGTGTTCTTTGGCTTCGACCCGGAGAAAGTGGTGCTGATGGGCGCCGAGCATCTGGAGCGGTTGATGCAGGATGAGCGCATTATCCGCCACCTGGGCAAGCTCAAGAGCGTGCCACGCAATGCGCAAATGGTGCTGGATATCGCCAAGGCGCATGGCAGTTTTGGTGCATTCATCGCCGATTGGCCAGTGACCGACATTGTCGGCTTGTGGAAGTACCTGGCCAAGCACGGCAACCAGTTGGGCGGG  >  NC_002947/1837538‑1837834

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