Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,765,224	A→C	12.4%	intergenic (‑1/+31)	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,224	0	A	C	12.4%	175.9 / 14.7	64	intergenic (‑1/+31)	PP_mr33/PP_2419	Cobalamin regulatory region/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base C (4/4);  total (30/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 = 9.96e-01

GGCCTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGC  >  NC_002947/2765083‑2765368                                                                                                                                              |                                                                                                                                                  ggCCTTCCCGGGCGTCAGGCACAGTGGCATCGATGGAGAAGGCTCGGCACCTACAGTTGCGGGGGCAGTTACGATTGGCCTCTATTGGCAGGCATGGGATTCCCTGGTAGTTCCGTGGTTGAAACCGGCGGCCCGAATGGTAGACCAT‑c                                                                                                                                           >  2:283057/1‑149 (MQ=255)    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                                                                                            <  3:239114/145‑1 (MQ=255)                                                    tACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGAt                                                                                            >  4:239114/1‑145 (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                                                                   <  4:189257/149‑1 (MQ=255)                                                                         gATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgg                                                                   <  5:234710/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                             >  5:89551/1‑147 (MQ=255)                                                                                                                 ccGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCcag                             <  6:89551/147‑1 (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                                                                      >  8:234698/3‑94 (MQ=25)                                                                                                                             ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                      >  2:121510/3‑94 (MQ=25)                                                                                                                             ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                      <  1:121510/92‑1 (MQ=25)                                                                                                                             ctggcggcGCG‑ATGGTCTACCATGCGGGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAgg                                                                      <  7:234698/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)                                                                                                                                              |                                                                                                                                                  GGCCTTCCCGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCAT‑CGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGC  >  NC_002947/2765083‑2765368

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