Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,882,252:1 +A 6.2% coding (1290/1356 nt) uacT → uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,882,2521.A 6.2% 383.5 / 18.7 97coding (1290/1356 nt)uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base . (46/45);  new base A (2/4);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

CCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCG  >  NC_002947/4882105‑4882401
                                                                                                                                                     |                                                                                                                                                     
ccgcAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGt                                                                                                                                                     >  2:1898570/1‑150 (MQ=255)
       cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGt                                                                                                                                                >  1:1964857/1‑148 (MQ=255)
       cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGt                                                                                                                                                <  2:1964857/148‑1 (MQ=255)
       cAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGGGGATGGAGCCCATCACCCCAGGCGCCTTCTCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTc                                                                                                                                              >  1:2378574/1‑150 (MQ=255)
          cTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACAGTACGTTCAAc                                                                                                                                           <  1:2116682/150‑1 (MQ=255)
          cTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACAGTACGTTCAAc                                                                                                                                           <  2:1778291/150‑1 (MQ=255)
           tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACa                                                                                                                                          >  1:2053318/1‑150 (MQ=255)
           tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACa                                                                                                                                          >  2:446078/1‑150 (MQ=255)
           tGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTAGACGAAGCCCATCCCCCACAGTTGAATGCTGTGTTCCACGGTTTTTGCCCTGGTGTTGAACAGTACGTTCAACa                                                                                                                                          <  1:414869/150‑1 (MQ=255)
                gttgCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCtt                                                                                                                                     <  2:530517/150‑1 (MQ=255)
                   gCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGc                                                                                                                                  <  1:446078/150‑1 (MQ=255)
                            gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCgg                                                                                                                                >  1:1317361/1‑143 (MQ=255)
                            gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCgg                                                                                                                                <  2:1317361/143‑1 (MQ=255)
                            gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGc                                                                                                                             >  2:1347143/1‑146 (MQ=255)
                            gTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGc                                                                                                                             <  1:1347143/146‑1 (MQ=255)
                             tCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGAcc                                                                                                                        >  1:706876/1‑150 (MQ=255)
                                 cTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCgcgc                                                                                                                    >  2:2452281/1‑150 (MQ=255)
                                  ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACAGTACGTTCAACa                                                                                                                                          >  2:20318/1‑127 (MQ=255)
                                  ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACAGTACGTTCAACa                                                                                                                                          <  1:20318/127‑1 (MQ=255)
                                  ttCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCa                                                                                                                   >  1:754463/1‑150 (MQ=255)
                                   tCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCaa                                                                                                                  <  1:2232292/150‑1 (MQ=255)
                                    cGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTAAACATCCTTGGCGGTGCCGACCGCGCAAt                                                                                                                 <  1:918145/150‑1 (MQ=255)
                                     gggTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATg                                                                                                                >  1:266602/1‑150 (MQ=255)
                                       gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAAc                                                                                                                                                        <  1:692758/109‑1 (MQ=255)
                                       gTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGca                                                                                                              <  1:1050712/150‑1 (MQ=255)
                                       gTTGATCCCCGTGGTGCGTCAGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAAc                                                                                                                                                        >  2:692758/1‑109 (MQ=255)
                                          gATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  1:1650689/150‑1 (MQ=255)
                                          gATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  2:1745698/150‑1 (MQ=255)
                                          gATCCCCGTGGTGCGTCCGGAGTTCCTTGCACAGACGCCACGGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACaa                                                                                                           <  2:123423/150‑1 (MQ=255)
                                           aTCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAac                                                                                                          <  2:2378859/150‑1 (MQ=255)
                                             ccccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                          <  1:496659/132‑1 (MQ=255)
                                             ccccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGa                                                                                                                          >  2:496659/1‑132 (MQ=255)
                                               ccGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgacg                                                                                                      >  2:1478319/1‑150 (MQ=255)
                                               ccGTGGTACGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAacgacg                                                                                                      >  1:1311987/1‑150 (MQ=255)
                                                cGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGt                                                                                                     >  1:2298572/1‑150 (MQ=255)
                                                  tgctgCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCt                                                                                                   <  1:1416191/147‑1 (MQ=255)
                                                  gggtgCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCt                                                                                                   <  2:1475925/149‑1 (MQ=255)
                                                         tCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCa                                                                                            >  2:2383008/1‑150 (MQ=255)
                                                            ggAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCa                                                                                         >  1:335694/1‑150 (MQ=255)
                                                            ggAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGAATCGCCATGGCCACGGTCAAACCCCGATTGTTGAACGTACTGTTCAACATCCTTGGCGGGGCCGACCGCGCAATGCACAACGACGTGTGTCACCAGcg                                                                                         >  1:535098/1‑149 (MQ=255)
                                                             gAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCAt                                                                                        <  2:266602/150‑1 (MQ=255)
                                                              aGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCCGCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCAtt                                                                                       <  2:535098/150‑1 (MQ=255)
                                                                 tCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAg                                                                                    <  1:1898570/150‑1 (MQ=255)
                                                                   tttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                                  >  2:760296/1‑150 (MQ=255)
                                                                   tttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGcc                                                                                  >  2:137794/1‑150 (MQ=255)
                                                                    ttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCg                                                                                                                                 >  1:2156858/1‑102 (MQ=255)
                                                                    ttGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCg                                                                                                                                 <  2:2156858/102‑1 (MQ=255)
                                                                       cccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGTGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgg                                                                              <  2:1167724/150‑1 (MQ=255)
                                                                       cccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgg                                                                              <  2:1284825/150‑1 (MQ=255)
                                                                        ccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAggg                                                                             >  1:1587276/1‑150 (MQ=255)
                                                                        ccAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAggg                                                                             >  2:252260/1‑150 (MQ=255)
                                                                         cAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggg                                                                            >  1:1703980/1‑150 (MQ=255)
                                                                           gaggCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggggg                                                                          <  2:285808/147‑1 (MQ=255)
                                                                           gATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAgggggg                                                                          <  1:1478319/150‑1 (MQ=255)
                                                                               ccACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGt                                                                      <  1:1558884/150‑1 (MQ=255)
                                                                                cacaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTTAACATCCTTGGCGGTGCCGACCGCGCCATGCACAACGACGTCTGTCACGAGCATTGAGCCAAGGGGGGCGGTg                                                                     >  1:1343318/1‑150 (MQ=255)
                                                                                cacaGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTg                                                                     >  1:1491991/1‑150 (MQ=255)
                                                                                 acaGTGGATGGAGCCAATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGc                                                                    >  2:394834/1‑150 (MQ=255)
                                                                                   aGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCAt                                                                  <  1:137794/150‑1 (MQ=255)
                                                                                      ggatggaGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGc                                                               <  2:754463/150‑1 (MQ=255)
                                                                                         tggaGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCAGTAGACGGTTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGc                                                            <  1:258090/150‑1 (MQ=255)
                                                                                           gaGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGccc                                                          <  2:2298572/150‑1 (MQ=255)
                                                                                              cccATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCg                                                                        <  2:1790927/133‑1 (MQ=255)
                                                                                              cccATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCg                                                                        >  1:1790927/1‑133 (MQ=255)
                                                                                              cccATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTc                                                       <  2:2031023/150‑1 (MQ=255)
                                                                                              cccATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTc                                                       <  2:706876/150‑1 (MQ=255)
                                                                                                cATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTACATGGCTGCCCCTCGt                                                     >  1:1304709/1‑150 (MQ=255)
                                                                                                 aTCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGtt                                                    <  2:978515/150‑1 (MQ=255)
                                                                                                    aCCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGAGCATGGCTGCCCCTCGtgccc                                                 >  1:97532/1‑146 (MQ=255)
                                                                                                        acaGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCaa                                                                                                                                            <  2:1800768/55‑1 (MQ=255)
                                                                                                        acaGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCaa                                                                                                                                            >  1:1800768/1‑55 (MQ=255)
                                                                                                           gCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCGCCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          <  2:743090/150‑1 (MQ=255)
                                                                                                               tATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTaa                                      >  1:2217405/1‑150 (MQ=255)
                                                                                                                aTCGCCATGGCCATGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                     <  2:710441/150‑1 (MQ=255)
                                                                                                                aTCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAAc                                     <  2:1587276/150‑1 (MQ=255)
                                                                                                                 tCGCCATGGCCACGGTCTGAGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                      <  1:101226/132‑1 (MQ=255)
                                                                                                                 tCGCCATGGCCACGGTCTGAGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCg                                                      >  2:101226/1‑132 (MQ=255)
                                                                                                                   gCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                <  2:233589/120‑1 (MQ=255)
                                                                                                                   gCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                >  1:233589/1‑120 (MQ=255)
                                                                                                                     cATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     <  2:2371670/97‑1 (MQ=255)
                                                                                                                     cATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     >  1:2371670/1‑97 (MQ=255)
                                                                                                                            cACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGt                          >  2:880052/1‑150 (MQ=255)
                                                                                                                                  cTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     <  2:863807/85‑1 (MQ=255)
                                                                                                                                  cTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGa                                                                                     >  1:863807/1‑85 (MQ=255)
                                                                                                                                   tgtgCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccg                   >  1:1199182/1‑150 (MQ=255)
                                                                                                                                     tgCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACccgcc                 >  1:874866/1‑150 (MQ=255)
                                                                                                                                       cccTGGTGTTGAACGTACTGTGCAACAGCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCgg               >  2:1403240/1‑150 (MQ=255)
                                                                                                                                       caaggatgtTGAACAGTACGTTCAACACCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGCGCGGTGCATGGCTGCCCCTCGCTCACCTCTGCAGTAACCGCTACCCCGTCCACCCGCCgg               >  2:518758/7‑150 (MQ=255)
                                                                                                                                        ccTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGa              >  1:1939550/1‑150 (MQ=255)
                                                                                                                                        ccTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGa              >  2:2021026/1‑150 (MQ=255)
                                                                                                                                         cTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGaa             <  2:2053318/150‑1 (MQ=255)
                                                                                                                                          tGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          >  2:539255/1‑120 (MQ=255)
                                                                                                                                          tGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCa                                          <  1:539255/120‑1 (MQ=255)
                                                                                                                                          tGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAAt            <  2:874866/150‑1 (MQ=255)
                                                                                                                                          tGGNGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAAt            >  2:1516874/1‑150 (MQ=255)
                                                                                                                                            gtgtTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                <  1:2450863/96‑1 (MQ=255)
                                                                                                                                            gtgtTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATgg                                                                >  2:2450863/1‑96 (MQ=255)
                                                                                                                                              gtTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGACTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGt        <  2:1311987/150‑1 (MQ=255)
                                                                                                                                                 gAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGc     >  2:2246585/1‑150 (MQ=255)
                                                                                                                                                    cGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCg  >  2:2311491/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCAACCAGCTGTTGGTTGCCGTGAGCGTCGGCTTCGGGTTGATCCCCGTGGTGCGTCCGGAGTTCTTTGCCCAGATGCCACAGTGGATGGAGCCCATCACCCACAGCGGTATCGCCATGGCCACGGTCTGTGCCCTGGTGTTGAACGTACTGTTCAACATCCTTGGCGGTGCCGACCGCGCAATGCACAACGACGTCTGTCACCAGCATTGAGCCAAGGGGGGCGGTGCATGGCTGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCG  >  NC_002947/4882105‑4882401

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