Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,765,233	C→G	9.8%	intergenic (‑10/+22)	PP_mr33 ← / ← PP_2419	Cobalamin regulatory region/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,765,233	0	C	G	9.8%	182.1 / 10.1	61	intergenic (‑10/+22)	PP_mr33/PP_2419	Cobalamin regulatory region/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base G (3/3);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATACCGTC  >  NC_002947/2765086‑2765377                                                                                                                                                     |                                                                                                                                                 cTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑cgcg                                                                                                                                                 >  6:234710/1‑149 (MQ=255)      cGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCGCG‑ATGGTCTACCATGCGggccgc                                                                                                                                             >  7:206741/1‑148 (MQ=255)      cGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCGCG‑ATGGTCTACCATGCGggccgc                                                                                                                                             <  8:206741/148‑1 (MQ=255)          cggcAGGCCCAGTGGCATCGATGTAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCAGATTGGCCTGTATTTGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGt                                                                                                                                        <  8:5547/149‑1 (MQ=255)                    aGTGGCATCGTTGGAGCAGGCTCGGCACCTACAGTTCCGGGGGCAGTTCCGATTGGCCTGTATTGGCATGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCtt                                                                                                                              <  2:291106/149‑1 (MQ=255)                         cATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                                         <  7:46957/149‑1 (MQ=255)                          aTCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                                         >  5:153907/1‑148 (MQ=255)                           tCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGaa                                                                                                                       <  2:240723/149‑1 (MQ=255)                                      ggCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGc                                                                                                             >  2:266178/1‑148 (MQ=255)                                               ccTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGTAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGt                                                                                                   <  5:23033/149‑1 (MQ=255)                                                cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     >  1:339983/1‑146 (MQ=255)                                                cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     <  2:339983/146‑1 (MQ=255)                                                cTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTc                                                                                                  >  1:207974/1‑149 (MQ=255)                                                 tACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     >  4:239114/1‑145 (MQ=255)                                                 tACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     <  3:239114/145‑1 (MQ=255)                                                         cGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     >  4:134757/1‑137 (MQ=255)                                                         cGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     <  3:134757/137‑1 (MQ=255)                                                           ggggCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAg                                                                                                                         <  7:221650/115‑1 (MQ=255)                                                             ggCAGTTCGGATTGGCATGTATTGGCAGGCATGGGATTCCCTGTTATTTCCGTGGTGGCCCACTGCGTCCCGCCCGTTAGAACAT‑CGCGCCGCCAGGTGAATACGCTTTTCAGAAGTACTTCAGCCAACTGATGTCGCGGCGACgctgc                                                                                     <  1:283057/149‑1 (MQ=255)                                                             ggCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACgctgc                                                                                     <  2:207974/149‑1 (MQ=255)                                                                 gTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCAGCtt                                                                                   >  7:297213/1‑147 (MQ=255)                                                                  ttCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     <  8:251227/128‑1 (MQ=255)                                                                  ttCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                                     >  7:251227/1‑128 (MQ=255)                                                                      gATTGGCCTGTATTGGCAGGCCTGGGATTGCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCa                                                                                                               <  8:313970/114‑1 (MQ=255)                                                                      gATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgg                                                                            <  5:234710/149‑1 (MQ=255)                                                                      gATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgg                                                                            <  4:189257/149‑1 (MQ=255)                                                                                tATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATg                                                                                                    <  3:43340/115‑1 (MQ=255)                                                                                tATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATg                                                                                                    >  4:43340/1‑115 (MQ=255)                                                                                     gCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTATCTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCAcc                                                             <  1:266178/149‑1 (MQ=255)                                                                                             gggATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAg                                                        >  6:44097/1‑146 (MQ=255)                                                                                             gggATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAg                                                        <  5:44097/146‑1 (MQ=255)                                                                                                aTTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCAGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAg                                                  <  8:297213/149‑1 (MQ=255)                                                                                                 ttCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgc                                                 >  7:43093/1‑149 (MQ=255)                                                                                                    ccTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgcgc                                               <  8:43093/148‑1 (MQ=255)                                                                                                      tGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAc                                            >  3:321144/1‑149 (MQ=255)                                                                                                      tGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAc                                            >  3:75048/1‑149 (MQ=255)                                                                                                        ttAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                                           <  4:321144/148‑1 (MQ=255)                                                                                                              ccGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcag                                      <  6:89551/147‑1 (MQ=255)                                                                                                              ccGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcag                                      >  5:89551/1‑147 (MQ=255)                                                                                                               cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                                           >  4:211898/1‑109 (MQ=255)                                                                                                               cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                                           <  3:211898/109‑1 (MQ=255)                                                                                                               cGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcagcag                                   >  5:230960/1‑149 (MQ=255)                                                                                                                     ggAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAgcagca                                 <  3:304329/145‑1 (MQ=255)                                                                                                                     ggAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACcgc                             >  6:37695/1‑149 (MQ=255)                                                                                                                      gAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACcgc                             <  6:153907/148‑1 (MQ=255)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                               <  7:234698/92‑1 (MQ=25)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                               >  8:234698/3‑94 (MQ=25)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                               >  2:121510/3‑94 (MQ=25)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                               <  1:121510/92‑1 (MQ=25)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                                           <  4:266411/128‑1 (MQ=38)                                                                                                                          ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAAcc                                           >  3:266411/3‑130 (MQ=38)                                                                                                                           cggcggCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg                       >  5:293113/1‑149 (MQ=255)                                                                                                                            ggcggcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGc                      >  8:145662/1‑149 (MQ=255)                                                                                                                             gcggcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGcc                          <  4:224077/144‑1 (MQ=255)                                                                                                                             gcggcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGcc                          >  3:224077/1‑144 (MQ=255)                                                                                                                               ggcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                                           <  2:142469/93‑1 (MQ=255)                                                                                                                               ggcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcggg                                                                           >  1:142469/1‑93 (MQ=255)                                                                                                                                gcCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCcaca                  >  4:265936/1‑149 (MQ=255)                                                                                                                                   cGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg                       <  4:148536/141‑1 (MQ=255)                                                                                                                                   cGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAg                       >  3:148536/1‑141 (MQ=255)                                                                                                                                       tGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGc           <  6:230960/149‑1 (MQ=255)                                                                                                                                          tAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCca                                       <  4:41280/118‑1 (MQ=255)                                                                                                                                          tAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCca                                       >  3:41280/1‑118 (MQ=255)                                                                                                                                                aT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATACCGTc  <  3:265936/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATACCGTC  >  NC_002947/2765086‑2765377

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