Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	355,218	C→A	5.7%	intergenic (‑47/+38)	cbcW ← / ← cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	355,218	0	C	A	5.7%	224.7 / 3.6	70	intergenic (‑47/+38)	cbcW/cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base A (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

CTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCT  >  NC_002947/355072‑355365                                                                                                                                                   |                                                                                                                                                     cTTGCGCGATGGCGTCGAAGTAATTCGCACCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCAc                                                                                                                                                    >  1:199572/1‑149 (MQ=255)            gCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGcccc                                                                                                                                         <  2:69709/149‑1 (MQ=255)              gTCGATGTAATTCGCGCCGTTCTGGCTCAGCCATTCGAAGAACGACGCTATGTACTGGCCCAGGGGTATTTTCTGATCTATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGa                                                                                                                                       <  4:25512/149‑1 (MQ=255)              gTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGcccc                                                                                                                                         >  7:288847/1‑147 (MQ=255)                 gAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCGGCCCCGAgg                                                                                                                                     >  7:304756/1‑148 (MQ=255)                   aGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTa                                                                                                                                   <  8:223676/148‑1 (MQ=255)                    gTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                      <  5:42590/144‑1 (MQ=255)                    gTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                      >  6:42590/1‑144 (MQ=255)                        ttCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTa                                                                                                                             >  5:67111/1‑149 (MQ=255)                           gcgcCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTg                                                                                                                          <  8:250183/149‑1 (MQ=255)                              ccGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGCCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGt                                                                                                                       <  5:311063/149‑1 (MQ=255)                               cGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTa                                                                                                                      <  5:225279/149‑1 (MQ=255)                                gTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAg                                                                                                                     >  4:314563/1‑149 (MQ=255)                                 ttCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTGTGATCGATAAGCATGATAGCGAGCTTCGACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGt                                                                                                                    >  8:239270/1‑149 (MQ=255)                                      ggTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTc                                                                                                                  <  3:314563/146‑1 (MQ=255)                                      ggTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCCGGGCAGGACCGTGCCTGCCCCgggtt                                                                                                                                    <  2:18040/128‑5 (MQ=255)                                      ggTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCCGGGCAGGACCGTGCCTGCCCCgggtt                                                                                                                                    >  1:18040/1‑124 (MQ=255)                                              aTTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGt                                                                                                                    <  8:316635/136‑1 (MQ=255)                                              aTTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGt                                                                                                                    >  7:316635/1‑136 (MQ=255)                                                  gacgaACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGc                                                                                                   <  1:53450/148‑1 (MQ=255)                                                     gaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGcccc                                                                                                                                         >  2:185475/1‑108 (MQ=255)                                                     gaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGcccc                                                                                                                                         <  1:185475/108‑1 (MQ=255)                                                     gaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTgg                                                                                               >  3:92915/1‑149 (MQ=255)                                                         gacgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCtt                                                                                           >  7:215034/1‑149 (MQ=255)                                                          acgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCtt                                                                                           >  3:254774/1‑148 (MQ=255)                                                          acgCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATAGAAATCAACACGGGCCAGGCACGGTCCTGCCCCGACGTAGCCCTACTGCGTATTTCGTTATTGCG‑TAAGCTAGGCCtt                                                                                           >  3:293228/1‑148 (MQ=255)                                                             cgATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGAcc                                                                                       <  6:296363/149‑1 (MQ=255)                                                                   aCTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTc                                                                                  >  2:249764/1‑148 (MQ=255)                                                                       gCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGc                                                                                                                         <  2:233225/106‑1 (MQ=255)                                                                       gCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGc                                                                                                                         >  1:233225/1‑106 (MQ=255)                                                                       gCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTccagg                                                                              <  1:12863/148‑1 (MQ=255)                                                                       gCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTccagg                                                                              <  1:166082/148‑1 (MQ=255)                                                                         ccAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTccaggc                                                                             >  4:352287/1‑147 (MQ=255)                                                                         ccAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTccaggc                                                                             >  8:291218/1‑147 (MQ=255)                                                                          cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCt                                                                                   <  5:121634/140‑1 (MQ=255)                                                                          cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCt                                                                                   >  6:121634/1‑140 (MQ=255)                                                                          cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTccaggccag                                                                          >  4:138630/1‑149 (MQ=255)                                                                             gggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                           >  5:176176/1‑114 (MQ=255)                                                                             gggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                           <  6:176176/114‑1 (MQ=255)                                                                              ggTAATTTCTGCACGATACGCATGATAGCGAGCTTCCACCTGCATAGATTGAACTCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCtt                                                                      <  8:304756/149‑1 (MQ=255)                                                                               gTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTg                                                                     >  2:350687/1‑149 (MQ=255)                                                                                  ttttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGAGATACTGCGTAGTTCCTTATTGCG‑TAAGCTAGGCCTTGACCGCCTCCAGGCCAGGCTTGCCa                                                                  >  5:197579/1‑149 (MQ=255)                                                                                    ttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCAt                                                                 <  7:291218/148‑1 (MQ=255)                                                                                      cTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCa                                                               <  3:352287/148‑1 (MQ=255)                                                                                         aTCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCTGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCTTAGTTCGTTATTGCG‑TAAGCTTGGACTTGACCGCCTCCAGGCCAGGCTTGCCATCAAc                                                             <  4:92915/147‑1 (MQ=255)                                                                                         aTCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGATAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAAc                                                             <  1:249764/148‑1 (MQ=255)                                                                                            gATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACggtg                                                         <  4:254774/148‑1 (MQ=255)                                                                                              tAAGCATGATAGCGAGCTTCCACCTTCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTa                                                      <  8:215034/149‑1 (MQ=255)                                                                                                  cATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGc                                                  >  3:139548/1‑149 (MQ=255)                                                                                                    tGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCg                                                >  7:104502/1‑149 (MQ=255)                                                                                                     gATAGCGAGCTTCCACCTGAATAGTTTGAAATCAACCCGGGGCAGGCTCTGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCATCGGTGGTAACGCCg                                                <  6:197579/148‑1 (MQ=255)                                                                                                          cGAGCTTCCACCTGCCTAGATTGAAATCAACCCGCGTCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGgccagc                                          <  8:288847/149‑1 (MQ=255)                                                                                                          cGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGgccagc                                          >  3:344059/1‑149 (MQ=255)                                                                                                          cGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGgccagc                                          >  3:338558/1‑149 (MQ=255)                                                                                                                     cTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAg                               <  4:105049/149‑1 (MQ=255)                                                                                                                         aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTg                           >  6:172237/1‑149 (MQ=255)                                                                                                                         aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTg                           >  2:262400/1‑149 (MQ=255)                                                                                                                         aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTg                           <  5:172237/149‑1 (MQ=255)                                                                                                                         aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTg                           >  4:286860/1‑149 (MQ=255)                                                                                                                              ttGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                      >  6:120663/1‑149 (MQ=255)                                                                                                                              ttGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                      >  7:160276/1‑149 (MQ=255)                                                                                                                                gAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                     <  1:262400/148‑1 (MQ=255)                                                                                                                                gAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                     <  4:139548/148‑1 (MQ=255)                                                                                                                                        cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTg            <  5:120663/149‑1 (MQ=255)                                                                                                                                         ccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa           <  2:199572/149‑1 (MQ=255)                                                                                                                                         ccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa           >  8:52113/1‑149 (MQ=255)                                                                                                                                           ggggCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGccagccag                                       <  6:347372/119‑1 (MQ=255)                                                                                                                                           ggggCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGccagccag                                       >  5:347372/1‑119 (MQ=255)                                                                                                                                           ggggCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATc         <  7:239270/149‑1 (MQ=255)                                                                                                                                           ggggCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATc         >  4:77667/1‑149 (MQ=255)                                                                                                                                                aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCAttt    >  4:114770/1‑149 (MQ=255)                                                                                                                                                  gcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCt  >  6:364836/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTTGCGCGATGGCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCT  >  NC_002947/355072‑355365

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