Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	355,223	T→G	5.6%	intergenic (‑52/+33)	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,223	0	T	G	5.6%	222.6 / 3.5	72	intergenic (‑52/+33)	cbcW/cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base G (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGGC  >  NC_002947/355083‑355369                                                                                                                                             |                                                                                                                                                    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                                                                                                                                          >  6:42590/1‑144 (MQ=255)         gTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAg                                                                                                                                          <  5:42590/144‑1 (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                                                                                                                                             <  1:185475/108‑1 (MQ=255)                                          gaacgaCGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGcccc                                                                                                                                             >  2:185475/1‑108 (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                               <  5:172237/149‑1 (MQ=255)                                                                                                              aTAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTg                               >  2:262400/1‑149 (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                         <  4:139548/148‑1 (MQ=255)                                                                                                                     gAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                         <  1:262400/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)                                                                                                                                           ggTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTggc  <  3:138630/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                    GCGTCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCG‑TAAGCTTGGCCTTGACCGCCTCCAGGCCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGGC  >  NC_002947/355083‑355369

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