Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	162,282	A→T	5.7%	intergenic (‑316/+58)	PP_0153 ← / ← scpC	hypothetical protein/propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	162,282	0	A	T	5.7%	303.0 / 4.8	98	intergenic (‑316/+58)	PP_0153/scpC	hypothetical protein/propionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (50/41);  new base T (3/3);  total (54/44)
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.20e-01

ACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCT  >  NC_002947/162136‑162428                                                                                                                                                     |                                                                                                                                                      aCTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGc                                                                                                                                                     >  3:367295/1‑149 (MQ=255) aCTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGc                                                                                                                                                     <  2:410685/149‑1 (MQ=255)     aTTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 >  6:252322/1‑149 (MQ=255)      tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 <  1:135042/148‑1 (MQ=255)      tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 >  5:308074/1‑148 (MQ=255)      tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 >  1:27038/1‑148 (MQ=255)         ggTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTTGTAGCGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  5:302628/148‑1 (MQ=255)         ggTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgcg                                                                                                                                                >  1:288384/1‑146 (MQ=255)          gagCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAag                                                                                                                                             <  1:430013/146‑1 (MQ=255)          gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCTTGCCCGCGaat                                                                                                                                             <  5:97764/148‑2 (MQ=255)          gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAaga                                                                                                                                            <  4:367295/149‑1 (MQ=255)              gTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGc‑‑c                                                                                                                                      <  5:32558/149‑1 (MQ=255)               ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAACAGGC‑‑Cg                                                                                                                                     >  2:349729/1‑149 (MQ=255)                tttGATTCTAAGAATTGGATTTGTCAGTCGGGAAGTGCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑Cg                                                                                                                                     <  6:428220/148‑1 (MQ=255)                 ttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGGTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑Cgg                                                                                                                                    <  5:368360/148‑1 (MQ=255)                  tGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑Cgg                                                                                                                                    >  5:26219/1‑147 (MQ=255)                  tGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCACGGGCAAGCCCGCGAAGAGGC‑‑Cgg                                                                                                                                    <  6:26219/147‑1 (MQ=255)                    aTTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGaa                                                                                                                                  >  5:280239/1‑147 (MQ=255)                     ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAg                                                                                                                               <  2:80927/149‑1 (MQ=255)                        gaagaaGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑g                                                                                                                            >  2:217946/1‑148 (MQ=255)                         aagaagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑gg                                                                                                                           >  4:295693/1‑148 (MQ=255)                          agaagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGaa                                                                                                                         <  6:77071/149‑1 (MQ=255)                          agaagGGGATTTGTCAGTAGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTCCGGCCTCTTCGCGGGCTTGCCCGCGAATGGGC‑‑CATAACAgtcaaa                                                                                                                          <  8:89505/149‑6 (MQ=25)                               gggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGAAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAgg                                                                                                                    >  3:291080/1‑149 (MQ=255)                                ggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAggg                                                                                                                   >  4:121881/1‑149 (MQ=255)                                 gATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTGGTACAGATTATGCTGCTAGTAGAGGGAGGTGGCTTAGTTATAAAAGGTTTTCAGTTTTATGCCAATTTCGCGGAAAAGCCCGCTAAGAGCC‑‑CGAAACAGC‑GAACGCAgg                                                                                                                    >  1:130291/1‑149 (MQ=255)                                 gATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGATGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATCGCACATTCGCGGGCACGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAgtgt                                                                                                                  >  6:21218/1‑146 (MQ=255)                                     tGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGgccgcc                                                                                                              >  2:359952/1‑149 (MQ=255)                                      gTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGccgccg                                                                                                             >  7:262036/1‑149 (MQ=255)                                       tCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGt                                                                                                            >  4:175076/1‑149 (MQ=255)                                         aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGtt                                                                                                           <  5:435669/148‑1 (MQ=255)                                         aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTc                                                                                                          >  5:151002/1‑149 (MQ=255)                                         aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCG‑AGAGGC‑‑CGGAAAAGG‑CGGAGCAGGGCCGCCGttgc                                                                                                         >  3:209240/1‑147 (MQ=255)                                             gggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTa                                                                                                      >  5:101742/1‑149 (MQ=255)                                               gAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTAc                                                                                                     >  3:70559/1‑148 (MQ=255)                                                   gTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACAt                                                                                                >  5:99444/1‑149 (MQ=255)                                                    tCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATc                                                                                               <  3:121881/149‑1 (MQ=255)                                                       gAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCg                                                                                             >  5:148244/1‑148 (MQ=255)                                                        aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGtt                                                                                                           <  8:113338/133‑1 (MQ=255)                                                        aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGtt                                                                                                           >  7:113338/1‑133 (MQ=255)                                                         aaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGAt                                                                                          >  7:335894/1‑149 (MQ=255)                                                               cAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCg                                                                                                        <  7:409808/129‑1 (MQ=255)                                                               cAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCg                                                                                                        >  8:409808/1‑129 (MQ=255)                                                                 gCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                                  >  3:363831/1‑149 (MQ=255)                                                                 gCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                                  <  5:252322/149‑1 (MQ=255)                                                                  cacaGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                                  >  1:162256/1‑148 (MQ=255)                                                                          taCAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGACAGCATgcg                                                                         >  5:211550/1‑149 (MQ=255)                                                                            cAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGcc                                                                       <  8:154421/149‑1 (MQ=255)                                                                              gATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCg                                                                     <  6:211550/149‑1 (MQ=255)                                                                                ttATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGt                                                                   <  8:262036/149‑1 (MQ=255)                                                                                     ctgctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttct                                                               <  4:77005/148‑1 (MQ=255)                                                                                      tgctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                             >  6:78548/1‑149 (MQ=255)                                                                                      tgctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                             >  4:131236/1‑149 (MQ=255)                                                                                       gctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                             >  1:116689/1‑148 (MQ=255)                                                                                        ctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                             <  6:280239/147‑1 (MQ=255)                                                                                        ctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCa                                                           <  3:152241/149‑1 (MQ=255)                                                                                        ctAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATATTCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTcc                                                            <  1:359952/148‑1 (MQ=255)                                                                                          agtagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAgg                                                                                                                    <  3:67251/90‑1 (MQ=255)                                                                                          agtagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAgg                                                                                                                    >  4:67251/1‑90 (MQ=255)                                                                                          agtagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGt                                                                                                       >  8:234038/1‑103 (MQ=255)                                                                                          agtagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGt                                                                                                       <  7:234038/103‑1 (MQ=255)                                                                                            tagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGtt                                                       >  8:66885/1‑149 (MQ=255)                                                                                            tagAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGtt                                                       >  7:173772/1‑149 (MQ=255)                                                                                              gagG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCTAGGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGtt                                                       >  1:137024/1‑149 (MQ=255)                                                                                                 g‑‑ggTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGt                                                  >  2:307142/1‑149 (MQ=255)                                                                                                        gCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAg                                             <  6:101742/149‑1 (MQ=255)                                                                                                          gTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGc                                            >  2:86912/1‑148 (MQ=255)                                                                                                             gTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                        >  6:347182/1‑149 (MQ=255)                                                                                                             gTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                        <  6:308074/149‑1 (MQ=255)                                                                                                             gTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                        <  2:288384/149‑1 (MQ=255)                                                                                                             gTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                        <  3:175076/149‑1 (MQ=255)                                                                                                              tGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGCAGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                        >  1:130736/1‑149 (MQ=255)                                                                                                                aTAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGTGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGc                                     >  2:121264/1‑149 (MQ=255)                                                                                                                aTAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAggg                                      <  8:173772/148‑1 (MQ=255)                                                                                                                aTAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAggg                                      <  4:291080/148‑1 (MQ=255)                                                                                                                 tAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCt                                    >  3:310886/1‑149 (MQ=255)                                                                                                                   gAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCtt                                   <  4:363831/148‑1 (MQ=255)                                                                                                                          ttGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAg                            >  1:317622/1‑148 (MQ=255)                                                                                                                            gACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAt                         >  3:71357/1‑149 (MQ=255)                                                                                                                               tGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGCAtgtg                      >  3:173521/1‑149 (MQ=255)                                                                                                                                 ttATGGCCCATTGCCGGGCATGCCGGC‑TAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGc                     <  4:209240/147‑1 (MQ=255)                                                                                                                                 ttATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCg                    >  6:223868/1‑149 (MQ=255)                                                                                                                                   aTGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgg                   <  6:148244/148‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg                 <  4:182619/149‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                  <  4:71357/148‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                  <  3:131236/148‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                  <  1:217946/148‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                  <  1:86912/148‑1 (MQ=255)                                                                                                                                    tGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                  <  2:162256/148‑1 (MQ=255)                                                                                                                                     ggcctcttCGCGGGCTTGCCCGCGAATGGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg                 >  8:266511/7‑148 (MQ=255)                                                                                                                                     ggCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCCCATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGt                >  1:411310/1‑149 (MQ=255)                                                                                                                                      gCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCACCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGt                <  8:267785/148‑1 (MQ=255)                                                                                                                                        ccATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg             >  3:347374/1‑149 (MQ=255)                                                                                                                                        ccATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg             >  6:40730/1‑149 (MQ=255)                                                                                                                                        ccATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg             <  5:223868/149‑1 (MQ=255)                                                                                                                                             cgcgGGCTTGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgccgc        >  2:232610/1‑149 (MQ=255)                                                                                                                                               cgGGCTTGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCa      >  4:182037/1‑149 (MQ=255)                                                                                                                                                   cAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCt  >  7:158531/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGG‑‑GGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGC‑‑CGGAACAGC‑GGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCT  >  NC_002947/162136‑162428

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