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

ACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTCC  >  NC_002947/4288883‑4289173                                                                                                                                                |                                                                                                                                                      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                                                                                                               >  2:306742/1‑149 (MQ=255)                                   tGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAg                                                                                                               >  3:38485/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                                                                                              <  3:287222/145‑1 (MQ=255)                                                       gCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATc                                                                                              >  4:287222/1‑145 (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                                                                    <  1:101576/147‑1 (MQ=255)                                                                               cTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGtt                                                                    <  8:40897/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                                >  4:78737/1‑130 (MQ=255)                                                                                                                                    ggCGCTTCGATC‑AAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGc                                <  3:78737/130‑1 (MQ=255)                                                                                                                                       gcTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCGGACACGGATCAAc                 <  6:215023/142‑1 (MQ=255)                                                                                                                                       gcTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCGGACACGGATCAAc                 >  5:215023/1‑142 (MQ=255)                                                                                                                                                gAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCCATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTcc  <  1:109212/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      ACGGAACATCTCGCTGACCGGGATCGTCAGGTACTGCAGCAACTGCATGAACATGCCCGGGTCGTGCAGGACGCGCTCCTGCAACGCAGACACCGTCAGACAATCGAACTGACGCAGCGCGTGGAGAATCCGGCGCTT‑GATCGAAGCGCCGGAATAATTGCGGAAATCGTAGCTGTACTTGAGGTAGAT‑CGCCTCGATCAGCAGTCGGATCTCGATGTCGGTGTTGCGTTCGCTAGTCAAATTCGCTCCAGTTGCGGCAGCCAGACACGGATCAACGAGAACAGGCGGTCC  >  NC_002947/4288883‑4289173

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