Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,790,734 C→G 5.5% G45G (GGC→GGG PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,790,7340CG 5.5% 287.7 / 6.4 91G45G (GGC→GGGPP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (38/48);  new base G (2/3);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

TGGTCAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTCGCGC  >  NC_002947/3790587‑3790882
                                                                                                                                                   |                                                                                                                                                    
tGGTCAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGccc                                                                                                                                                    <  1:495996/150‑1 (MQ=255)
tGGTCAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGccc                                                                                                                                                    >  2:1774666/1‑150 (MQ=255)
   tCAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCTGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTg                                                                                                                                                 <  2:1100113/150‑1 (MQ=255)
    cAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTgg                                                                                                                                                <  1:1154191/150‑1 (MQ=255)
     aaGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGc                                                                                                                                               >  2:946546/1‑150 (MQ=255)
     aaGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGc                                                                                                                                               <  1:1673404/150‑1 (MQ=255)
       gCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCga                                                                                                                                             >  1:1035316/1‑150 (MQ=255)
       gCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCga                                                                                                                                             >  1:954106/1‑150 (MQ=255)
       gCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCga                                                                                                                                             <  2:1463766/150‑1 (MQ=255)
       gCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGCACTGGTTACAGCGGCCTGCCCGCGCAGCGCCAAGGGCCTTGGCgc                                                                                                                                             <  1:1024868/150‑2 (MQ=255)
        ccGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGGCCTTGGCgct                                                                                                                                            <  1:2204088/150‑3 (MQ=255)
        ccGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCgag                                                                                                                                            >  1:912057/1‑150 (MQ=255)
         cGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCgaga                                                                                                                                           >  1:1557235/1‑150 (MQ=255)
         cGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCgaga                                                                                                                                           <  1:1704247/150‑1 (MQ=255)
          gCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGacac                                                                                                                                          <  1:2400084/150‑1 (MQ=255)
          gCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGAc                                                                                                                                          <  2:1257361/150‑1 (MQ=255)
          gCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGAc                                                                                                                                          <  1:784889/150‑1 (MQ=255)
            aaTGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGCGACgg                                                                                                                                        >  1:357930/1‑150 (MQ=255)
            aaTGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGACTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACgg                                                                                                                                        >  2:1861306/1‑150 (MQ=255)
                aTCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATg                                                                                                                                    >  2:1344893/1‑150 (MQ=255)
                aTCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATg                                                                                                                                    >  1:354137/1‑150 (MQ=255)
                   gCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTg                                                                                                                                 >  2:144636/1‑150 (MQ=255)
                   gCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTg                                                                                                                                 <  1:1839709/150‑1 (MQ=255)
                    cccGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGc                                                                                                                                >  2:783819/1‑150 (MQ=255)
                    cccGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGc                                                                                                                                >  2:605404/1‑150 (MQ=255)
                     ccGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCa                                                                                                                               <  2:2060243/150‑1 (MQ=255)
                        ggCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGt                                                                                                                            <  2:904110/150‑1 (MQ=255)
                              cGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCggg                                                                                                                      <  1:617051/150‑1 (MQ=255)
                               ggCTGCGCGGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGggc                                                                                                                     <  2:6464/150‑1 (MQ=255)
                                  tGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGggcggc                                                                                                                  <  1:1492202/150‑1 (MQ=255)
                                  tGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGACTATGCCGTTTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGggcggc                                                                                                                  >  1:1665835/1‑150 (MQ=255)
                                    cgcAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGa                                                                                                                <  2:2280625/150‑1 (MQ=255)
                                    cgcAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGa                                                                                                                <  1:643077/150‑1 (MQ=255)
                                     gcagcgCTCGCCACCGGCCTGCTCGCCGCGGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCGGCCCGCGCAGCGCTAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAg                                                                                                               >  1:1311218/1‑150 (MQ=255)
                                     gcagcaCTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAg                                                                                                               >  2:201128/1‑150 (MQ=255)
                                      cagcaCTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGt                                                                                                              <  1:201128/150‑1 (MQ=255)
                                       agcaCTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGtt                                                                                                             <  1:1491872/150‑1 (MQ=255)
                                             cGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCtt                                                                                                       <  1:783819/150‑1 (MQ=255)
                                               ccACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTg                                                                                                     <  1:144636/150‑1 (MQ=255)
                                               ccACCGGCATGCTCGCCGACGCCTATGCCGTGTCGCCGCAAGCCCCGGCGCCGCCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTg                                                                                                     <  2:954106/150‑1 (MQ=255)
                                                 aCCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCga                                                                                                                                             <  1:520622/108‑1 (MQ=255)
                                                 aCCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCga                                                                                                                                             >  2:520622/1‑108 (MQ=255)
                                                 aCCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGt                                                                                                   >  1:1364617/1‑150 (MQ=255)
                                                   cGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTc                                                                                                 >  1:2021102/1‑150 (MQ=255)
                                                     gCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCAGCGCAGCGCCAAGGCCCTTGGAGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCaa                                                                                               <  1:1975517/150‑1 (MQ=255)
                                                          cTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGc                                                                                          <  1:1344893/150‑1 (MQ=255)
                                                           tCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTCCAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCt                                                                                         >  1:696961/1‑150 (MQ=255)
                                                           tCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGGCCTTGGCgctgc                                                                                                                                          <  1:39197/101‑5 (MQ=255)
                                                           tCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGGCCTTGGCgctgc                                                                                                                                          >  2:39197/1‑97 (MQ=255)
                                                             gccgccgccTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCttt                                                                                       <  2:2265925/150‑1 (MQ=255)
                                                             ccctccgccgATGCCGTGTGGCCGCAAGACCCGGAGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCttt                                                                                       <  2:1600582/146‑1 (MQ=255)
                                                               cgccgccTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTa                                                                                     >  2:1902829/1‑150 (MQ=255)
                                                                    ccTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTacgacg                                                                                >  2:1405668/1‑150 (MQ=255)
                                                                         gCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACgaagaa                                                                           <  1:1110162/150‑1 (MQ=255)
                                                                         gCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTCGACTGGCGGCATGGTGCACGTGCCGGGCGGCGAGTTCCGCTTTGGTTCAAGCCGCTGTTACGACgaagaa                                                                           >  1:151433/1‑150 (MQ=255)
                                                                          ccGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGaagaag                                                                          <  1:1902829/150‑1 (MQ=255)
                                                                           cGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCtt                                                                                                                                                  >  2:2107368/1‑77 (MQ=255)
                                                                           cGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCtt                                                                                                                                                  <  1:2107368/77‑1 (MQ=255)
                                                                               tgGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCg                                                                     >  2:1864268/1‑150 (MQ=255)
                                                                                 gCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGcc                                                                   >  1:1008270/1‑150 (MQ=255)
                                                                                   cGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTc                                                                 >  1:1712083/1‑150 (MQ=255)
                                                                                    gCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCa                                                                >  2:2117141/1‑150 (MQ=255)
                                                                                      aaGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCAcc                                                              >  2:2178578/1‑150 (MQ=255)
                                                                                             ggcgccggcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCaa                                                       <  2:1035316/150‑1 (MQ=255)
                                                                                                  cggcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGgtgt                                                  <  1:946546/150‑1 (MQ=255)
                                                                                                   ggcgcTGCTGTGTACGGGCGACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGCCGGCGGGGCGCATGTGCCGGGCGGCGAGTTCAGCCTTGGTTTACGCATCTTTGACGTCGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTc                                                 <  1:1150497/150‑1 (MQ=255)
                                                                                                    gcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCaa                                                                                               >  1:2380122/1‑103 (MQ=255)
                                                                                                    gcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCaa                                                                                               >  1:2380164/1‑103 (MQ=255)
                                                                                                    gcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCaa                                                                                               <  2:2380164/103‑1 (MQ=255)
                                                                                                    gcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCaa                                                                                               <  2:2380122/103‑1 (MQ=255)
                                                                                                    gcgcTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTcc                                                <  1:2117141/150‑1 (MQ=255)
                                                                                                      gctgctGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCgg                                              >  1:604610/1‑150 (MQ=255)
                                                                                                            gtgtACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTg                                        <  1:1405668/150‑1 (MQ=255)
                                                                                                              gtACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGCTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGa                                      <  1:169633/150‑1 (MQ=255)
                                                                                                               tACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGAt                                     >  1:669688/1‑150 (MQ=255)
                                                                                                                  gggCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGggcggc                                                                                                                  >  2:1921923/1‑70 (MQ=255)
                                                                                                                  gggCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGggcggc                                                                                                                  <  1:1921923/70‑1 (MQ=255)
                                                                                                                   ggCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGAc                                 <  1:1864268/150‑1 (MQ=255)
                                                                                                                   ggCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGAc                                 <  2:1346033/150‑1 (MQ=255)
                                                                                                                   ggCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGAc                                 <  1:605404/150‑1 (MQ=255)
                                                                                                                     cTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACgg                               <  2:1364617/150‑2 (MQ=255)
                                                                                                                      tACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTg                              <  1:1886315/150‑1 (MQ=255)
                                                                                                                      tACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTg                              <  2:2008878/150‑1 (MQ=255)
                                                                                                                         aGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCAt                           >  1:2302899/1‑150 (MQ=255)
                                                                                                                             gCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCgg                       <  2:354137/150‑1 (MQ=255)
                                                                                                                              ccTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGt                      <  2:2302899/150‑1 (MQ=255)
                                                                                                                               cTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTc                     >  1:1613178/1‑150 (MQ=255)
                                                                                                                                    cgcgCAGCGCCAAGGCCCTTGGCGAGACGGCATGATGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCaa                >  1:1769860/1‑150 (MQ=255)
                                                                                                                                          tcgcCAAGGGCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCa          >  1:1302453/2‑150 (MQ=255)
                                                                                                                                               aaGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTcg     >  2:1937174/1‑150 (MQ=255)
                                                                                                                                                  gCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTcgcgc  <  2:1309654/150‑1 (MQ=255)
                                                                                                                                                  gCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTcgcgc  >  1:283151/1‑150 (MQ=255)
                                                                                                                                                  gCCCTTGGCGAGACGGCATGGTGCACGGGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTcgcgc  <  2:912057/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGTCAAGCCGCAATGATCGCCCGGGCAGCCGGCTGCGCAGCACTCGCCACCGGCCTGCTCGCCGCCGCCTATGCCGTGTGGCCGCAAGCCCCGGCGCCGGCGCTGCTGTGTACGGGCTACAGCGGCCTGCCCGCGCAGCGCCAAGGCCCTTGGCGAGACGGCATGGTGCACGTGCCGGGCGGCGAGTTCAGCTTTGGTTCAAGCCGCTTTTACGACGAAGAAGGCCCGCCTCACCCCGCCAAGGTGTCCGGCTTCTGGATTGACGTGCATCCGGTCACCAACGCCCAGTTCGCGC  >  NC_002947/3790587‑3790882

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: