Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,289,020:1	+C	7.2%	coding (103/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,020	1	.	C	7.2%	268.1 / 15.3	69	coding (103/822 nt)	PP_3760	methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (28/36);  new base C (3/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

ACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGA  >  NC_002947/4288883‑4289160                                                                                                                                           |                                                                                                                                              aCGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAgcgc                                                                                                                                    >  8:81723/1‑149 (MQ=255)    gAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCgg                                                                                                                                 <  2:120276/149‑1 (MQ=255)        aTCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGaat                                                                                                                              >  6:202309/1‑148 (MQ=255)             gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGc                                                                                                                        <  7:81723/149‑1 (MQ=255)             gCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCAT‑GATCGAAGCGCCGGAATAATTGc                                                                                                                        <  7:208722/149‑1 (MQ=255)               tGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                   <  4:50541/136‑1 (MQ=255)               tGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTTCGATC‑AAGCGcc                                                                                                                                   >  3:50541/1‑136 (MQ=255)                   cGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAAt                                                                                                                  >  2:65336/1‑149 (MQ=255)                      gATCGTCAGGTACTGCGGCAACTGCATGAACACGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGa                                                                                                                                >  5:99356/1‑132 (MQ=255)                      gATCGTCAGGTACTGCGGCAACTGCATGAACACGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGa                                                                                                                                <  6:99356/132‑1 (MQ=255)                       aTCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTa                                                                                                              <  2:61657/149‑1 (MQ=255)                        tCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAg                                                                                                             <  1:204406/149‑1 (MQ=255)                           tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTg                                                                                                          >  6:165698/1‑149 (MQ=255)                           tCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTg                                                                                                          >  3:324501/1‑149 (MQ=255)                              ggTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTAc                                                                                                       <  2:326001/149‑1 (MQ=255)                                  cTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACtt                                                                                                     >  5:108787/1‑147 (MQ=255)                                   tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                  >  3:38485/1‑149 (MQ=255)                                   tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                  >  2:306742/1‑149 (MQ=255)                                   tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACACAGCGCGTGGAGAATCCGGCGCTT‑GATCGaa                                                                                                                                        >  7:155131/1‑111 (MQ=255)                                   tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACACAGCGCGTGGAGAATCCGGCGCTT‑GATCGaa                                                                                                                                        <  8:155131/111‑1 (MQ=255)                                     cagcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGTGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                  <  5:81976/147‑1 (MQ=255)                                       gcaACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAgg                                                                                                 >  2:101576/1‑146 (MQ=255)                                              aaTGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGAAGACACCGTCAGACAATAGAACTGACGCAGCGCGTGGAGAATCCTGCGCTT‑GATCGAAGCGCCGGAAAAATAGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGcat                                                                                      >  7:62187/2‑147 (MQ=255)                                               aTGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTc                                                                                     <  8:13793/149‑1 (MQ=255)                                                 gAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGa                                                                                   >  7:40897/1‑149 (MQ=255)                                                     aTGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATcag                                                                               <  4:240997/149‑1 (MQ=255)                                                       gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATc                                                                                 >  4:287222/1‑145 (MQ=255)                                                       gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATc                                                                                 <  3:287222/145‑1 (MQ=255)                                                           gggTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCg                                                                         >  7:42095/1‑149 (MQ=255)                                                            ggTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCgg                                                                        <  2:266364/149‑1 (MQ=255)                                                                   cAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGa                                                                 >  5:197292/1‑149 (MQ=255)                                                                       aCGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTc                                                             <  1:304813/149‑1 (MQ=255)                                                                          cgcTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGgt                                                          <  8:175828/149‑1 (MQ=255)                                                                          cgcTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGgt                                                          <  5:165698/149‑1 (MQ=255)                                                                           gcTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGgtg                                                         >  5:267384/1‑149 (MQ=255)                                                                            cTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACTCATCGCGTGTAGAATCCGGCTCTT‑GATCGAAGCGCCTGAATAATTGCGGAAATCGTTGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCTGTCGTATCTCGATGTCGgtgt                                                        <  6:64688/149‑1 (MQ=255)                                                                             tCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                       >  7:262944/1‑149 (MQ=255)                                                                              ccTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GAACCAAGCGCCGGAATAATTGCGGAAATCGTAGAGGTACATGAGCTAGAT‑CGCCTAGATCAGCAGTCGGATCTCGATGAAGGTGTTg                                                      >  2:158521/1‑149 (MQ=255)                                                                               cTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                       <  8:40897/147‑1 (MQ=255)                                                                               cTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                       <  1:101576/147‑1 (MQ=255)                                                                                  cAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGtt                                                  >  2:109212/1‑149 (MQ=255)                                                                                    aCGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCg                                                <  6:108787/149‑1 (MQ=255)                                                                                     cGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGATCCGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCg                                                <  1:65336/149‑1 (MQ=255)                                                                                     cGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGAt                                                                                  >  8:241943/1‑114 (MQ=255)                                                                                     cGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGAt                                                                                  <  7:241943/114‑1 (MQ=255)                                                                                       cAGACACCGTCAGACAATCGAACTGACTCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                             <  1:306742/149‑1 (MQ=255)                                                                                        aGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAg                                            <  5:202309/149‑1 (MQ=255)                                                                                         gACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAg                                            >  4:139818/1‑148 (MQ=255)                                                                                            acCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGTTAGTCaa                                        <  4:38485/149‑1 (MQ=255)                                                                                              cGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATATTTGCGGAAATCGTAGCTGTACTTGAGGTATAT‑CGGCTCGATCAGGAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAAt                                      <  1:158521/149‑1 (MQ=255)                                                                                                  aGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGc                                  >  6:161457/1‑149 (MQ=255)                                                                                                      aaTCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                             <  1:242125/134‑1 (MQ=255)                                                                                                      aaTCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTa                                             >  2:242125/1‑134 (MQ=255)                                                                                                      aaTCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTc                                          >  3:277753/1‑137 (MQ=255)                                                                                                      aaTCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTc                                          <  4:277753/137‑1 (MQ=255)                                                                                                           aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGAtctc                                                                   >  8:157881/1‑107 (MQ=255)                                                                                                           aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGAtctc                                                                   <  7:157881/107‑1 (MQ=255)                                                                                                           aaCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGc                         <  7:59145/149‑1 (MQ=255)                                                                                                            aCTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCg                        <  8:262944/149‑1 (MQ=255)                                                                                                                 cGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGc                   <  6:197292/149‑1 (MQ=255)                                                                                                                 cGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGc                   <  8:42095/149‑1 (MQ=255)                                                                                                                        cgTGGAGAATCAGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGacac            <  4:324501/149‑1 (MQ=255)                                                                                                                                aTCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAAc    >  5:335557/1‑149 (MQ=255)                                                                                                                                   cGGCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAAc    >  1:25299/1‑146 (MQ=255)                                                                                                                                   cGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACga  <  3:139818/148‑1 (MQ=255)                                                                                                                                    ggCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGc                   <  3:78737/130‑1 (MQ=255)                                                                                                                                    ggCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGc                   >  4:78737/1‑130 (MQ=255)                                                                                                                                       gcTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCGGACACGGATCAAc    <  6:215023/142‑1 (MQ=255)                                                                                                                                       gcTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCGGACACGGATCAAc    >  5:215023/1‑142 (MQ=255)                                                                                                                                           |                                                                                                                                              ACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGA  >  NC_002947/4288883‑4289160

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