Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,694,434	A→G	5.5%	Y770C (TAC→TGC) ‡	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,434	0	A	G	5.5%	214.4 / 2.1	72	Y770C (TAC→TGC) ‡	PP_4154	class 3 aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base G (2/2);  total (35/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 = 9.23e-01

GGCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGG  >  NC_002947/4694297‑4694580                                                                                                                                           |                                                                                                                                                   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                                                                                                                          <  8:257649/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                          <  8:169793/149‑1 (MQ=255)                cgACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTcc                                                                                                                          <  7:283330/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                                                                                                                  <  1:109304/148‑1 (MQ=255)                         gCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgc                                                                                                                  >  2:109304/1‑148 (MQ=255)                         gCAAGACGTCGACGCCAAGCTCGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCcgcg                                                                                                                 >  3:179721/1‑149 (MQ=255)                            aGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCgcggcg                                                                                                              >  4:204397/1‑149 (MQ=255)                                gTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                          <  4:149647/149‑1 (MQ=255)                                gTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                          <  3:265727/149‑1 (MQ=255)                                   gACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGtgtg                                                                                                       >  5:312193/1‑149 (MQ=255)                                   gACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGtgtg                                                                                                       >  3:33092/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                                                                                                <  3:182032/149‑1 (MQ=255)                                          aGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGcc                                                                                                >  1:306073/1‑149 (MQ=255)                                             tggctgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCcgac                                                                                             >  4:164760/1‑149 (MQ=255)                                                 tgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                          >  2:11098/1‑132 (MQ=255)                                                 tgACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgt                                                                                                          <  1:11098/132‑1 (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                                                                                       >  7:137533/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                       >  2:209355/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                       >  8:120147/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTg                                                                                       >  5:146708/1‑147 (MQ=255)                                                     cTCGACGCCCGTGGCCGTCAGCTGGCCGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGaa                                                                                          <  3:261292/144‑1 (MQ=255)                                                             ccGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgtgt                                                                                                        <  8:262307/122‑1 (MQ=255)                                                             ccGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGgtgt                                                                                                        >  7:262307/1‑122 (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                                                               <  6:312193/147‑1 (MQ=255)                                                                             gCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGc                                                               <  4:245704/147‑1 (MQ=255)                                                                                   atcatcTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                       >  2:322183/1‑149 (MQ=255)                                                                                           cGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                   <  5:200131/145‑1 (MQ=255)                                                                                           cGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTgg                                                   <  3:164760/145‑1 (MQ=255)                                                                                               ccGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGa                                           <  2:286483/149‑1 (MQ=255)                                                                                               ccGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGa                                           <  7:120147/149‑1 (MQ=255)                                                                                                  gtgtACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                       <  7:92491/134‑1 (MQ=255)                                                                                                  gtgtACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACtt                                                       >  8:92491/1‑134 (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:305143/148‑1 (MQ=255)                                                                                                                ccGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCgtggt                            <  1:173170/148‑1 (MQ=255)                                                                                                                       tATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGAca                    <  3:272450/149‑1 (MQ=255)                                                                                                                           tcGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGATATGGCCGCATGGGCGAGTAGTGGTGGGGGTACGAGGACCAGGGGGTGGTCCGGGAcagcat                >  4:36897/1‑149 (MQ=255)                                                                                                                               tGCGCGCAGGTAGCCTGCGGGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTAc             >  7:173148/1‑148 (MQ=255)                                                                                                                               tGCGCGCAGGTAGCCTGCGGGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTAc             >  1:101559/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)                                                                                                                                           |                                                                                                                                                   GGCCGACAGCGCCGCCGACTACCTGCAAGACGTCGACGCCAAGCTGGCTGACCTCGACGCCCGTGGCCGTCAGCTGGCGGGCATCATCTGCGAACCGGTGTACGGCAATG‑CCGGCGGTATCTCGCTGCCCGCAGGCTACCTGCGCGCAGCCTATGCCAAGGTCCGCGCCCGCGGCGGGGTGTGCATCGCCGACGAAGTGCAGGTAGGCTATGGCCGCCTGGGCGAGTACTTCTGGGGCTTCGAGGAACAGGGCGTGGTGCCGGACATCATTACCATGGCCAAGG  >  NC_002947/4694297‑4694580

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