Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,440,461:1 +A 5.3% intergenic (+26/‑50) pcaF‑II → / → PP_2138 beta‑ketoadipyl‑CoA thiolase subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,440,4611.A 5.3% 316.0 / 8.7 76intergenic (+26/‑50)pcaF‑II/PP_2138beta‑ketoadipyl‑CoA thiolase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (33/39);  new base A (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TCCGGGGCGCGGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTG  >  NC_002947/2440313‑2440606
                                                                                                                                                     |                                                                                                                                                  
tCCGGGGCGCGGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGc                                                                                                                                                     >  6:252974/1‑149 (MQ=255)
    gggCGCGGAGTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGAGCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAg                                                                                                                                                <  6:300445/149‑1 (MQ=255)
     ggCGCGGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGt                                                                                                                                               >  8:40655/1‑149 (MQ=255)
     ggCGCGGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGt                                                                                                                                               <  2:204072/149‑1 (MQ=255)
       cgcgGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGc                                                                                                                                             >  4:101885/1‑149 (MQ=255)
           gATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCg                                                                                                                                         <  2:127443/149‑1 (MQ=255)
           gATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCg                                                                                                                                         <  6:363491/149‑1 (MQ=255)
            aTTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCACTAGGCCCCgg                                                                                                                                        <  1:205432/149‑1 (MQ=255)
            aTTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCgg                                                                                                                                        >  8:302960/1‑149 (MQ=255)
                  ggCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGgttttgt                                                                                                                                  <  5:185535/149‑1 (MQ=255)
                     aCCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGgttttgtttt                                                                                                                               <  5:17368/149‑1 (MQ=255)
                     aCCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGgttttgtttt                                                                                                                               >  7:110375/1‑149 (MQ=255)
                        cTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTc                                                                                                                            <  5:237889/149‑1 (MQ=255)
                          gCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAACCGGGGCACTAGGCCCCggctgc                                                                                                                                    >  8:102889/1‑135 (MQ=255)
                          gCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAACCGGGGCACTAGGCCCCggctgc                                                                                                                                    <  7:102889/139‑5 (MQ=255)
                              aaCGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGt                                                                                                                        <  7:302960/147‑1 (MQ=255)
                                  tCATGAAGCAAAATCGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGttttttt                                                                                                                  <  4:205553/149‑1 (MQ=255)
                                    aTGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGttttttt                                                                                                                  >  6:190821/1‑147 (MQ=255)
                                    aTGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTc                                                                                                                 <  7:123591/148‑1 (MQ=255)
                                    aTGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCa                                                                                                                <  6:406816/149‑1 (MQ=255)
                                    aTGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCAGGTTTTGTTTTTTCAGGttttttt                                                                                                                  >  6:265288/1‑147 (MQ=255)
                                      gAAGCAAAATGGTGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGgttttgtttt                                                                                                                               >  5:202250/1‑132 (MQ=255)
                                      gAAGCAAAATGGTGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGgttttgtttt                                                                                                                               <  6:202250/132‑1 (MQ=255)
                                      gAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCag                                                                                                               >  1:311998/1‑148 (MQ=255)
                                       aaGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCagga                                                                                                             <  2:311998/149‑1 (MQ=255)
                                       aaGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCag                                                                                                               >  3:209630/1‑147 (MQ=255)
                                        aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggag                                                                                                            <  4:209630/149‑1 (MQ=255)
                                        aGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCagga                                                                                                             <  7:40655/148‑1 (MQ=255)
                                         gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggagg                                                                                                           >  5:186803/1‑149 (MQ=255)
                                         gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggag                                                                                                            >  8:348185/1‑148 (MQ=255)
                                         gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggag                                                                                                            <  7:348185/148‑1 (MQ=255)
                                         gCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCCGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCaggagg                                                                                                           >  5:24830/1‑149 (MQ=255)
                                            aaaTGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGgtgt                                                                                                        >  1:247245/1‑149 (MQ=255)
                                               tGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCa                                                                                                     <  6:186803/149‑1 (MQ=255)
                                               tGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCCGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCa                                                                                                     <  6:24830/149‑1 (MQ=255)
                                                    gTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATg                                                                                                <  5:221948/149‑1 (MQ=255)
                                                        gCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCagg                                                                                                              >  1:103104/1‑131 (MQ=255)
                                                        gCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCagg                                                                                                              <  2:103104/131‑1 (MQ=255)
                                                            ggTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACGGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCA‑‑ATGCAGATACAAc                                                                                     <  4:122109/149‑1 (MQ=255)
                                                             gtgtTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACaa                                                                                      >  5:226332/1‑149 (MQ=255)
                                                                ttGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAAc                                                                                     >  3:2404/1‑147 (MQ=255)
                                                                    gACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATAAAACCAGGca                                                                               >  2:284732/1‑148 (MQ=255)
                                                                     aCCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCg                                                                               >  1:270386/1‑148 (MQ=255)
                                                                         tgtgCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATAc                                                                          >  7:26109/1‑149 (MQ=255)
                                                                              gTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCa                                                                     <  3:101885/149‑1 (MQ=255)
                                                                               tCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCag                                                                                                               <  7:251738/107‑1 (MQ=255)
                                                                               tCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCag                                                                                                               >  8:251738/1‑107 (MQ=255)
                                                                               tCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTCAGGAGgt                                                                                                          >  1:53315/1‑111 (MQ=255)
                                                                               tCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTCAGGAGgt                                                                                                          <  2:53315/111‑1 (MQ=255)
                                                                                 ggCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCAtt                                                                   <  5:252974/148‑1 (MQ=255)
                                                                                   ccTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTAc                                                                 >  2:91395/1‑148 (MQ=255)
                                                                                   ccTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTAc                                                                 >  7:161935/1‑148 (MQ=255)
                                                                                       ggCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAgg                                                            >  6:181155/1‑149 (MQ=255)
                                                                                             ggTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGa                                                      <  8:161935/149‑1 (MQ=255)
                                                                                              gTATCACCACTCTCTTCCAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCAAAGAGCCCCGGTTATGTTTTAACACGATTTTTTCAGGACGAGTGCAACATGCAGAAACAACCAGCGCTATTCCGGCAATAAAATGGGCCGGAg                                                     >  4:209493/1‑149 (MQ=255)
                                                                                                   accacTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCg                                                <  2:247245/149‑1 (MQ=255)
                                                                                                     cacTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCgtg                                              <  4:135127/149‑1 (MQ=255)
                                                                                                            ttCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACaaa                                                              >  3:219515/1‑126 (MQ=255)
                                                                                                            ttCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACaaa                                                              <  4:219515/126‑1 (MQ=255)
                                                                                                            ttCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAg                                       >  4:359275/1‑149 (MQ=255)
                                                                                                              cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTAcc                                                 <  6:24120/137‑1 (MQ=255)
                                                                                                              cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTAcc                                                 >  5:24120/1‑137 (MQ=255)
                                                                                                              cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCa                                        <  8:110375/146‑1 (MQ=255)
                                                                                                              cGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAgt                                      <  4:83877/148‑1 (MQ=255)
                                                                                                                     gTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGcgcg                              <  2:270386/149‑1 (MQ=255)
                                                                                                                     atttGATCGCGTAGCGGGACATCATCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTAc                                                                 <  8:199365/111‑1 (MQ=255)
                                                                                                                         gATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCAc                          >  1:202863/1‑149 (MQ=255)
                                                                                                                         gATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCAc                          >  5:385799/1‑149 (MQ=255)
                                                                                                                               gTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGa                     >  7:52973/1‑148 (MQ=255)
                                                                                                                                  gCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAAcc                 <  8:26109/149‑1 (MQ=255)
                                                                                                                                     ggACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGaa              <  8:52973/149‑1 (MQ=255)
                                                                                                                                     ggACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAAAAATCAGCATAAAACCGGGCTTATACCGGAATTACAAAGGCCCGGCGTACCGTGAATAAATTGTTCCCCGGACCTCCCAAACCg                >  1:120299/1‑148 (MQ=255)
                                                                                                                                            aaccGGGGCACTAGGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGa       >  2:156720/3‑149 (MQ=255)
                                                                                                                                             gcCGGGGCCTAGTGCCCCTGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGAt      <  5:265288/149‑1 (MQ=255)
                                                                                                                                             gcCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCATTCCGAAACCGAAGAGTGGAt      <  6:226332/149‑1 (MQ=255)
                                                                                                                                              ccGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATg     <  1:91395/149‑1 (MQ=255)
                                                                                                                                                 gggCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGgtg  >  6:403329/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCCGGGGCGCGGATTTCCGGCACCCTGCTCAACGTCATGAAGCAAAATGGCGGTACGCTGGGTGTTGCGACCATGTGCGTCGGCCTGGGCCAAGGTATCACCACTGTCTTCGAACGCGTCTGATCGCGTAGCGGGACAGCAGCCGGGGCCTAGTGCCCCGGTTTTGTTTTTTCAGGTTTTTTTCAGGAGGTGTGCAACATGCAGATACAACCAGGCGTATACCGGCATTACAAAGGGCCTGAGTACCGTGTCTTCAGTGTTGCGCGGCACTCCGAAACCGAAGAGTGGATGGTG  >  NC_002947/2440313‑2440606

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