Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,575,435 T→A 7.9% intergenic (+52/‑136) pcaF‑I → / → pcaT beta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,575,4350TA 7.9% 518.3 / 12.0 194intergenic (+52/‑136)pcaF‑I/pcaTbeta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (73/100);  new base A (15/2);  total (88/106)
Fisher's exact test for biased strand distribution p-value = 2.56e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

CGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTCAACCTA  >  NC_002947/1575286‑1575584
                                                                                                                                                      |                                                                                                                                                     
   tGCACCAGCTGGAGATGAGTGGCGGTCGCTAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCAc                                                                                                                                                    <  2:7143011/150‑1 (MQ=255)
   tGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCAc                                                                                                                                                    <  2:1484735/150‑1 (MQ=255)
       ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                                <  2:8783549/150‑1 (MQ=255)
       ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                                <  2:2935762/150‑1 (MQ=255)
       ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                                <  2:743270/150‑1 (MQ=255)
       ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                                >  2:6404103/1‑150 (MQ=255)
        cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                                               <  1:821316/150‑1 (MQ=255)
        cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                                               >  2:4291093/1‑150 (MQ=255)
        cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtggg                                                                                                                                               >  2:5755737/1‑148 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGCGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                              <  2:3838738/150‑1 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGGATCCACGCAATCACCTGTgg                                                                                                                                              <  1:853506/150‑1 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                              <  2:6144906/150‑1 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                              <  1:5403422/150‑1 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                              <  2:7529735/150‑1 (MQ=255)
         aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCACCAAAAGGGATACACGCAAACAccaggtg                                                                                                                                              >  1:586545/1‑145 (MQ=255)
          gCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCAAAACGGAAACACCCAAAAACCAGTggg                                                                                                                                             >  1:6702613/1‑150 (MQ=255)
             ggAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGc                                                                                                                                          <  1:5786942/150‑1 (MQ=255)
              gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                                                         <  2:771081/150‑1 (MQ=255)
               agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                                                        <  1:2278894/150‑1 (MQ=255)
               agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                                                        <  1:3432368/150‑1 (MQ=255)
               agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                                                        <  2:3668557/150‑1 (MQ=255)
               agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                                                        <  2:5475029/150‑1 (MQ=255)
                gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCggg                                                                                                                                       <  1:5000497/150‑1 (MQ=255)
                 aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                      <  2:5438920/150‑1 (MQ=255)
                 aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                      <  1:1434686/150‑1 (MQ=255)
                 aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                      <  1:4918570/150‑1 (MQ=255)
                 aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                      <  1:4104429/150‑1 (MQ=255)
                 aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                      <  1:2729734/150‑1 (MQ=255)
                     gTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCACCAACACGGATCCACGCAAACACCAGTGGGAGCGGGCGTcc                                                                                                                                  >  1:5295071/1‑150 (MQ=255)
                      acatggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTAGCGGGGAAGCCCCCACCAAAAACGATACAACCAAAAAACAGTGGGAGCGGGGTTccc                                                                                                                                 >  1:3128041/5‑150 (MQ=255)
                       ggCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTcccc                                                                                                                                <  1:5938292/150‑1 (MQ=18)
                        atggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCGTCGCGGGGGCGCCCCCGCCAACACGGGGACACGCCAAAACCAGTGGGCGGCGGGTTccccc                                                                                                                               >  1:4499712/3‑149 (MQ=255)
                         tggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCAACAGGGATCCACGCAAAAACCGGGGGGAGCGGGCGTCCCcgc                                                                                                                              >  1:5589952/2‑150 (MQ=255)
                         tggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCGGTGGGAGCGGGCTTCCCcgc                                                                                                                              >  2:3328199/2‑150 (MQ=255)
                          ggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAAAAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTACCcgcg                                                                                                                             >  2:2321067/1‑150 (MQ=255)
                          ggTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCTCCAACACGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCcgcg                                                                                                                             >  2:1721982/1‑150 (MQ=255)
                              gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCGCCAACACGGATACACGCAAAAAACAGTGGGTGCGGGCGTCCCCCCGCaga                                                                                                                         >  1:3199226/1‑150 (MQ=255)
                              gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                                                         <  2:1981198/150‑1 (MQ=255)
                              gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCAGTGGGCGCGGGCGTCCCCCCGAaga                                                                                                                         >  2:1836527/1‑150 (MQ=255)
                               cAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAagag                                                                                                                        >  2:7239409/1‑150 (MQ=255)
                                aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGCTCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAgg                                                                                                                       >  2:8095143/1‑150 (MQ=255)
                                            aCCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTga                                                                                                           <  2:4856682/150‑1 (MQ=255)
                                              cATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACCCAATCACCGGTGGGAGCGGGCTTCCCCGCGCAGAGGCCGCTCCTGCgagc                                                                                                         >  2:7632855/1‑149 (MQ=255)
                                              cATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCTCCAACAGGGATCCACCCAATAACCAGTGGGTGCGGGCTTCCCCCCGAAGAGGCCACACCGGTgaga                                                                                                         >  1:2120381/1‑150 (MQ=255)
                                                           gtCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCACTCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCAtt                                                                                            <  2:6504105/150‑1 (MQ=255)
                                                           gtCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCAtt                                                                                            <  2:3603405/150‑1 (MQ=255)
                                                              ggCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACCTAAGGCCCAGTGcc                                                                                         >  2:5655715/1‑150 (MQ=255)
                                                               gCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCTCCAAAACGGCTCCACCCAATCACCAGTGGGAGCCGGCTTCCCCCCGCAGAGGCCTCACCGGTGAGACCTAAGGGCCAGTGCCa                                                                                        >  2:417878/1‑150 (MQ=255)
                                                                 cAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGGGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGc                                                                                      <  2:4536195/150‑1 (MQ=255)
                                                                  aaGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGCAGCCCGCTCCAACAGGGATACACGCCATCACCGGTGGGCGCGGGCGTCCCCGCGAAGAGGCCGCTCCAGTGAGCCATAAGGCCCAATGCCAGcc                                                                                     >  1:4984356/1‑150 (MQ=255)
                                                                   aGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCTGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGccc                                                                                    <  2:5870286/150‑1 (MQ=255)
                                                                     gTCTGGCGTTGGCCAGCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCCCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCTAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGa                                                                                  <  2:2290050/150‑1 (MQ=255)
                                                                        tGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAATGCCCATTGCCAGCCCGATAg                                                                               >  2:3451474/1‑150 (MQ=255)
                                                                        tGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAg                                                                               >  2:7836164/1‑150 (MQ=255)
                                                                             ttGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGCAGCCCGCTCCAACAGGGATCCCCGCCATCCCCGGTGGGCGCGGGCGTCCCCGCGCCGCGGCCTGACCGGTGAGACATAAGGCCCCGTGCCAGCCCGAGAGCGGac                                                                          >  2:4710890/1‑149 (MQ=255)
                                                                              tGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAAc                                                                         <  1:7167858/150‑1 (MQ=255)
                                                                               ggCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCCTTGCCAGCCCGATAGCGGCACg                                                                        >  2:5269619/1‑150 (MQ=255)
                                                                                gCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCAGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGa                                                                       <  2:3413529/150‑1 (MQ=255)
                                                                                 ccATCGTGGGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAg                                                                      <  2:2611322/150‑1 (MQ=255)
                                                                                   aTCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCTCCAACAGGGGTCCACGCAATCACCAGTGGGTGCGGGCGTCCCCGCGGAGAGGCCGCTCCAGTGAGCCCAAAGGCCCCTTGCCCGCCCGATAGCGGAACGAgtg                                                                    >  1:7927876/1‑150 (MQ=255)
                                                                                      gggggggTTTCACCTGTACCGGCCTCTTGGGGGGCACGCCCGCTCCAACTGGGACCCACGAAATCCCCTGTGGGAGGGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgtgt                                                                 <  1:3305477/146‑1 (MQ=255)
                                                                                            ttttcaacTGTACCGCCCTATTCGCGGTCAACCCCCAACAAACAGCCAACAACCAAAACCCCTGTGGGAGCGGCCACCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAg                                                           <  1:964025/143‑1 (MQ=255)
                                                                                             tttGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAgg                                                          <  1:604675/150‑1 (MQ=255)
                                                                                              ttGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGAAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGt                                                         <  2:3462395/150‑1 (MQ=255)
                                                                                                gACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAt                                                       <  2:1021039/150‑1 (MQ=255)
                                                                                                gACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAt                                                       <  2:5182138/150‑1 (MQ=255)
                                                                                                gACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCAGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAt                                                       >  2:735558/1‑150 (MQ=255)
                                                                                                   cTGTACCGGCCTCTTCGCGGGGAAGCCCCCTCCAACAGGGATACACGCAATAACCTGTGGGAGCGGGCTTACCCGCGAAGAGGCCACTACAGTTAGACAAAAAGCCCAATTCCAGCCCCATAGCGGAACGAGTGTGTTAAAAGGTAGGTc                                                    >  2:1188614/1‑150 (MQ=255)
                                                                                                      tACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAg                                                 <  2:6559300/150‑1 (MQ=255)
                                                                                                            ccTCTTGGGGGGGACGCCCGCGCCACCTGGGATCCACGCAATCACGGGGGGGTGGGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTAc                                           <  1:1133114/150‑1 (MQ=255)
                                                                                                               ctTCGGGGGGAACCCCGCCCAAACTGGGATCCCCGAAACCACGTGTGGGCGCGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa                                        <  1:6733544/150‑1 (MQ=255)
                                                                                                               ctTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa                                        <  2:3442299/150‑1 (MQ=255)
                                                                                                                  cgcgGGGACGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGa                                     <  2:4892183/150‑1 (MQ=255)
                                                                                                                     ggggAAGCCCGCTCCACCAGGGAGCCACGCAATCACCTGTGTGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGAccc                                  <  1:5831739/150‑1 (MQ=255)
                                                                                                                      gggAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCt                                 >  2:2336264/1‑150 (MQ=255)
                                                                                                                        gAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTcc                               >  1:7287914/1‑150 (MQ=255)
                                                                                                                          aGCCCGCTCCCACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                             <  2:6705250/150‑1 (MQ=255)
                                                                                                                          aGCCCGCTCCAACTGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                             <  1:6404103/150‑1 (MQ=255)
                                                                                                                          aGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                             <  2:4455775/150‑1 (MQ=255)
                                                                                                                            cccGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           <  2:3705154/150‑1 (MQ=255)
                                                                                                                            cccGCTCCAACAGGGATCCACGCAATCACGTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           <  2:7042728/150‑1 (MQ=255)
                                                                                                                            cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           >  1:7584607/1‑150 (MQ=255)
                                                                                                                            cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           >  2:3032838/1‑150 (MQ=255)
                                                                                                                            cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           >  2:6877865/1‑150 (MQ=255)
                                                                                                                            cccGCTCCAACAGGCATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                           <  2:3128041/150‑1 (MQ=255)
                                                                                                                             ccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAg                          >  1:5155579/1‑150 (MQ=255)
                                                                                                                                 tCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTaaa                      >  1:7981222/1‑150 (MQ=255)
                                                                                                                                  ccAACAGGGATCCACGCAATCACCTGGGGGAGCGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAc                     <  1:704484/150‑1 (MQ=255)
                                                                                                                                   cAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGt                         >  1:3825170/1‑145 (MQ=255)
                                                                                                                                    aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAccc                   >  1:1768554/1‑150 (MQ=255)
                                                                                                                                      cAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATCGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGt                 >  2:7784432/1‑150 (MQ=255)
                                                                                                                                       aGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTc                >  1:1936/1‑150 (MQ=255)
                                                                                                                                       aGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTc                >  1:1191917/1‑150 (MQ=255)
                                                                                                                                         ggATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCAt              >  2:2117085/1‑150 (MQ=255)
                                                                                                                                           gATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATg             >  1:6795227/1‑150 (MQ=255)
                                                                                                                                             tCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGAc           >  2:1225518/1‑150 (MQ=255)
                                                                                                                                               cACGCACTCACCTGGGGGAGCGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCt         <  2:586545/150‑1 (MQ=255)
                                                                                                                                                aCGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTc        >  2:4034838/1‑150 (MQ=255)
                                                                                                                                                aCGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTc        >  2:8561518/1‑150 (MQ=255)
                                                                                                                                                    aaTCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTCAAcc    >  1:2591509/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTCAACCTA  >  NC_002947/1575286‑1575584

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