Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,289,025	Δ1 bp	7.1%	coding (98/822 nt)	PP_3760 ←	methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,289,025	0	G	.	7.1%	315.6 / 16.5	98	coding (98/822 nt)	PP_3760	methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/53);  new base . (2/5);  total (40/58)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGG  >  NC_002947/4288880‑4289170                                                                                                                                                      |                                                                                                                                                  gTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAg                                                                                                                                                >  3:60479/1‑149 (MQ=255)  tCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAgc                                                                                                                                               >  7:114147/1‑149 (MQ=255)     cGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                            <  8:372981/149‑1 (MQ=255)     cGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                            <  6:319520/149‑1 (MQ=255)       gAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCg                                                                                                                                           >  4:378893/1‑148 (MQ=255)        aagaTCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                         <  2:423573/146‑1 (MQ=255)        aaCATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCgg                                                                                                                                          >  5:420749/1‑148 (MQ=255)        aaCATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGa                                                                                                                                         <  8:121345/149‑1 (MQ=255)          cATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAGCCTCACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                         <  1:348666/149‑1 (MQ=255)           aTCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaat                                                                                                                                       >  8:89715/1‑148 (MQ=255)           aTCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaat                                                                                                                                       >  5:101210/1‑148 (MQ=255)                gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                     <  5:221333/145‑1 (MQ=255)                gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                     <  5:437282/145‑1 (MQ=255)                gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGc                                                                                                                                 <  3:166991/149‑1 (MQ=255)                 cTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                                <  3:251729/149‑1 (MQ=255)                 cTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCCGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                                <  4:123584/149‑1 (MQ=255)                    aCCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGa                                                                                                                              >  4:407745/1‑148 (MQ=255)                      cGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAAt                                                                                                                           >  7:280186/1‑149 (MQ=255)                      cGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAAt                                                                                                                           >  6:272729/1‑149 (MQ=255)                          aTCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                                       <  7:404628/149‑1 (MQ=255)                          aTCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                                       <  5:272729/149‑1 (MQ=255)                          aTCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                                       <  6:420749/149‑1 (MQ=255)                           tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                      <  1:498/149‑1 (MQ=255)                           tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                      >  5:63261/1‑149 (MQ=255)                           tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                      <  3:232877/149‑1 (MQ=255)                             gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGTTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGATTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGc                                                                                                                     <  6:85369/148‑1 (MQ=255)                             gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCt                                                                                                                    >  2:315739/1‑149 (MQ=255)                              tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGGCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCt                                                                                                                    >  1:119546/1‑149 (MQ=255)                                 ggTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTAc                                                                                                                <  3:378893/149‑1 (MQ=255)                                 ggTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCCTCAGACAATCGAACTCACGCAGCGCGTGGAGAATCCGGCGCTA‑GATCGAAGCGCCGGAATAATAGCGGAAATCGTAGCAGTAc                                                                                                                >  1:123928/1‑149 (MQ=255)                                  gTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACt                                                                                                               >  1:122674/1‑149 (MQ=255)                                     cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCgg                                                                                                                               <  6:393921/130‑1 (MQ=255)                                     cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCgg                                                                                                                               >  5:393921/1‑130 (MQ=255)                                     cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                            <  4:251112/149‑1 (MQ=255)                                      tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                           >  6:128518/1‑149 (MQ=255)                                      tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                           >  5:110115/1‑149 (MQ=255)                                        cagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAgg                                                                                                          <  4:60479/148‑1 (MQ=255)                                            aaCTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAt                                                                                                     <  6:207069/149‑1 (MQ=255)                                             aCTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATc                                                                                                    <  6:110115/149‑1 (MQ=255)                                                 cATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGcc                                                                                                 <  4:223650/148‑1 (MQ=255)                                                 cATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                                <  8:126085/149‑1 (MQ=255)                                                 cATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCCGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                                >  6:52027/1‑149 (MQ=255)                                                    gAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTg                                                                                                             <  3:148099/133‑1 (MQ=255)                                                    gAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTg                                                                                                             >  4:148099/1‑133 (MQ=255)                                                    gAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                             >  1:270201/1‑149 (MQ=255)                                                    gAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                             >  6:122220/1‑149 (MQ=255)                                                      aCATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAAACGTAGCTGTACTTGAGGTAGATCGCCTCGAt                                                                                            >  1:345006/1‑148 (MQ=255)                                                        aTGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATcag                                                                                         >  8:395492/1‑149 (MQ=255)                                                        aTGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATcag                                                                                         <  1:392993/149‑1 (MQ=255)                                                         tGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                      <  5:367561/119‑1 (MQ=255)                                                         tGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                      >  6:367561/1‑119 (MQ=255)                                                            ccGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGt                                                                                     >  5:86622/1‑149 (MQ=255)                                                            ccGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGt                                                                                     >  3:387176/1‑149 (MQ=255)                                                                 tCGTGCAGGACGCGCTCCT‑GGAACGCAG‑‑ACAGCGTCAGTCAATCGAACTGACGCTGCGCGTGGAGAATCCGTCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTTGTTGGCGTCGATCAGCAGTCGGAt                                                                                <  2:298012/149‑1 (MQ=255)                                                                    tGCAGGACGCTCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTGGGAtctc                                                                             <  8:114147/149‑1 (MQ=255)                                                                       aGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTCGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGAt                                                                          <  1:90279/149‑1 (MQ=255)                                                                         gACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGt                                                                        <  8:282649/149‑1 (MQ=255)                                                                          aCGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTc                                                                       <  8:280186/149‑1 (MQ=255)                                                                          aCGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTCTCGCTGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGTGGTAGATGGCCTCGATCAGCAGTCGGATCTCGATGTc                                                                       <  2:123928/149‑1 (MQ=255)                                                                               cTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGgtgt                                                                  <  6:429075/149‑1 (MQ=255)                                                                                tCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                                 >  2:328642/1‑149 (MQ=255)                                                                                  cT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGc                                                               >  3:420454/1‑149 (MQ=255)                                                                                  cT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGc                                                               <  6:101210/149‑1 (MQ=255)                                                                                   t‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCg                                                              >  1:400347/1‑149 (MQ=255)                                                                                        aCGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCg                                                          >  2:397304/1‑149 (MQ=255)                                                                                         cGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGc                                                         <  3:407745/149‑1 (MQ=255)                                                                                           cAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                                       <  1:328642/149‑1 (MQ=255)                                                                                                    cGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                                <  5:122220/149‑1 (MQ=255)                                                                                                    cGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                                <  2:270201/149‑1 (MQ=255)                                                                                                     gTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAtt                                               <  8:23446/149‑1 (MQ=255)                                                                                                      tCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTc                                              >  8:386919/1‑149 (MQ=255)                                                                                                        aGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCg                                             <  7:395492/148‑1 (MQ=255)                                                                                                          aCAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                          <  2:119546/149‑1 (MQ=255)                                                                                                           cAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGAAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                          <  1:315739/148‑1 (MQ=255)                                                                                                             aTCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAg                                       >  5:373984/1‑149 (MQ=255)                                                                                                                      aCGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAg                              >  1:6966/1‑149 (MQ=255)                                                                                                                      aCGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAg                              <  5:52027/149‑1 (MQ=255)                                                                                                                       cGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGACATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAg                              <  7:296638/148‑1 (MQ=255)                                                                                                                        gCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGcc                            <  6:63261/149‑1 (MQ=255)                                                                                                                              cgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGacac                      <  6:86622/149‑1 (MQ=255)                                                                                                                                tGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGacact                     >  1:263907/1‑147 (MQ=255)                                                                                                                                caattattCCGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACCCg                     >  2:302730/8‑148 (MQ=255)                                                                                                                                  gagaATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGa                   <  7:89715/148‑1 (MQ=255)                                                                                                                                      aTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAAc              <  2:400347/149‑1 (MQ=255)                                                                                                                                      aTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAACCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAAc              <  2:345006/149‑1 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                           <  4:269929/135‑1 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                           >  3:269929/1‑135 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAg                          >  2:212050/1‑136 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAg                          <  1:212050/136‑1 (MQ=255)                                                                                                                                         cGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                                <  4:349422/112‑1 (MQ=255)                                                                                                                                         cGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGAGCAAc              >  4:4783/1‑146 (MQ=255)                                                                                                                                         cGGAGCGT‑GATCGAAGCGCCTGAATAATTGCGGAAATCGTAGCTGTTCTTGAGTTAGATCGCGTCGATCAGCAGTCGTATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTTCGGCAGCCAGACACTGATCAACga            <  2:263907/148‑1 (MQ=255)                                                                                                                                              cTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCATATTCGCTCCAGTTGCGGCAGCCTGACTCGGATCTAGGATAACAg      >  3:323843/1‑149 (MQ=255)                                                                                                                                              cTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAg      <  4:420454/149‑1 (MQ=255)                                                                                                                                                  gATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAgg     <  4:387176/147‑1 (MQ=255)                                                                                                                                                   aTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc    >  3:273066/1‑147 (MQ=255)                                                                                                                                                   aTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc    <  4:273066/147‑1 (MQ=255)                                                                                                                                                   aTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCgg  <  5:128518/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCT‑GCAACGCAG‑‑ACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGG  >  NC_002947/4288880‑4289170

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