Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,249,705 A→C 8.9% K428Q (AAG→CAG)  ureC → urease subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,249,7050AC 8.9% 310.8 / 17.3 109K428Q (AAG→CAG) ureCurease subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (52/47);  new base C (5/5);  total (57/52)
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.91e-01

Only 100 of 111 total aligned reads displayed.

GCACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTT  >  NC_002947/3249557‑3249853
                                                                                                                                                    |                                                                                                                                                    
gcACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCaa                                                                                                                                                     >  1:604529/1‑150 (MQ=255)
gcACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGAACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCaa                                                                                                                                                     >  1:1948331/1‑150 (MQ=255)
  aCCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGc                                                                                                                                                   >  1:1001944/1‑150 (MQ=255)
   ccAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCt                                                                                                                                                  <  2:1287237/150‑1 (MQ=255)
    cAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTg                                                                                                                                                 >  1:909479/1‑150 (MQ=255)
    cAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTg                                                                                                                                                 >  2:1368875/1‑150 (MQ=255)
      gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                               >  2:748308/1‑150 (MQ=255)
      gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                               >  1:1165214/1‑150 (MQ=255)
      gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTAGGTGGAGGTCGGCAAGCTGGc                                                                                                                                               >  1:532452/1‑150 (MQ=255)
         gAAGCCGCGGCGCGGGGCGCTGGCACCGGACACCCCGTAAAGCGACAACTTCCGGGTTAAGCGCTAAATCGCCAAGGACACCATCCACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGa                                                                                                                                            <  2:1214873/150‑1 (MQ=255)
           aGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGAcc                                                                                                                                          <  1:2318146/150‑1 (MQ=255)
           aGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAcc                                                                                                                                          <  1:1898922/150‑1 (MQ=255)
            gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAcc                                                                                                                                          <  1:2216655/149‑1 (MQ=255)
            gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAcc                                                                                                                                          >  2:2216655/1‑149 (MQ=255)
            gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCt                                                                                                                                         <  1:905716/150‑1 (MQ=255)
              tGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTgg                                                                                                                                       <  1:746676/150‑1 (MQ=255)
                    gcgGGGCGCTGGCACCGGACACCCCGTACAGCGACAACTGCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgt                                                                                                                                 <  2:172158/150‑1 (MQ=255)
                          cGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGc                                                                                                                           >  2:1910894/1‑150 (MQ=255)
                             tGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAg                                                                                                                        <  1:748308/150‑1 (MQ=255)
                             tGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCAAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGGCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAg                                                                                                                        >  2:1290517/1‑150 (MQ=255)
                               gCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCg                                                                                                                      >  2:1798477/1‑150 (MQ=255)
                                  ccGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCCTGCTTGCTGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTc                                                                                                                   >  1:762387/1‑150 (MQ=255)
                                     gACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCttt                                                                                                                <  2:1165214/150‑1 (MQ=255)
                                          cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGcc                                                                                                                          >  2:2139205/1‑135 (MQ=255)
                                          cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGcc                                                                                                                          <  1:2139205/135‑1 (MQ=255)
                                           ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                <  1:291720/128‑1 (MQ=255)
                                           ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGGCGGCAAGCTGGCCGACCTGGTGCtgtg                                                                                                                                >  2:291720/1‑128 (MQ=255)
                                           ccGCACAGCGACAACTAACGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCACCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTc                                                                                                          >  2:1602819/1‑150 (MQ=255)
                                             gTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGAAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCaa                                                                                                        <  2:2269160/150‑1 (MQ=255)
                                                    gACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGATGTGGTCGCAGGAGGTCGTTGCTTTCAAGCCGGcc                                                                                                 >  1:1020357/1‑150 (MQ=255)
                                                       aaCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCt                                                                                                                                                  >  1:1503586/1‑98 (MQ=255)
                                                       aaCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCt                                                                                                                                                  <  2:1503586/98‑1 (MQ=255)
                                                         cTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                            <  1:506761/150‑1 (MQ=255)
                                                          ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAGCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTCGCGCACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                           >  1:1954136/1‑150 (MQ=255)
                                                          ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                           >  1:1564390/1‑150 (MQ=255)
                                                             cGGGTCAAGCGCTACATCGCCAAGTACACCATCTACCCGTCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTc                                                                                        <  1:860037/150‑1 (MQ=255)
                                                                gTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                             <  2:2052645/110‑1 (MQ=255)
                                                                gTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                             >  1:2052645/1‑110 (MQ=255)
                                                                 tCATGCCCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAgg                                                                                    >  2:1298214/1‑150 (MQ=255)
                                                                     gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCtt                                                                                                                 <  2:1920486/117‑1 (MQ=255)
                                                                     gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCtt                                                                                                                 >  1:1920486/1‑117 (MQ=255)
                                                                     gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCgg                                                                                >  1:1112503/1‑150 (MQ=255)
                                                                       gcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGga                                                                              <  1:613667/150‑1 (MQ=255)
                                                                        cTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgat                                                                             >  2:2285896/1‑150 (MQ=255)
                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACTAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTACTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatgccc                                                                         <  2:1611209/150‑4 (MQ=255)
                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                            >  1:1556952/1‑131 (MQ=255)
                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGTATGATc                                                                         <  2:532452/150‑1 (MQ=255)
                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                         <  1:927805/150‑1 (MQ=255)
                                                                            aTCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAACCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                            <  2:1556952/131‑1 (MQ=255)
                                                                                ccAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCa                                                                     <  1:1910894/150‑1 (MQ=255)
                                                                                ccAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGAGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGAGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCACGGGCGGGCTGATCGtaa                                                                     >  2:134084/1‑148 (MQ=255)
                                                                                 cAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCAc                                                                    <  1:1368875/150‑1 (MQ=255)
                                                                                    gTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCaa                                                                                      >  2:82874/1‑129 (MQ=255)
                                                                                    gTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCaa                                                                                      <  1:82874/129‑1 (MQ=255)
                                                                                     tACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGc                                                                                                              <  1:2189334/104‑1 (MQ=255)
                                                                                     tACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGc                                                                                                              >  2:2189334/1‑104 (MQ=255)
                                                                                       cacCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                           >  1:1862716/1‑121 (MQ=255)
                                                                                       cacCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTg                                                                                           <  2:1862716/121‑1 (MQ=255)
                                                                                       cacCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatg                                                                            >  1:218513/1‑136 (MQ=255)
                                                                                       cacCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGgatg                                                                            <  2:218513/136‑1 (MQ=255)
                                                                                       cacCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGcc                                                              >  1:1261968/1‑150 (MQ=255)
                                                                                          cATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTTCTCAAGGGCGGGATGATCGTCACCGCGCCCAt                                                           >  2:1108807/1‑150 (MQ=255)
                                                                                           aTCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCAtc                                                          <  2:1653395/150‑2 (MQ=255)
                                                                                            tCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCAGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                         >  1:1718304/1‑150 (MQ=255)
                                                                                              caCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGc                                                       <  1:1108807/149‑1 (MQ=255)
                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCgg                                                                                >  1:2071726/1‑125 (MQ=255)
                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCgg                                                                                <  2:2071726/125‑1 (MQ=255)
                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATggga                                                       <  2:604529/150‑2 (MQ=255)
                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                         <  1:1912200/148‑1 (MQ=255)
                                                                                              aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                         >  2:1912200/1‑148 (MQ=255)
                                                                                                  cGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                            >  2:1034201/1‑109 (MQ=255)
                                                                                                  cGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                            <  1:1034201/109‑1 (MQ=255)
                                                                                                    gCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGga                                                                              >  2:1778500/1‑121 (MQ=255)
                                                                                                    gCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGga                                                                              <  1:1778500/121‑1 (MQ=255)
                                                                                                    gCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                         >  1:432884/1‑126 (MQ=255)
                                                                                                    gCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATc                                                                         <  2:432884/126‑1 (MQ=255)
                                                                                                       cTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAAc                                              <  1:1162147/150‑1 (MQ=255)
                                                                                                         cACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACgg                                            <  1:1074861/150‑1 (MQ=255)
                                                                                                            ccATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTc                                         <  2:1001944/150‑1 (MQ=255)
                                                                                                             cATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTcc                                        <  1:2285896/150‑1 (MQ=255)
                                                                                                                    tCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGa                                                     >  2:323712/1‑130 (MQ=255)
                                                                                                                    tCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATacccc                                 >  2:1089540/1‑149 (MQ=255)
                                                                                                                      ggTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATacccacc                               <  1:715296/150‑1 (MQ=255)
                                                                                                                        tCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCAcccc                             >  2:1996331/1‑150 (MQ=255)
                                                                                                                          aCGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCAcccccc                           >  2:2303718/1‑150 (MQ=255)
                                                                                                                           cgtaGTGGGCTCGGTTGAGCTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGTCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGAGATCAACGGCTCCATACCCACCCCCCa                          <  2:1020357/147‑1 (MQ=255)
                                                                                                                              aGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                         <  2:217823/116‑1 (MQ=255)
                                                                                                                              aGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATgg                                                         >  1:217823/1‑116 (MQ=255)
                                                                                                                                   gCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc                  >  2:1694528/1‑150 (MQ=255)
                                                                                                                                    cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCa                 >  1:1107244/1‑150 (MQ=255)
                                                                                                                                         tGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATc            <  2:2401839/150‑1 (MQ=255)
                                                                                                                                           caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCAc                                                                    >  2:1980227/2‑92 (MQ=255)
                                                                                                                                           caGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGACAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCAc                                                                    <  1:1980227/91‑1 (MQ=255)
                                                                                                                                            aGGTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc         >  1:776156/1‑150 (MQ=255)
                                                                                                                                             ggTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATggg                                                        >  2:2054632/1‑102 (MQ=255)
                                                                                                                                             ggTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATggg                                                        <  1:2054632/102‑1 (MQ=255)
                                                                                                                                               tCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCgg                                                                                >  1:1081606/1‑76 (MQ=255)
                                                                                                                                               tCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGa      >  2:820727/1‑150 (MQ=255)
                                                                                                                                                cGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGAt     >  1:2408406/1‑150 (MQ=255)
                                                                                                                                                   cAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGtt  >  2:651808/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACCAGATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAAGGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCCGATGTT  >  NC_002947/3249557‑3249853

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