Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 162,282 A→T 7.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_002947162,2820AT 7.7% 644.0 / 33.2 207intergenic (‑316/+58)PP_0153/scpChypothetical protein/propionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (118/73);  new base T (10/6);  total (128/79)
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.14e-01

Only 100 of 218 total aligned reads displayed.

CGCACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTGGC  >  NC_002947/162133‑162431
                                                                                                                                                     |                                                                                                                                                     
 gCACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAg                                                                                                                                                      >  1:2815542/1‑150 (MQ=255)
    cTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGccc                                                                                                                                                   >  1:3523578/1‑150 (MQ=255)
        tttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGa                                                                                                                                               >  1:5095411/1‑150 (MQ=255)
         ttGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  1:2082673/150‑1 (MQ=255)
           ggTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAaga                                                                                                                                            >  1:598961/1‑150 (MQ=255)
           ggTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAaga                                                                                                                                            >  2:3414401/1‑150 (MQ=255)
            gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCTCTTCGCGGGCTTGCCCGCGaatgg                                                                                                                                           <  1:3079201/150‑4 (MQ=255)
            gTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCTTGCCCGCGaatgg                                                                                                                                           <  1:545844/150‑4 (MQ=255)
             tGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAgg                                                                                                                                          <  1:243386/150‑1 (MQ=255)
                 ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCgg                                                                                                                                      <  1:5309693/150‑1 (MQ=255)
                 ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCgg                                                                                                                                      >  1:317647/1‑150 (MQ=255)
                 ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCgg                                                                                                                                      >  2:1103387/1‑150 (MQ=255)
                 ttttGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCgg                                                                                                                                      >  1:2283175/1‑150 (MQ=255)
                       ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 >  1:2464353/1‑133 (MQ=255)
                       ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCcgc                                                                                                                                                 <  2:2464353/133‑1 (MQ=255)
                       ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAaga                                                                                                                                            <  2:5008917/138‑1 (MQ=255)
                       ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                                >  1:1920706/1‑150 (MQ=255)
                       ttCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGc                                                                                                                                >  1:1790374/1‑150 (MQ=255)
                        tCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCg                                                                                                                               >  1:4552317/1‑150 (MQ=255)
                         cGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCgg                                                                                                                              <  2:1060113/150‑1 (MQ=255)
                          gaagaaGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTTGTGATAGACGGGTTGGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGa                                                                                                                             <  1:706045/150‑1 (MQ=255)
                          gaagaaGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGa                                                                                                                             >  1:4479346/1‑150 (MQ=255)
                           aagaagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGaa                                                                                                                            >  2:490367/1‑150 (MQ=255)
                              aagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCa                                                                                                                         >  2:2732608/1‑150 (MQ=255)
                              aagGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCa                                                                                                                         >  2:2631973/1‑150 (MQ=255)
                                   gATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGgcc                                                                                                                    >  2:2138996/1‑150 (MQ=255)
                                          cAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGaa                                                                                                                                              >  2:4402321/1‑117 (MQ=255)
                                          cAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  1:4402321/117‑1 (MQ=255)
                                           aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTcc                                                                                                            >  2:1584575/1‑150 (MQ=255)
                                           aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTcc                                                                                                            >  1:3971078/1‑150 (MQ=255)
                                           aGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTcc                                                                                                            >  1:2981447/1‑150 (MQ=255)
                                               gggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCTCTTCGCGGGCTTGCCCGCGaatgg                                                                                                                                           >  1:4295828/1‑112 (MQ=39)
                                               gggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCTCTTCGCGGGCTTGCCCGCGaatgg                                                                                                                                           <  2:4295828/115‑4 (MQ=39)
                                               gggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTAc                                                                                                        <  2:3133425/150‑1 (MQ=255)
                                               gggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTAc                                                                                                        <  2:317647/150‑1 (MQ=255)
                                                ggAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTAcc                                                                                                       >  1:3282190/1‑150 (MQ=255)
                                                 gAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACca                                                                                                      >  2:2377687/1‑150 (MQ=255)
                                                 gAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACca                                                                                                      >  1:1604235/1‑150 (MQ=255)
                                                     gTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATc                                                                                                  >  1:1588215/1‑150 (MQ=255)
                                                      tCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATcc                                                                                                 >  2:3138075/1‑150 (MQ=255)
                                                             gAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAg                                                                                          >  2:3052821/1‑150 (MQ=255)
                                                              aaCCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGc                                                                                         <  2:5095608/150‑1 (MQ=255)
                                                                          agtaCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATgc                                                                             <  1:4002890/150‑1 (MQ=255)
                                                                           gtaCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCAGGCCAGCATgcg                                                                            >  1:4564946/1‑150 (MQ=255)
                                                                            taCAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATgcgc                                                                           >  1:5237422/1‑150 (MQ=255)
                                                                                  ttATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGtt                                                                     <  2:759298/150‑1 (MQ=255)
                                                                                       ctgctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTttcttc                                                                <  1:2188848/150‑1 (MQ=255)
                                                                                         gctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCa                                                              <  2:2806862/150‑1 (MQ=255)
                                                                                         gctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCa                                                              <  2:4693191/150‑1 (MQ=255)
                                                                                         gctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCa                                                              <  2:2109652/150‑1 (MQ=255)
                                                                                         gctAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCa                                                              <  2:177862/150‑1 (MQ=255)
                                                                                           tagtagAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAgg                                                            >  1:4871364/1‑150 (MQ=255)
                                                                                              tagAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                                       >  1:3739485/1‑88 (MQ=255)
                                                                                              tagAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAgg                                                                                                                       <  2:3739485/88‑1 (MQ=255)
                                                                                                gagGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTCTGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTc                                                       >  1:4934841/1‑150 (MQ=255)
                                                                                                gagGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTc                                                       >  2:3578733/1‑150 (MQ=255)
                                                                                                    ggTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGc                                                   <  2:3516811/150‑1 (MQ=255)
                                                                                                        gCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGTAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGc                                               >  1:1370272/1‑150 (MQ=255)
                                                                                                        gCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGTAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGc                                               >  1:1987014/1‑150 (MQ=255)
                                                                                                          gTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGc                                                   >  2:447715/1‑144 (MQ=255)
                                                                                                            aGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                           >  1:928879/1‑150 (MQ=255)
                                                                                                            aGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAg                                           <  1:1526800/150‑1 (MQ=255)
                                                                                                                 tAGAAGGGTTGGACTGTTCTGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCtt                                      <  2:5143941/150‑1 (MQ=255)
                                                                                                                  aGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTc                                     >  2:3969757/1‑150 (MQ=255)
                                                                                                                    aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGc                                                                                 >  1:4794502/1‑104 (MQ=255)
                                                                                                                    aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTcgc                                   >  2:5186367/1‑150 (MQ=255)
                                                                                                                    aaGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGAACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTcgc                                   >  1:4069640/1‑150 (MQ=255)
                                                                                                                      gggTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCgcgc                                 <  1:3052887/150‑1 (MQ=255)
                                                                                                                          ttGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGa                             <  2:3630177/150‑1 (MQ=255)
                                                                                                                               tGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGc                        >  2:451449/1‑150 (MQ=255)
                                                                                                                                 ttATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgg                      >  2:3816185/1‑150 (MQ=255)
                                                                                                                                 ttATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgg                      >  1:4732656/1‑150 (MQ=255)
                                                                                                                                 ttATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGAACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgg                      <  2:4069640/150‑1 (MQ=255)
                                                                                                                                  tATGGCCCCTTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                     <  2:5272870/150‑1 (MQ=255)
                                                                                                                                  tATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                     <  2:701831/150‑1 (MQ=255)
                                                                                                                                  tATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCggg                     <  1:1074287/150‑1 (MQ=255)
                                                                                                                                   aTGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg                    <  1:1849304/150‑1 (MQ=255)
                                                                                                                                   aTGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCgggg                    <  1:875051/150‑1 (MQ=255)
                                                                                                                                     ggcctcttCGCGGGCTTGCCCGCGAATGGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                  >  2:4749997/7‑150 (MQ=255)
                                                                                                                                     ggcctcttCGCGGGCTTGCCCGCGAATGGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                  >  1:1839495/7‑150 (MQ=255)
                                                                                                                                     ggcctcttCGCGGGCTTGCCCGCGAATGGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                  >  2:2690627/7‑150 (MQ=255)
                                                                                                                                     ggCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTa                  <  1:3138075/150‑1 (MQ=255)
                                                                                                                                      gCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTAt                 >  2:872494/1‑150 (MQ=255)
                                                                                                                                       cccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGcc                                                  <  1:1489663/116‑1 (MQ=255)
                                                                                                                                       cccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATg                >  1:2612212/1‑150 (MQ=255)
                                                                                                                                       cccATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCAGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCCGGGCTTCGCGCAGGATGTGCGGGGTATg                >  1:3002261/1‑150 (MQ=255)
                                                                                                                                         tcttCGCGGGCTTGCCCGCGGAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgc              >  2:1866722/3‑150 (MQ=255)
                                                                                                                                         tcttCGCGGGCTTGCCCGCGGAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgc              >  2:547477/3‑150 (MQ=255)
                                                                                                                                          aTTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgcc             <  1:3414401/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgccg            >  2:2955880/1‑150 (MQ=255)
                                                                                                                                            tCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgccgc           <  2:3396152/150‑1 (MQ=255)
                                                                                                                                            tCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGgccgc           >  2:1464539/1‑150 (MQ=255)
                                                                                                                                              gcgGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCTCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCa         >  1:3918312/1‑150 (MQ=255)
                                                                                                                                                 ggCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGc      <  2:2172877/150‑1 (MQ=255)
                                                                                                                                                  gCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAg                                                             <  1:3900128/94‑1 (MQ=255)
                                                                                                                                                  gCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAg                                                             >  2:3900128/1‑94 (MQ=255)
                                                                                                                                                    aaGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                                     >  1:5354598/1‑68 (MQ=255)
                                                                                                                                                    aaGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGc                                                                                     <  2:5354598/68‑1 (MQ=255)
                                                                                                                                                     aGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTGGc  <  1:4465536/150‑1 (MQ=255)
                                                                                                                                                     aGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTGGc  <  1:4940371/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCACTTATTTGGTGCGTTTTGATTCGAAGAAGGGGATTTGTCAGTCGGGAAGGTCTGAAAGAACCAGCACAGTAGTACAGATTATGCTGCTAGTAGAGGGGTGGCGTAGTGATAGAAGGGTTTGACTGTTATGGCCCATTCGCGGGCAAGCCCGCGAAGAGGCCGGAACAGCGGAAGCAGGGCCGCCGTTCCGTACCACATCCGGATCAGCTGGCCAGCATGCGCCCGGTTTCTTCCAGGTTTTCGTGCCAGCTCAGGGCTTCGCGCAGGATGTGCGGGGTATGGCCGCCACGCTGGC  >  NC_002947/162133‑162431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: