Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,800,110	T→C	6.6%	I110V (ATC→GTC) ‡	PP_4227 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,800,110	0	T	C	6.6%	261.3 / 6.7	91	I110V (ATC→GTC) ‡	PP_4227	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/38);  new base C (2/4);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

ATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCAG  >  NC_002947/4799964‑4800257                                                                                                                                                   |                                                                                                                                                    aTACTACTTTCACAGTTCTCCACCCTCGAACCATATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGGTGCTGGCGAGGAtc                                                                                                                                                    <  7:281120/148‑1 (MQ=255)  tactacTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGAtctc                                                                                                                                                  >  8:2074/1‑149 (MQ=255)   actactTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCg                                                                                                                                                 >  2:257307/1‑149 (MQ=255)      aattGCACACATCTCCACCCTCCAACCCTATAGCTTTTCTGACTTCTAAAAACACTATTGTCGGGTAAATAGCGGCTGATGAGATGCATGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                              <  5:103484/147‑1 (MQ=37)       ctTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCAc                                                                                                                                             >  6:175393/1‑149 (MQ=255)       ctTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCag                                                                                                                                             <  7:237790/149‑2 (MQ=255)        ttGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACg                                                                                                                                            >  5:266051/1‑149 (MQ=255)           cacaGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCa                                                                                                                                          <  7:106853/148‑1 (MQ=255)            acaGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGc                                                                                                                                        <  3:3026/149‑1 (MQ=255)              aGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGc                                                                                                                                        >  1:55113/1‑147 (MQ=255)               gATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAg                                                                                                                                     >  7:202393/1‑149 (MQ=255)                   tcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctga                                                                                                                                 <  6:266051/149‑1 (MQ=255)                   tcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                              <  5:257440/136‑1 (MQ=255)                   tcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                              >  6:257440/1‑136 (MQ=255)                   tcCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAAAACAGAATTGGCCGCGAAAAAGCGCCGGATGAGATGCACTGCTTCGTCACGGCGTCCATGACGCCATCAAGCATGTGCTTGCTGGCGAGGATCTCGCCACCCAGCGAGcaga                                                                                                                                 >  3:234951/1‑149 (MQ=37)                     cACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTTGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgacc                                                                                                                               <  4:187754/149‑1 (MQ=255)                     cACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGTGAGGATCTCGCCACGCAGCGAGctgac                                                                                                                                <  7:124527/148‑1 (MQ=255)                     cACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGctgacc                                                                                                                               >  4:282817/1‑149 (MQ=255)                     cACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGct                                                                                                                                   <  8:164026/145‑1 (MQ=255)                            gAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGTGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAg                                                                                                                         <  1:297554/148‑1 (MQ=255)                              aCCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGct                                                                                                                                   <  4:218308/136‑1 (MQ=255)                              aCCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGct                                                                                                                                   >  3:218308/1‑136 (MQ=255)                                  tataGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgccatga                                                                                                                  >  1:307133/1‑148 (MQ=255)                                       cTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgc                                                                                                                       <  3:88099/139‑1 (MQ=255)                                       cTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGgc                                                                                                                       >  4:88099/1‑139 (MQ=255)                                       cTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATcc                                                                                                             <  6:309847/149‑1 (MQ=255)                                                 ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTg                                                                                                    >  4:152467/1‑148 (MQ=255)                                                 ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTg                                                                                                    >  4:28725/1‑148 (MQ=255)                                                 ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGa                                                                                                   >  4:55636/1‑149 (MQ=255)                                                 ttCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGAACGGCTTAGTCCCGGCGTCAATGACGACATCGTGGGTGTGGTTGCTGGCGAGGATCTCGCCCCGCAGCGAGCTGACCTGACAGGCCCTTCCATCCGGTACGCTg                                                                                                    >  4:60496/1‑148 (MQ=25)                                                  tCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGa                                                                                                   >  6:143133/1‑148 (MQ=255)                                                   cTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAAc                                                                                                 >  3:57970/1‑149 (MQ=255)                                                    tAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGt                                                                                                          >  2:249258/1‑139 (MQ=255)                                                    tAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGt                                                                                                          <  1:249258/139‑1 (MQ=255)                                                      aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTgg                                                                                              >  3:353716/1‑149 (MQ=255)                                                      aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTgg                                                                                              >  3:216821/1‑149 (MQ=255)                                                      aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTgg                                                                                              >  1:42361/1‑149 (MQ=255)                                                      aaCACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTg                                                                                               <  5:175393/148‑1 (MQ=255)                                                       acacAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGt                                                                                             <  1:257307/149‑1 (MQ=255)                                                       acacAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGt                                                                                             <  3:152467/149‑1 (MQ=255)                                                        cacaGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTc                                                                                            >  7:222859/1‑149 (MQ=255)                                                         acaGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcg                                                                                           >  3:270354/1‑149 (MQ=255)                                                           aGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggc                                                                                         >  4:356991/1‑149 (MQ=255)                                                            gTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGGTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                        <  5:375257/149‑1 (MQ=255)                                                             tATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGGTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcgg                                                                                       <  1:274836/149‑1 (MQ=255)                                                               ttGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGaa                                                                                     <  3:210187/149‑1 (MQ=255)                                                                  gCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGCCCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGaaaaa                                                                                  >  5:278794/1‑149 (MQ=255)                                                                  gCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCCTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGAAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGCCGGaaaaa                                                                                  >  3:47797/1‑149 (MQ=255)                                                                    cGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATg                                                                                >  4:140747/1‑149 (MQ=255)                                                                     gggAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCt                                                                                                     <  4:364485/127‑1 (MQ=255)                                                                     gggAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCt                                                                                                     >  3:364485/1‑127 (MQ=255)                                                                      tgAAAAATCGGCGGATGAGATTGACGGCTTCGTCGCGGCGTCCATTACGCCATCGTGCTTGTTGTTGCTGGCGAGGATCTCGACACGCAGCTAGCTGAAGTGACAGGCCATGCCATTCGGTACGCTGAACTGGTCGGCGGAAAAATgc                                                                               <  4:234951/147‑1 (MQ=18)                                                                       gAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATgcgc                                                                             >  1:113592/1‑149 (MQ=255)                                                                           aaGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTcc                                                                         >  2:173141/1‑149 (MQ=255)                                                                            aGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTccc                                                                        <  2:42361/149‑1 (MQ=255)                                                                                  ggATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCt                                                                  >  4:164517/1‑149 (MQ=255)                                                                                   gATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCtt                                                                 >  3:4537/1‑149 (MQ=255)                                                                                     tGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTgg                                                               <  4:57970/149‑1 (MQ=255)                                                                                     tGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTgg                                                               <  3:28725/149‑1 (MQ=255)                                                                                       agagGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCg                                                             <  3:140747/149‑1 (MQ=255)                                                                                       agaTGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAACTGCGCCTCCCCGCCCTTGGc                                                              <  3:60496/148‑1 (MQ=255)                                                                                       agaTGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGc                                                              <  5:143133/148‑1 (MQ=255)                                                                                       agaTGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGc                                                              <  3:55636/148‑1 (MQ=255)                                                                                           gCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                        <  5:288207/118‑1 (MQ=255)                                                                                           gCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTcggcg                                                                                        >  6:288207/1‑118 (MQ=255)                                                                                             aCGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGc                                                       >  8:296962/1‑149 (MQ=255)                                                                                               ggCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCgg                                                     <  2:55113/149‑1 (MQ=255)                                                                                               ggCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCgg                                                     <  8:202393/149‑1 (MQ=255)                                                                                                 cTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCt                                                   >  8:72843/1‑149 (MQ=255)                                                                                                  ttCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                              >  5:277354/1‑57 (MQ=18)                                                                                                  ttCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGATCCTCGCCa                                                                                                                                              <  6:277354/57‑1 (MQ=18)                                                                                                      tCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCTAGTATCTCGCCACGCAGCGAGCTTACCTGTCAGGCCATGCCATACGGTACGCTGAACTGGTCGTCGGAAAAATGCGCCTCCCCGCACTTGGCGGTCAGCGGCtgctgc                                              <  2:164324/149‑1 (MQ=14)                                                                                                        ccGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGTGGCTGCGGCaa                                            >  7:81845/1‑149 (MQ=255)                                                                                                        ccGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCaa                                            <  1:173141/149‑1 (MQ=40)                                                                                                         cGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAAc                                           >  5:33580/1‑149 (MQ=40)                                                                                                             gTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGt                                       >  7:45955/1‑149 (MQ=255)                                                                                                              tCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTc                                      <  4:353716/149‑1 (MQ=255)                                                                                                               ccATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCg                                     <  8:45955/149‑1 (MQ=255)                                                                                                                cATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGa                                    >  5:276526/1‑149 (MQ=255)                                                                                                                     cGCCATCGTTCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGc                               <  4:47797/149‑1 (MQ=255)                                                                                                                          tCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGc                          >  8:128461/1‑149 (MQ=39)                                                                                                                                ttGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGGCATTGCATACGGTACGCTGAACTGGTGGGCGGAAAAATGCGCCTACCCGCCCTTGGCGGTAAGCGGCTGCGGCAACAGGTCGATAAGCCCGGCGACCa                     >  1:90410/1‑148 (MQ=14)                                                                                                                                    ggTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                >  4:282660/1‑149 (MQ=39)                                                                                                                                    ggTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                >  3:12235/1‑149 (MQ=39)                                                                                                                                    ggTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                <  2:113592/149‑1 (MQ=39)                                                                                                                                     gTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgcg               >  2:274632/1‑149 (MQ=39)                                                                                                                                     gTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATAcgc                <  3:164517/148‑1 (MQ=39)                                                                                                                                      ttGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACgcgc              <  3:282817/149‑1 (MQ=39)                                                                                                                                          tGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGAc          >  1:36446/1‑149 (MQ=39)                                                                                                                                            gCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGAcgc        >  4:129337/1‑149 (MQ=39)                                                                                                                                            gCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGAcgc        <  6:33580/149‑1 (MQ=39)                                                                                                                                              gAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGAGCCGAcgccg      >  4:41506/1‑149 (MQ=255)                                                                                                                                                  aTCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCAg  >  8:126268/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATACTACTTGCACAGATCTCCACCCTCGAACCCTATAGCTTTTCTGACTTCTAAACACAGTATTGGCCGGGAAAAAGCGGCGGATGAGATGCACGGCTTCGTCCCGGCGTCCATGACGCCATCGTGCTTGTGGTTGCTGGCGAGGATCTCGCCACGCAGCGAGCTGACCTGACAGGCCATGCCATCCGGTACGCTGAACTGGTCGGCGGAAAAATGCGCCTCCCCGCCCTTGGCGGTCAGCGGCTGCGGCAACTGGTCGATCAGCCCGGCGACCATACGCGCCGACGCCGCCAG  >  NC_002947/4799964‑4800257

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