Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,494,900 G→A 7.1% H362H (CAC→CAT cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,494,9000GA 7.1% 988.3 / 42.6 322H362H (CAC→CATcobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (152/147);  new base A (15/8);  total (167/155)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

TCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCAGGCGACC  >  NC_002947/5494751‑5495049
                                                                                                                                                     |                                                                                                                                                     
tCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcag                                                                                                                                                       >  2:3045614/1‑150 (MQ=255)
   cTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgc                                                                                                                                                    <  1:6656520/150‑1 (MQ=255)
    tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                                                                                                   >  2:7371513/1‑150 (MQ=255)
      gATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagc                                                                                                                                                 <  1:8307257/150‑1 (MQ=255)
       aTTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGGCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagca                                                                                                                                                >  1:5767943/1‑150 (MQ=255)
        ttCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgcacgcc                                                                                                                                               >  1:3693880/1‑146 (MQ=255)
          cGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCt                                                                                                                                             <  1:3216564/150‑1 (MQ=255)
           gTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTg                                                                                                                                            <  1:1144919/150‑1 (MQ=255)
                 tGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgcc                                                                                                                                      <  2:2521949/150‑1 (MQ=255)
                  gCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg                                                                                                                                     >  1:7240234/1‑150 (MQ=255)
                   cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                                                                                                    <  2:3332982/150‑1 (MQ=255)
                   cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                                                                                                    >  1:7474094/1‑150 (MQ=255)
                    ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                                                                                                   >  2:3118411/1‑150 (MQ=255)
                       gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGaa                                                                                                                                >  1:9066111/1‑150 (MQ=255)
                          cACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACCCCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACgg                                                                                                                             <  2:4296589/150‑1 (MQ=255)
                              gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGTAGCACCTGGTCGCCGCGGAACGGgcca                                                                                                                         <  2:8835117/150‑1 (MQ=255)
                              gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                                                                                                         <  2:489906/150‑1 (MQ=255)
                                   ggCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCAt                                                                                                                    <  2:1634746/150‑1 (MQ=255)
                                     ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                                                                                                  <  2:5577289/150‑1 (MQ=255)
                                     ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                                                                                                  >  2:2718367/1‑150 (MQ=255)
                                     ccATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgc                                                                                                                  >  1:7129146/1‑150 (MQ=255)
                                       aTGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATgcgc                                                                                                                <  2:646547/150‑1 (MQ=255)
                                        tGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCa                                                                                                               <  2:4519620/150‑1 (MQ=255)
                                        tGGCAAATACGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCa                                                                                                               >  1:7954723/1‑150 (MQ=255)
                                        tGGCAAATACGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTACATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCa                                                                                                               >  1:7069412/1‑150 (MQ=255)
                                            gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                                                                                                           >  2:5003288/1‑150 (MQ=255)
                                            gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTa                                                                                                           >  1:5270167/1‑150 (MQ=255)
                                             aataCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATag                                                                                                          <  2:8358807/147‑2 (MQ=255)
                                               taCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGt                                                                                                        <  1:4427159/149‑1 (MQ=255)
                                                  gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                                                                                                     >  2:7856871/1‑150 (MQ=255)
                                                  gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                                                                                                     <  2:2778246/150‑1 (MQ=255)
                                                  gCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                                                                                                     >  2:3241412/1‑150 (MQ=255)
                                                    ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGt                                                                                                   >  2:6177682/1‑150 (MQ=255)
                                                      gggTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTaa                                                                                                 >  2:7270463/1‑150 (MQ=255)
                                                       ggTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAAc                                                                                                <  2:7699576/150‑1 (MQ=255)
                                                          cGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGa                                                                                             <  2:282845/150‑1 (MQ=255)
                                                          cGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGa                                                                                             <  2:7647114/150‑1 (MQ=255)
                                                           gCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAg                                                                                            >  1:8994380/1‑150 (MQ=255)
                                                            ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc                                                                                           >  1:1971685/1‑150 (MQ=255)
                                                             cTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTTTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGca                                                                                          >  2:4205000/1‑150 (MQ=255)
                                                              tGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGcac                                                                                         >  1:1538430/1‑150 (MQ=255)
                                                                     accaccGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCg                                                                                  <  2:2968135/150‑1 (MQ=255)
                                                                      ccaccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCgg                                                                                 <  2:2193355/150‑1 (MQ=255)
                                                                       caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCTGc                                                                                <  2:6003278/150‑1 (MQ=255)
                                                                       caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                                                                                >  1:2026112/1‑150 (MQ=255)
                                                                       caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                                                                                >  2:5276049/1‑150 (MQ=255)
                                                                       caccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGc                                                                                <  1:3054988/150‑1 (MQ=255)
                                                                         ccgccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCt                                                                              <  2:6227848/150‑1 (MQ=255)
                                                                           gccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTgc                                                                            <  2:4919263/150‑1 (MQ=255)
                                                                            ccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTgcg                                                                           >  2:4792945/1‑150 (MQ=255)
                                                                               ccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCaa                                                                        <  1:8619116/150‑1 (MQ=255)
                                                                                  aaCGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGtt                                                                     >  2:6291521/1‑150 (MQ=255)
                                                                                    cGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCg                                                                   <  1:3603101/150‑1 (MQ=255)
                                                                                    cGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCg                                                                   <  1:6006229/150‑1 (MQ=255)
                                                                                    cGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCg                                                                   <  1:6447818/150‑1 (MQ=255)
                                                                                     ggCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGAAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCgg                                                                  <  1:5005503/150‑1 (MQ=255)
                                                                                        ctgaagccgcgAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATgcg                                                                                                                 <  1:2718936/95‑1 (MQ=255)
                                                                                        cTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCt                                                               >  1:6299897/1‑150 (MQ=255)
                                                                                         ttGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCTTCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCtt                                                              <  1:7531317/150‑1 (MQ=255)
                                                                                         ttGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCtt                                                              <  2:4720740/150‑1 (MQ=255)
                                                                                           gggCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGa                                                            >  2:5148419/1‑150 (MQ=255)
                                                                                              ccgccAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGAccg                                                         >  2:8931183/1‑150 (MQ=255)
                                                                                               cgccAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGAccgc                                                        <  2:7549628/150‑1 (MQ=255)
                                                                                                 ccAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGAccgccg                                                      >  1:2209401/1‑150 (MQ=255)
                                                                                                       cGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGATCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCa                                                <  1:878872/150‑1 (MQ=255)
                                                                                                        gAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCAt                                               >  1:542408/1‑150 (MQ=255)
                                                                                                        gAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCAt                                               <  2:7129146/150‑1 (MQ=255)
                                                                                                        gAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCAt                                               <  1:621708/150‑1 (MQ=255)
                                                                                                          aGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCa                                             <  1:2773478/150‑1 (MQ=255)
                                                                                                          aGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCa                                             <  2:2680503/150‑1 (MQ=255)
                                                                                                          aGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCa                                             <  2:4478685/150‑1 (MQ=255)
                                                                                                              aTGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTc                                         >  2:813566/1‑150 (MQ=255)
                                                                                                                 tCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTccgc                                      >  2:112853/1‑150 (MQ=255)
                                                                                                                 tCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTccgc                                      >  1:4868373/1‑150 (MQ=255)
                                                                                                                   gggCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGGCCGCCGGCACCATCAGTTccgccg                                    >  2:1824868/1‑150 (MQ=255)
                                                                                                                     gcgcgcTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCa                                  <  1:7641905/150‑1 (MQ=255)
                                                                                                                         gcTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAg                              >  2:3913404/1‑150 (MQ=255)
                                                                                                                         gcTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAg                              >  1:4256826/1‑150 (MQ=255)
                                                                                                                           tGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGACATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAggg                            >  1:8552369/1‑150 (MQ=255)
                                                                                                                            gCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGc                           >  2:843196/1‑150 (MQ=255)
                                                                                                                            gCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGc                           <  1:3945255/150‑1 (MQ=255)
                                                                                                                               tctcTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCaa                        <  1:3341201/150‑1 (MQ=255)
                                                                                                                               tctcTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCaa                        <  1:4962257/150‑1 (MQ=255)
                                                                                                                                 tctTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAgg                      >  1:5513064/1‑150 (MQ=255)
                                                                                                                                 tctTCGCGGTTGTCGGTGCAGAGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAgg                      <  1:728668/150‑1 (MQ=255)
                                                                                                                                   ttCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGcc                    <  1:4360919/150‑1 (MQ=255)
                                                                                                                                      gcgGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGAt                 >  2:3517960/1‑150 (MQ=255)
                                                                                                                                      gcgGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGAt                 >  2:8261509/1‑150 (MQ=255)
                                                                                                                                      gcgGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGAt                 >  1:6371100/1‑150 (MQ=255)
                                                                                                                                      gcgGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGAt                 <  2:6045414/150‑1 (MQ=255)
                                                                                                                                      gcgGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGAt                 <  2:5767943/150‑1 (MQ=255)
                                                                                                                                        ggTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGa               <  1:3118032/150‑1 (MQ=255)
                                                                                                                                          ttGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGAcc             <  2:1986240/150‑1 (MQ=255)
                                                                                                                                          ttGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGAcc             >  1:1490408/1‑150 (MQ=255)
                                                                                                                                            gTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGc           >  2:1463704/1‑150 (MQ=255)
                                                                                                                                              cGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCa         <  2:681808/150‑1 (MQ=255)
                                                                                                                                              cGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCa         <  1:2350336/150‑1 (MQ=255)
                                                                                                                                               ggTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCAg        <  2:2286445/150‑1 (MQ=255)
                                                                                                                                                    agtgcagCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCACCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCAGGCgcc   >  2:6965989/1‑148 (MQ=255)
                                                                                                                                                     gtgcagCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCAGGCGAcc  <  1:4128882/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAAGTTCGGCCTTGACCGCCGGCACCATCAGTTCCGCCGCACCAGGGCCAAGGCCGATGACCGCCAGGCGACC  >  NC_002947/5494751‑5495049

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