Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,357,137	C→A	6.5%	G81W (GGG→TGG)	PP_2071 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,357,137	0	C	A	6.5%	186.6 / 3.4	62	G81W (GGG→TGG)	PP_2071	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base A (2/2);  total (36/26)
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.86e-01

ATCCAGACATCGCTCTCTGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCATCTCGGC  >  NC_002947/2356995‑2357284                                                                                                                                                |                                                                                                                                                    aTCCAGACATCGCTCTCTGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGATAGACGCCGCGGGGCACAACATAGAGTT‑CCCCGGCGtt                                                                                                                                               >  6:24772/1‑149 (MQ=255)     aGACATCGCTCTCTGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTcaca                                                                                                                                           >  3:290837/1‑149 (MQ=255)         aTCGCTCTCTGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCa                                                                                                                                       >  2:109579/1‑149 (MQ=255)             ctctctGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCa                                                                                                                                   <  4:226280/149‑1 (MQ=255)                 ctgctgTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCa                                                                                                                               >  1:105886/1‑149 (MQ=255)                         cGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAg                                                                                                                       >  3:351915/1‑149 (MQ=255)                             tCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTc                                                                                                                    >  6:134656/1‑148 (MQ=255)                               cccccTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGa                                                                                                                  >  1:197665/1‑148 (MQ=255)                               cacccTTCATGACCGGTATTCAGAATGCCCAGCTTTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCTGTCCGGTGATCGACGACGCGTGGCACCACATAGAGTT‑CACCGCGTTTCACATGCACGCAGCCATCACTCAGGTcta                                                                                                                  >  8:163265/3‑146 (MQ=255)                                 cccTTCATGACCGGTATTCAGAATGCCCCGCGTTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCg                                                                                                                   <  4:164543/145‑1 (MQ=255)                                 cccTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCg                                                                                                                   >  3:164543/1‑145 (MQ=255)                                       aTGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAgg                                                                                                          <  7:206411/148‑1 (MQ=255)                                        tGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGt                                                                                                        >  8:314222/1‑149 (MQ=255)                                         gACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGc                                                                                                                                        >  4:228795/1‑116 (MQ=255)                                            cGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCt                                                                                                    <  8:246784/149‑1 (MQ=255)                                            cGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCt                                                                                                    <  1:261956/149‑1 (MQ=255)                                            cGGTATTCAGAATGCCCCGCGGTTCGAGCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCCGCCATCCGGCAGGTCGATACGCAGGGTGCCt                                                                                                    >  6:347067/1‑149 (MQ=255)                                              gTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTc                                                                                                  <  5:341895/149‑1 (MQ=255)                                                    aGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCAc                                                                                            <  8:308976/149‑1 (MQ=255)                                                      aaTGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCAc                                                                                            >  5:201930/1‑147 (MQ=255)                                                      aaTGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCa                                                                                          <  2:105886/149‑1 (MQ=255)                                                         gCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAgg                                                                                        >  5:45498/1‑148 (MQ=255)                                                           cccGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGCCGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAAt                                                                                     >  6:81616/1‑149 (MQ=255)                                                           cccGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGt                                                                                                        <  4:326390/130‑1 (MQ=255)                                                           cccGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGt                                                                                                        >  3:326390/1‑130 (MQ=255)                                                           cccGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAAt                                                                                     >  7:219040/1‑149 (MQ=255)                                                              gcgGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCt                                                                                  <  4:351915/149‑1 (MQ=255)                                                                  ttCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCAGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTtcgt                                                                              >  5:429263/1‑149 (MQ=255)                                                                  ttCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCAGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTtcgt                                                                              >  5:435026/1‑149 (MQ=255)                                                                    cGATCATCATCATTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTtcgtc                                                                             <  8:349757/148‑1 (MQ=255)                                                                               aGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAAGAATGCTTCGTCGGTCTCCGGGTg                                                                 <  7:270376/149‑1 (MQ=255)                                                                                 tttCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGTAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTgcg                                                               >  7:48866/1‑149 (MQ=255)                                                                                 tttCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTgcg                                                               >  4:204769/1‑149 (MQ=255)                                                                                      cttcgcCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGc                                                          >  2:306692/1‑149 (MQ=255)                                                                                          gcCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGc                                                          >  2:304521/1‑145 (MQ=255)                                                                                            cTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGAt                                                     >  3:269304/1‑148 (MQ=255)                                                                                             ttGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGtctc                                                                       >  4:406062/1‑129 (MQ=255)                                                                                             ttGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGtctc                                                                       <  3:406062/129‑1 (MQ=255)                                                                                               gcgcAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGa                                                   <  5:134656/147‑1 (MQ=255)                                                                                                  cAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCAACCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACAAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGAGGAATTCg                                              >  6:299892/1‑149 (MQ=255)                                                                                                     cGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATc                                                                                                                             <  1:168682/67‑1 (MQ=255)                                                                                                     cGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATc                                                                                                                             >  2:168682/1‑67 (MQ=255)                                                                                                      ggTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCtt                                          <  3:204769/149‑1 (MQ=255)                                                                                                      ggTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCtt                                          >  3:183634/1‑149 (MQ=255)                                                                                                      ggTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTG‑CCCCGGGGTTCACATGCACGCAGCCATCCTGCAGGTCGATACGCAGGGTGCCTTCCATCACCAGGAATGCTTCGTCTGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCtt                                          <  4:283133/149‑1 (MQ=255)                                                                                                        tCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCTGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTTCGCGTGCCATATGAATTCGCCTTCg                                        <  8:219040/149‑1 (MQ=255)                                                                                                          cGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGAt                                      >  2:4250/1‑149 (MQ=255)                                                                                                          cGGTGCTCGACGCCGCGGGGCACCACATAAAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCCTCCAGAACCAGAAATGCTTCGTAGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGAt                                      >  7:228935/1‑149 (MQ=255)                                                                                                               cTCGACGCCGCGGGGCACCACATAGAGTTCCCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGaa                                                                                      >  1:267644/1‑97 (MQ=255)                                                                                                               cTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGaa                                                                                      <  2:267644/96‑1 (MQ=255)                                                                                                               cTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCa                                 <  6:45498/149‑1 (MQ=255)                                                                                                                  gACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACg                               <  5:24772/148‑1 (MQ=255)                                                                                                                  gACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGa                              >  5:191172/1‑149 (MQ=255)                                                                                                                     gccgcGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCt                           >  1:389295/1‑149 (MQ=255)                                                                                                                       cgcgGGGCACCACATAGAGTT‑CCCAGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCTGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTg                         <  6:435026/149‑1 (MQ=255)                                                                                                                       cgcgGGGCACCACATAGAGTT‑CCCAGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTg                         <  6:429263/149‑1 (MQ=255)                                                                                                                          ggggCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAAc                      >  2:70001/1‑149 (MQ=255)                                                                                                                            gtcaccaaATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCATCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTg                    <  7:163265/147‑1 (MQ=255)                                                                                                                               accacATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTa                 <  1:109579/149‑1 (MQ=255)                                                                                                                                   caTAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCg             >  4:389230/1‑149 (MQ=255)                                                                                                                                     tAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGTGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGtt           <  3:20398/149‑1 (MQ=255)                                                                                                                                      agagTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGtt           >  4:426216/1‑148 (MQ=255)                                                                                                                                          gtgccccGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCAtc       <  2:153949/146‑1 (MQ=255)                                                                                                                                             ccccGGGGTTCACCTGCACGCAGGCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGACGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTc          <  8:281766/143‑1 (MQ=255)                                                                                                                                             ccccGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCATCTCg    >  2:133984/1‑149 (MQ=255)                                                                                                                                             caccGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGACGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTc          >  7:281766/3‑143 (MQ=255)                                                                                                                                              cccGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCATCTCg    <  5:85546/148‑1 (MQ=255)                                                                                                                                               ccGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCATCTCGGc  <  4:269304/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    ATCCAGACATCGCTCTCTGCTGTGCGATTCCCCCCTTCATGACCGGTATTCAGAATGCCCCGCGGTTCGATCATCATCAGTTTCGCTTCGCCTTGCGCAGCGGTCCGGTGCTCGACGCCGCGGGGCACCACATAGAGTT‑CCCCGGGGTTCACATGCACGCAGCCATCCGGCAGGTCGATACGCAGGGTGCCTTCCAGCACCAGGAATGCTTCGTCGGTCTCCGGGTGCGCGTGCCAGATGAATTCGCCTTCGATGCGCACGACCTTGAACTGGTAATCGTTCATCTCGGC  >  NC_002947/2356995‑2357284

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