Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,871,418 T→A 7.8% L217Q (CTG→CAG)  cobC → adenosylcobinamide‑phosphate synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,871,4180TA 7.8% 264.6 / 12.2 90L217Q (CTG→CAG) cobCadenosylcobinamide‑phosphate synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base A (3/4);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.65e-01

GGACAACACCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGC  >  NC_002947/1871269‑1871562
                                                                                                                                                     |                                                                                                                                                
ggACAACACCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGTGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAgct                                                                                                                                                  >  2:1560085/1‑150 (MQ=255)
ggACAACACCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAgct                                                                                                                                                  >  2:2205248/1‑150 (MQ=255)
  acaacaCCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAgctgc                                                                                                                                                >  1:1480078/1‑150 (MQ=255)
      cacCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCgg                                                                                                                                            >  2:93577/1‑150 (MQ=255)
       acCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGc                                                                                                                                           >  2:1849329/1‑150 (MQ=255)
         cccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGccc                                                                                                                                         >  2:1336416/1‑150 (MQ=255)
         cccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGccc                                                                                                                                         <  1:2192314/150‑1 (MQ=255)
         cccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGcca                                                                                                                                         <  2:239320/150‑2 (MQ=255)
          ccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGcccc                                                                                                                                        >  2:2359585/1‑150 (MQ=255)
          ccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGcccc                                                                                                                                        >  1:855703/1‑150 (MQ=255)
          ccGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGccag                                                                                                                                        <  2:1031841/150‑3 (MQ=255)
            ggcggACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTg                                                                                                                                      >  2:876949/1‑150 (MQ=255)
                       tggtggACTGGGGCGAACGCGCGGGGCTGATCGTGCTGCTCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGACGAACACAGCTTGCTGCTGCGCCTGACCGTGCTGCTCGGCCCCTTGACCGTCAACg                                                                                                                           <  2:252474/150‑1 (MQ=255)
                          tggACTGCGCCGAACGCCCGGGGCTGCTCTTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGcc                                                                                                                        <  1:1578039/150‑1 (MQ=255)
                              cTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGTGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGAc                                                                                                                    <  2:1745208/150‑1 (MQ=255)
                              cTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGAc                                                                                                                    <  1:2205248/150‑1 (MQ=255)
                              cTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGAc                                                                                                                    <  1:1009342/150‑1 (MQ=255)
                               tGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCCAcc                                                                                                                   >  1:646887/1‑150 (MQ=255)
                                  gcCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCgg                                                                                                                <  1:93577/150‑1 (MQ=255)
                                    cGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGt                                                                                                              <  1:2359585/150‑1 (MQ=255)
                                       aCGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCt                                                                                                           <  1:2171122/150‑1 (MQ=255)
                                        cGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTg                                                                                                          <  2:1430299/150‑1 (MQ=255)
                                         gCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTgg                                                                                                         <  1:2278272/150‑1 (MQ=255)
                                          cccGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGc                                                                                                        >  1:2267654/1‑150 (MQ=255)
                                            cGGGGCTGATCGTGCTGCGCTCGTTTGGNAAGTTCTTNGGCTTGGCCNGCGTGCGCCTGGGCTTCGTTGCCGCCAGACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTgaggc                                                                                                      >  1:772861/1‑146 (MQ=255)
                                            cGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGccc                                                                                                      <  1:1849329/150‑1 (MQ=255)
                                                gCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGc                                                                                                                         <  2:329842/127‑1 (MQ=255)
                                                gCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGc                                                                                                                         >  1:329842/1‑127 (MQ=255)
                                                gCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGc                                                                                                  <  1:2152454/150‑1 (MQ=255)
                                                    aTCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGc                                                                                              <  2:646887/150‑1 (MQ=255)
                                                    aTCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGATTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGGTGGCCCAGGCCAGc                                                                                              >  1:315940/1‑150 (MQ=255)
                                                    aTCGTGCTGAGATCGTTGTGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGc                                                                                              <  2:1267339/150‑1 (MQ=255)
                                                      cGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCTGGGCCAGcac                                                                                            <  2:1974367/150‑3 (MQ=255)
                                                       gTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCTGGGCCAGcgcc                                                                                           <  2:2129160/150‑4 (MQ=255)
                                                        tgctgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGt                                                                                          <  2:103720/150‑1 (MQ=255)
                                                         gctgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCt                                                                                                                                       <  2:680498/104‑1 (MQ=255)
                                                         gctgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCt                                                                                                                                       >  1:680498/1‑104 (MQ=255)
                                                          ctgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCt                                                                                        <  2:1212900/150‑1 (MQ=255)
                                                           tgcGCTCGTTTGGCAAGTTCTTCGGCTTGGGCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGccag                                                                                                                                        >  1:1596603/1‑99 (MQ=255)
                                                           tgcGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGccag                                                                                                                                        <  2:1596603/101‑3 (MQ=255)
                                                               cTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTTCTGGCCCAGGCCAGCCTGTCtgatga                                                                                   >  2:883220/1‑150 (MQ=255)
                                                               cTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCtgatga                                                                                   <  2:402538/150‑1 (MQ=255)
                                                                tCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATgag                                                                                  >  2:1840354/1‑150 (MQ=255)
                                                                      ggCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCa                                                                                                <  2:1537638/130‑1 (MQ=255)
                                                                      ggCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCa                                                                                                >  1:1537638/1‑130 (MQ=255)
                                                                      ggCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGcc                                                                            <  2:772861/150‑1 (MQ=255)
                                                                         aaGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAg                                                                         >  2:2260385/1‑150 (MQ=255)
                                                                          aGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgc                                                                        <  2:855703/150‑1 (MQ=255)
                                                                            ttcttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCa                                                                                                <  1:1842083/124‑1 (MQ=255)
                                                                            ttcttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCa                                                                                                >  2:1842083/1‑124 (MQ=255)
                                                                               ttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgcgcgcg                                                                   >  1:608193/1‑150 (MQ=255)
                                                                               ttcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGAACCGTCAACGGCCCGACCCGGGTGCTAGCCCAGGCGAGCCTGTCTGATGAGACTGCCCAgcgcgcg                                                                   >  1:213136/1‑150 (MQ=255)
                                                                                tcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGccag                                                                                                                                        <  1:2346169/80‑3 (MQ=255)
                                                                                tcGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCAGCTCGGccag                                                                                                                                        >  2:2346169/1‑78 (MQ=255)
                                                                                  ggCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTAGCCAGCGCGCGCAg                                                                >  2:1043764/1‑150 (MQ=255)
                                                                                  ggCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTCGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAg                                                                <  1:2034773/150‑1 (MQ=255)
                                                                                        cccGGCAGGCGCCTGGGCTTCGTTGCCGCCGAAGACACCCTGCTGCTCCCCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCaa                                                          <  2:315940/149‑1 (MQ=255)
                                                                                          cGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACg                                                        >  2:1145154/1‑150 (MQ=255)
                                                                                              gTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCtgt                                                    <  1:775443/150‑1 (MQ=255)
                                                                                               tGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCtgtg                                                   <  1:876949/150‑1 (MQ=255)
                                                                                               tGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCtgtg                                                   <  2:385134/150‑1 (MQ=255)
                                                                                                  gcCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGtgctg                                                >  2:616403/1‑150 (MQ=255)
                                                                                                   ccTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTTGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGt                                                                                          <  1:1530144/107‑1 (MQ=255)
                                                                                                   ccTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTTGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGt                                                                                          >  2:1530144/1‑107 (MQ=255)
                                                                                                   ccTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGtgctgc                                               <  1:2260385/150‑1 (MQ=255)
                                                                                                    cTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgc                                                                        <  1:2113743/124‑1 (MQ=255)
                                                                                                    cTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAgc                                                                        >  2:2113743/1‑124 (MQ=255)
                                                                                                        gCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTTCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTgcca                                          <  1:883220/150‑1 (MQ=255)
                                                                                                        gCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTgcca                                          <  1:1145154/150‑1 (MQ=255)
                                                                                                          ttCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTgccagc                                        >  1:643080/1‑150 (MQ=255)
                                                                                                          ttCGTTGCCGACGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTgccagc                                        >  2:1192564/1‑150 (MQ=255)
                                                                                                              ttGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGccc                                                                                                      >  2:1476074/1‑84 (MQ=255)
                                                                                                              ttGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGccc                                                                                                      <  1:1476074/84‑1 (MQ=255)
                                                                                                              ttGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCtgatga                                                                                   <  2:726196/103‑1 (MQ=255)
                                                                                                              ttGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCtgatga                                                                                   >  1:726196/1‑103 (MQ=255)
                                                                                                              ttGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCcagccagc                                    >  2:1571690/1‑150 (MQ=255)
                                                                                                                  cgccgAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAg                                                                                               <  2:17318/87‑1 (MQ=255)
                                                                                                                  cgccgAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAg                                                                                               >  1:17318/1‑87 (MQ=255)
                                                                                                                        acacAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCa                          >  2:413162/1‑150 (MQ=255)
                                                                                                                        acacAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCa                          >  1:2326652/1‑150 (MQ=255)
                                                                                                                               ttGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTg                   >  2:526345/1‑150 (MQ=255)
                                                                                                                                tgctgctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTgc                  <  2:2326652/150‑1 (MQ=255)
                                                                                                                                   tgctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCg                                   >  2:664109/1‑130 (MQ=255)
                                                                                                                                   tgctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCg                                   <  1:664109/130‑1 (MQ=255)
                                                                                                                                    gctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAg              <  1:1717378/150‑1 (MQ=255)
                                                                                                                                    gctgcGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCAAGGACAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGACAGGCGCACGTGTTGCGCAg              >  1:1328969/1‑150 (MQ=255)
                                                                                                                                         gcCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCg         >  1:1829969/1‑150 (MQ=255)
                                                                                                                                          ccTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCAGACCCGGGTGCTGGACAAGGCCAGCCTGTCTGATGAGACTGCCTAGCCTGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACCTGTTGCGCAGTGCCgg        >  2:50212/1‑150 (MQ=255)
                                                                                                                                           tgggCCGAGCAGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGc       >  1:1756055/3‑150 (MQ=255)
                                                                                                                                             ggCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCt     <  1:663120/150‑1 (MQ=255)
                                                                                                                                             ggCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCt     <  1:1298402/150‑1 (MQ=255)
                                                                                                                                              gCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTg    <  1:468340/150‑1 (MQ=255)
                                                                                                                                                cGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGgc  >  1:1437257/1‑150 (MQ=255)
                                                                                                                                                   gctgctCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCCGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCg                                   >  2:770885/1‑114 (MQ=255)
                                                                                                                                                   gctgctCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCCGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCg                                   <  1:770885/114‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGACAACACCCCGGCGGACAGCGTGGTGGACTGCGCCGAACGCCCGGGGCTGATCGTGCTGCGCTCGTTTGGCAAGTTCTTCGGCTTGGCCGGCGTGCGCCTGGGCTTCGTTGCCGCCGAACACAGCTTGCTGCTGCGCCTGGCCGAGCTGCTCGGCCCCTGGACCGTCAACGGCCCGACCCGGGTGCTGGCCCAGGCCAGCCTGTCTGATGAGACTGCCCAGCGCGCGCAGGCCCAACGCTGTGCTGCTGCCAGCCAGCGGCTTGCGCACGTGTTGCGCAGTGCCGGCCTGGC  >  NC_002947/1871269‑1871562

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: