Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	162,281	A→T	13.3%	intergenic (‑315/+59)	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,281	0	A	T	13.3%	167.1 / 13.3	66	intergenic (‑315/+59)	PP_0153/scpC	hypothetical protein/propionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base T (7/2);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTG  >  NC_002947/162135‑162429                                                                                                                                                   |                                                                                                                                                      cACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAg                                                                                                                                                     >  3:147303/1‑149 (MQ=255)       tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                <  4:147303/148‑1 (MQ=255)           gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAa‑ga                                                                                                                                          <  5:87922/149‑1 (MQ=255)             gCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAgg                                                                                                                                        >  8:161767/1‑149 (MQ=255)             gCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAgg                                                                                                                                        >  1:76939/1‑149 (MQ=255)              cGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGc                                                                                                                                       <  4:154483/149‑1 (MQ=255)              cGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGc                                                                                                                                       <  1:54657/149‑1 (MQ=255)               gTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGcc                                                                                                                                      >  1:266212/1‑149 (MQ=255)                ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCg                                                                                                                                     <  2:266212/149‑1 (MQ=255)                  ttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCTCTTCGCGGGCTTGCCCGCGAA‑TGGGccat                                                                                                                                    <  2:125545/148‑3 (MQ=39)                     aTTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACa                                                                                                                                <  3:105636/149‑1 (MQ=255)                     aTTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACa                                                                                                                                <  3:152562/149‑1 (MQ=255)                      ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAg                                                                                                                               >  4:39407/1‑149 (MQ=255)                         gaacaaGGGGATTTGTCAGTCGGGAAGCTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCg                                                                                                                             >  2:258827/1‑148 (MQ=255)                           agaagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGaa                                                                                                                          <  7:161767/149‑1 (MQ=255)                                 ggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCaa                                                                                                                                                      >  2:189137/1‑116 (MQ=255)                                 ggATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCaa                                                                                                                                                      <  1:189137/116‑1 (MQ=255)                                  gATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGgc                                                                                                                   >  5:256180/1‑149 (MQ=255)                                       gTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGccgccg                                                                                                              >  8:10037/1‑149 (MQ=255)                                          aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTTGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCAATTCGAGGGCAAGCCCGCGAA‑GAGTCCGGAACAGCGTAAGCAGGGCCGCCGTTc                                                                                                           >  8:144055/1‑149 (MQ=255)                                               ggAAGGTCTGAAAGAACCAGCACAGTAGTACAGAATATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTAc                                                                                                      <  3:39407/149‑1 (MQ=255)                                                gAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTAcc                                                                                                     <  2:76939/149‑1 (MQ=255)                                                       tGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCg                                                                                              >  2:20969/1‑149 (MQ=255)                                                         aaaGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCAATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACGACATCCGGa                                                                                            >  8:224529/1‑149 (MQ=255)                                                             aaCCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAg                                                                                        >  2:113987/1‑149 (MQ=255)                                                               ccAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGG‑TTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTg                                                                                     >  6:223841/1‑149 (MQ=255)                                                                cAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTg                                                                                     >  8:119797/1‑149 (MQ=255)                                                                    acaGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCa                                                                                 <  7:261000/149‑1 (MQ=255)                                                                      agtagtAAAGATTATGCTGCTAGTTGAGGGGTGGCTTTGTGATAGAAGGGTTTGACTGTTATTGCCCTTTCGCGGGCAAGACCTCGAA‑GAGGCCGGAACAGCGTAAGCAGGGCCGCCGTTACGTACCACATCCGGATCAGCTGGc                                                                                   <  7:144055/145‑1 (MQ=255)                                                                           taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATgcg                                                                          >  5:110780/1‑149 (MQ=255)                                                                                aTTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCgg                                                                     >  3:4812/1‑149 (MQ=255)                                                                                 ttATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGt                                                                    <  3:82756/149‑1 (MQ=255)                                                                                     gctgctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttct                                                                >  3:105514/1‑149 (MQ=255)                                                                                      ctgctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttctt                                                               >  2:184775/1‑149 (MQ=255)                                                                                       tgctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                              >  5:127197/1‑149 (MQ=255)                                                                                            gtagAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAg                                                           <  1:184775/147‑1 (MQ=255)                                                                                            gtagAGGGGGGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAg                                                           <  4:68310/147‑1 (MQ=255)                                                                                                agGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCAGTACCAAATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTc                                                     >  5:309026/1‑149 (MQ=255)                                                                                                  gggTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGt                                                   <  6:108818/149‑1 (MQ=255)                                                                                                  gggTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGACCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGt                                                   >  4:95802/1‑149 (MQ=255)                                                                                                    gTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGc                                                 <  1:113987/149‑1 (MQ=255)                                                                                                          tAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCt                                            >  6:144043/1‑148 (MQ=255)                                                                                                            gTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                         <  6:127197/149‑1 (MQ=255)                                                                                                             tGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTcc                                                                                                          >  6:72331/1‑83 (MQ=255)                                                                                                                   aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTGCCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTcg                                  >  4:119065/1‑149 (MQ=255)                                                                                                                   aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCtt                                    >  4:203302/1‑147 (MQ=255)                                                                                                                   aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCtt                                    <  3:203302/147‑1 (MQ=255)                                                                                                                      ggTTTGACTGTTTTGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCgcgc                               <  1:258827/149‑1 (MQ=255)                                                                                                                       gTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCa                              <  7:108955/149‑1 (MQ=255)                                                                                                                            aCTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAt                          <  6:256180/148‑1 (MQ=255)                                                                                                                             cTGTTCCGGCCTCTTCGCGGGCTTGCCCGCGAATGTGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCTTGCGCCCGGTTTCTTCGAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAtg                         >  1:201931/1‑149 (MQ=38)                                                                                                                              tGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGAtgtg                       >  1:5227/1‑149 (MQ=255)                                                                                                                                ttATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCCGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCg                     >  8:238006/1‑149 (MQ=255)                                                                                                                                  aTGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgg                    <  1:20969/148‑1 (MQ=255)                                                                                                                                   tGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                   <  7:10037/148‑1 (MQ=255)                                                                                                                                   tGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                   <  7:119797/148‑1 (MQ=255)                                                                                                                                    ggcctcttCGCGGGCTTGCCCGCGAA‑TAGGCCGGACCAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg                  >  1:46505/7‑148 (MQ=255)                                                                                                                                     gcctcttCGCGGGCTTGCCCGCGAA‑TAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTGGCACATCCGGATCAGCTGGCCAGCATGCGCGCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                >  3:243005/6‑149 (MQ=255)                                                                                                                                     gCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                >  6:79913/1‑149 (MQ=255)                                                                                                                                      cccATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTAt               >  7:98453/1‑149 (MQ=255)                                                                                                                                      cccATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTAt               >  4:242671/1‑149 (MQ=255)                                                                                                                                        tcttCGCGGGCTTGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTACGGGGTATgg             >  6:279036/3‑149 (MQ=255)                                                                                                                                         cttCGCGGGCTCGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGACGTTCCGTACCACATACGGATCAGCTGGCCAGACTGCGCCCGGTTTCTTCAATGTTTTCGTGACAGCTCAGGGCGTCGCGCAGCATGTGAGGGGTATGgc            >  8:278139/2‑149 (MQ=255)                                                                                                                                         aTTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgc            >  1:9684/1‑149 (MQ=255)                                                                                                                                          ttCGCGGGCTTGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgcc           >  8:28079/1‑149 (MQ=255)                                                                                                                                           tCGCGGGCTTGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgccg          >  1:10316/1‑149 (MQ=255)                                                                                                                                            cgcgGGCTTGCCCGCGAA‑TGGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgcc           <  6:226720/147‑1 (MQ=255)                                                                                                                                             gcgGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGcc                                                >  8:14566/1‑109 (MQ=255)                                                                                                                                                   aaGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                   <  7:259309/132‑1 (MQ=255)                                                                                                                                                   aaGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                   >  8:259309/1‑132 (MQ=255)                                                                                                                                                   aaGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTg  <  2:9684/149‑1 (MQ=255)                                                                                                                                                   aaGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTg  <  4:4812/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAA‑GAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTG  >  NC_002947/162135‑162429

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