Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,289,069	Δ1 bp	10.3%	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,069	0	G	.	10.3%	300.9 / 28.9	97	coding (98/822 nt)	PP_3760	methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base . (5/5);  total (50/47)
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.91e-01

GGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTCCA  >  NC_002947/4288922‑4289218                                                                                                                                                      |                                                                                                                                                      gggTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGaa                                                                                                                                                     <  2:597058/150‑1 (MQ=255)      aCGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                                <  2:1037098/150‑1 (MQ=255)      aCGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGcc                                                                                                                                                <  2:2214582/150‑1 (MQ=255)      aCGGAACATCTCGCTGACCGGGATCGCTTCCTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGcc                                                                                                                                                >  1:1871190/1‑150 (MQ=255)        ggAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCgg                                                                                                                                              >  2:1811533/1‑150 (MQ=255)         gAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                             <  2:1253382/150‑1 (MQ=255)              tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                                >  1:795371/1‑142 (MQ=255)              tctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                                <  2:795371/142‑1 (MQ=255)               ctcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGTCGGAATAAtt                                                                                                                                       >  2:1454863/1‑150 (MQ=255)                tcGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTg                                                                                                                                      >  2:668988/1‑150 (MQ=255)                 cGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGc                                                                                                                                     <  1:1888555/150‑1 (MQ=255)                  gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGGGCATCTACCCTAGCACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                                    <  2:1871190/150‑1 (MQ=255)                  gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                                    <  2:1409268/150‑1 (MQ=255)                  gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCg                                                                                                                                    <  1:1302849/150‑1 (MQ=255)                      aCCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCCAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGaaa                                                                                                                                >  2:37313/1‑150 (MQ=255)                           gATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                                           <  1:1811533/150‑1 (MQ=255)                             tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGc                                                                                                                         >  1:258265/1‑150 (MQ=255)                               gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                             >  2:2408391/1‑128 (MQ=255)                               gTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGCCGGa                                                                                                                                             <  1:2408391/128‑1 (MQ=255)                                tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                      <  2:1539787/150‑1 (MQ=255)                                tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                      >  1:863004/1‑150 (MQ=255)                                acAGATACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGCGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                      >  2:705798/2‑150 (MQ=255)                                aaaGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGt                                                                                                                      >  1:84818/3‑150 (MQ=255)                                      cgtggagcgACTGCATGAACATGCCGGGGTCGTGCAGGACGCGC‑TCCTGCAACGCGGACCCCGTCAGTCAATCGAACTGACGCAGCGCGTGGAGAGTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCGGTACTTGa                                                                                                                <  2:47575/145‑1 (MQ=255)                                      aCTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                                <  1:432363/150‑1 (MQ=255)                                      aCTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                                >  2:462556/1‑150 (MQ=255)                                       cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGA‑TCGTGTAAATCGCACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                               <  2:64372/150‑1 (MQ=255)                                        tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                                <  1:73828/148‑1 (MQ=255)                                        tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGa                                                                                                                >  2:73828/1‑148 (MQ=255)                                        tGCAGAAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAgg                                                                                                              >  2:1722686/1‑150 (MQ=255)                                          cagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTa                                                                                                            >  2:2086801/1‑150 (MQ=255)                                                cTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                                <  1:869277/108‑1 (MQ=255)                                                cTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                                >  2:869277/1‑108 (MQ=255)                                                   cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                                    <  2:2086157/149‑1 (MQ=255)                                                   cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCt                                                                                                    >  1:2086157/1‑149 (MQ=255)                                                   cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                                   >  2:1311578/1‑150 (MQ=255)                                                   cATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                                   >  1:1086988/1‑150 (MQ=255)                                                    gtGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGCTCGCCTCg                                                                                                  >  1:348916/2‑150 (MQ=255)                                                     tGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GGTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                                 >  1:836908/1‑150 (MQ=255)                                                     tGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGa                                                                                                 <  2:912961/150‑1 (MQ=255)                                                      gAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGAt                                                                                                >  1:1606948/1‑150 (MQ=255)                                                       aaCATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATc                                                                                               <  2:258265/150‑1 (MQ=255)                                                              ccGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTc                                                                                        <  2:1606948/150‑1 (MQ=255)                                                               cGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGATAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                         <  1:2268620/100‑1 (MQ=255)                                                               cGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaataa                                                                                                                                         >  2:2268620/1‑100 (MQ=255)                                                                    cGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATGGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                                  <  2:84818/150‑1 (MQ=255)                                                                      tGCAGGA‑GCGCGTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCg                                                                                >  2:607726/1‑150 (MQ=255)                                                                       gCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGa                                                                               >  1:1901974/1‑150 (MQ=255)                                                                            gagcgcgtCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCCCGATGTc                                                                           >  1:1243419/8‑150 (MQ=255)                                                                            aCGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GGTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCg                                                                          <  2:836908/150‑1 (MQ=255)                                                                                 c‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                                     <  2:1086988/150‑1 (MQ=255)                                                                                      tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACtt                                                                                                                  >  1:83381/1‑101 (MQ=255)                                                                                      tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACtt                                                                                                                  <  2:83381/101‑1 (MQ=255)                                                                                      tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                                   <  1:138981/116‑1 (MQ=255)                                                                                      tGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTc                                                                                                   >  2:138981/1‑116 (MQ=255)                                                                                            gCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                                           <  1:462556/150‑1 (MQ=255)                                                                                            gCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                                           <  1:1311578/150‑1 (MQ=255)                                                                                              aGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTAGCTAGt                                                         >  2:1116579/1‑150 (MQ=255)                                                                                                   ccGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                                    <  1:2086801/150‑1 (MQ=255)                                                                                                        aGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCt                                               >  1:2012687/1‑150 (MQ=255)                                                                                                         gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGa                                                                                     <  1:1122601/111‑1 (MQ=255)                                                                                                         gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                              <  1:579632/150‑1 (MQ=255)                                                                                                         gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTc                                              >  1:2361938/1‑150 (MQ=255)                                                                                                         gACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTAGGa                                                                                     >  2:1122601/1‑111 (MQ=255)                                                                                                              tCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                          <  2:1133041/69‑1 (MQ=255)                                                                                                              tCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                                          >  1:1133041/1‑69 (MQ=255)                                                                                                               cGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTg                                        >  1:1693592/1‑150 (MQ=255)                                                                                                                 aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCg                                      >  1:1583187/1‑150 (MQ=255)                                                                                                                 aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCg                                      >  1:2072852/1‑150 (MQ=255)                                                                                                                  aCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCAAGTTGCgg                                     <  2:1901974/150‑1 (MQ=255)                                                                                                                       cGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGcc                                <  2:348916/150‑1 (MQ=255)                                                                                                                        gCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGTGGCAGCTAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCCCCAGTTGCGGCAGCCa                               >  1:150731/1‑150 (MQ=255)                                                                                                                        gCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                               <  2:2012687/150‑1 (MQ=255)                                                                                                                           gcgcGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGac                            >  2:1416789/1‑150 (MQ=255)                                                                                                                             gcgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGAt                                                                              >  2:19076/1‑98 (MQ=255)                                                                                                                             gcgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGAt                                                                              <  1:19076/98‑1 (MQ=255)                                                                                                                              cgTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACg                         <  1:1688057/150‑1 (MQ=255)                                                                                                                                tGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGTTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTACAGTTACGGCCACCAGACACGGa                       >  2:357723/1‑150 (MQ=255)                                                                                                                                 ggAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGAt                      <  2:1693592/150‑1 (MQ=255)                                                                                                                                  gagaATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATc                     >  1:859968/1‑150 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACga                <  1:1416789/150‑1 (MQ=255)                                                                                                                                       tCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACga                <  1:1559345/150‑1 (MQ=255)                                                                                                                                         cGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCCCCAGTTGCGGCAGCCAGACACGGATCAACgaga              >  2:479256/1‑150 (MQ=255)                                                                                                                                          ggCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGaa             >  1:713170/1‑150 (MQ=255)                                                                                                                                          ggCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGaa             <  1:705798/150‑1 (MQ=255)                                                                                                                                              cTT‑GATCGAAGCGTCGGAATAATTGCGGAAATCGTATCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAgg         <  1:1454863/150‑1 (MQ=255)                                                                                                                                              cTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAgg         <  2:2072852/150‑1 (MQ=255)                                                                                                                                               tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                                  <  1:938243/76‑1 (MQ=255)                                                                                                                                               tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGAtct                                                                                  >  2:938243/1‑76 (MQ=255)                                                                                                                                               tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc        >  2:993092/1‑150 (MQ=255)                                                                                                                                               tt‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGc        >  2:492447/1‑150 (MQ=255)                                                                                                                                                  gATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCgg      >  2:763180/1‑150 (MQ=255)                                                                                                                                                  gATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCgg      <  2:1334375/150‑1 (MQ=255)                                                                                                                                                   aTCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGt     >  1:1155099/1‑150 (MQ=255)                                                                                                                                                    tCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                               <  2:527027/123‑1 (MQ=255)                                                                                                                                                    tCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCa                               >  1:527027/1‑123 (MQ=255)                                                                                                                                                    tCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTc    >  1:1503762/1‑150 (MQ=255)                                                                                                                                                      gAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGATCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTCCa  >  2:1983623/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GGGTCACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGC‑TCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTCCA  >  NC_002947/4288922‑4289218

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