Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,694,441	C→G	5.5%	R772R (CGC→CGG)	PP_4154 →	class 3 aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,694,441	0	C	G	5.5%	237.1 / 2.9	73	R772R (CGC→CGG)	PP_4154	class 3 aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (34/35);  new base G (2/2);  total (36/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 = 2.63e-01

GGCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATGGG  >  NC_002947/4694297‑4694587                                                                                                                                                  |                                                                                                                                                   ggCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCGCGCAGGTAGCCTGCGGGCAg                                                                                                                                                <  1:303347/149‑1 (MQ=255) ggCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCGCGCAGGTAGCCTGCGGGCAg                                                                                                                                                <  7:161287/149‑1 (MQ=255)    cGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGcc                                                                                                                                              >  2:173170/1‑148 (MQ=255)      aCAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTAt                                                                                                                                           <  4:188699/149‑1 (MQ=255)            ccgccgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAg                                                                                                                                     <  5:68602/149‑1 (MQ=255)             cgccgACTACCTGCAAGAAGTGGACGCCAATCTGGCTGACCTCGACGCCCGTGGCCGTGAGCGGGCGGGCATCATCTGCGAACCGGTGTACGGCAATGCCCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAg                                                                                                                                     <  1:77542/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                                 <  7:283330/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                                 <  8:257649/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                                 <  8:169793/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                                 <  7:15698/149‑1 (MQ=255)                    tACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCgcgc                                                                                                                             >  8:24891/1‑149 (MQ=255)                      ccTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGccc                                                                                                                           >  6:200131/1‑149 (MQ=255)                       cTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcg                                                                                                                          >  6:31111/1‑149 (MQ=255)                         gCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgcg                                                                                                                        >  1:286483/1‑149 (MQ=255)                         gCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgc                                                                                                                         >  2:109304/1‑148 (MQ=255)                         gCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgc                                                                                                                         <  1:109304/148‑1 (MQ=255)                         gCAAGACGTCGACGCCAAGCTCGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgcg                                                                                                                        >  3:179721/1‑149 (MQ=255)                            aGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCgcggcg                                                                                                                     >  4:204397/1‑149 (MQ=255)                                gTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                                 <  3:265727/149‑1 (MQ=255)                                gTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                                 <  4:149647/149‑1 (MQ=255)                                   gACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGtgtg                                                                                                              >  3:33092/1‑149 (MQ=255)                                   gACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGtgtg                                                                                                              >  5:312193/1‑149 (MQ=255)                                    aCGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGtgtg                                                                                                              >  1:81296/1‑148 (MQ=255)                                       ccAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATc                                                                                                          >  4:193491/1‑149 (MQ=255)                                        cAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATg                                                                                                                                          >  8:190105/1‑116 (MQ=255)                                         aaGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGc                                                                                                        <  3:193491/149‑1 (MQ=255)                                          aGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGcc                                                                                                       >  1:306073/1‑149 (MQ=255)                                          aGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGcc                                                                                                       <  3:182032/149‑1 (MQ=255)                                             tggctgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCcgac                                                                                                    >  4:164760/1‑149 (MQ=255)                                                 tgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                                 <  1:11098/132‑1 (MQ=255)                                                 tgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                                 >  2:11098/1‑132 (MQ=255)                                                 tgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAg                                                                                                >  6:73190/1‑149 (MQ=255)                                                   aCCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAg                                                                                                <  3:5790/147‑1 (MQ=255)                                                   aCCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGt                                                                                               <  5:73190/148‑1 (MQ=255)                                                   aCCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  2:290290/1‑149 (MQ=255)                                                    ccTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  4:112189/1‑148 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGaa                                                                                                 >  4:261292/1‑144 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  8:120147/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  5:146708/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  7:137533/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                              >  2:209355/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCCGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGaa                                                                                                 <  3:261292/144‑1 (MQ=255)                                                             ccGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgtgt                                                                                                               >  7:262307/1‑122 (MQ=255)                                                             ccGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgtgt                                                                                                               <  8:262307/122‑1 (MQ=255)                                                                   ccGTCAGCTGGCGGGCATCATATGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGgcc                                                                              >  2:305143/1‑149 (MQ=255)                                                                       cAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCt                                                                          >  7:187365/1‑149 (MQ=255)                                                                        aGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTg                                                                         >  3:172730/1‑149 (MQ=255)                                                                           tGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGc                                                                      <  2:306073/149‑1 (MQ=255)                                                                             gCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGc                                                                      <  4:245704/147‑1 (MQ=255)                                                                             gCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGc                                                                      <  6:312193/147‑1 (MQ=255)                                                                                   atcatcTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                              >  2:322183/1‑149 (MQ=255)                                                                                           cGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                          <  3:164760/145‑1 (MQ=255)                                                                                           cGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                          <  5:200131/145‑1 (MQ=255)                                                                                               ccGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGa                                                  <  7:120147/149‑1 (MQ=255)                                                                                               ccGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGa                                                  <  2:286483/149‑1 (MQ=255)                                                                                                  gtgtACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                              >  8:92491/1‑134 (MQ=255)                                                                                                  gtgtACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                              <  7:92491/134‑1 (MQ=255)                                                                                                       cGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAgg                                          <  5:31111/149‑1 (MQ=255)                                                                                                       cGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAg                                           <  4:226790/148‑1 (MQ=255)                                                                                                           aaTG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgt                                      >  4:272450/1‑149 (MQ=255)                                                                                                                ccGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgtggt                                   <  1:173170/148‑1 (MQ=255)                                                                                                                ccGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgtggt                                   <  1:305143/148‑1 (MQ=255)                                                                                                                       tATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGAca                           <  3:272450/149‑1 (MQ=255)                                                                                                                           tcGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGATATGGCCGCATGGGCGAGTAGTGGTGGGGGTACGAGGACCAGGGGGTGGTCCGGGAcagcat                       >  4:36897/1‑149 (MQ=255)                                                                                                                               tGCGCGCAGGTAGCCTGCGGGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTAc                    >  1:101559/1‑148 (MQ=255)                                                                                                                               tGCGCGCAGGTAGCCTGCGGGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTAc                    >  7:173148/1‑148 (MQ=255)                                                                                                                               tGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgtgg                                    >  2:24374/1‑132 (MQ=255)                                                                                                                               tGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgtgg                                    <  1:24374/132‑1 (MQ=255)                                                                                                                                       ggCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCaa           <  4:33092/149‑1 (MQ=255)                                                                                                                                        gCTACCTGCGCGGAGCCTATGCCAAGGTCCGGGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGGAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGTGGCTTCGAGCTACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAg          <  7:104278/149‑1 (MQ=255)                                                                                                                                         cTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAgg         <  1:290290/149‑1 (MQ=255)                                                                                                                                         cTACCTGCGCGCAGCCTATGCCAAGGGCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAgg         <  8:187365/149‑1 (MQ=255)                                                                                                                                            ccTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCa      >  2:241797/1‑149 (MQ=255)                                                                                                                                               gcgcgcAGCCTATGCCAAGGTCCGCTCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATgg   <  4:179721/149‑1 (MQ=255)                                                                                                                                                cgcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                          >  2:60266/1‑93 (MQ=255)                                                                                                                                                cgcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                          <  1:60266/93‑1 (MQ=255)                                                                                                                                                cgcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATgg   <  3:204397/148‑1 (MQ=255)                                                                                                                                                cgcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATgg   <  3:112189/148‑1 (MQ=255)                                                                                                                                                 gcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATggg  >  8:4068/1‑148 (MQ=255)                                                                                                                                                 gcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATgg   >  5:116779/1‑147 (MQ=255)                                                                                                                                                 gcgcAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATgg   <  6:116779/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGGGCATGGG  >  NC_002947/4694297‑4694587

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