Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,494,831 C→A 7.4% S385S (TCG→TCT) ‡ cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,494,8310CA 7.4% 548.2 / 21.2 167S385S (TCG→TCT) ‡cobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (71/83);  new base A (5/8);  total (76/91)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

TGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAA  >  NC_002947/5494682‑5494979
                                                                                                                                                     |                                                                                                                                                    
tgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccacc                                                                                                                                                      >  1:3399111/1‑150 (MQ=255)
 gggCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACAccaccgccaccg                                                                                                                                                     >  1:2030795/1‑150 (MQ=255)
   gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGaa                                                                                                                                                   <  2:3727493/150‑1 (MQ=255)
    cggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAAc                                                                                                                                                  >  1:2145854/1‑150 (MQ=255)
    cggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCAGGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAAc                                                                                                                                                  <  2:639752/150‑1 (MQ=255)
     ggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACg                                                                                                                                                 >  2:5250962/1‑150 (MQ=255)
       cggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGc                                                                                                                                               <  2:3339091/150‑1 (MQ=255)
       cggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGc                                                                                                                                               <  2:3338689/150‑1 (MQ=255)
            gCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                                          <  2:1535728/150‑1 (MQ=255)
            gCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTgg                                                                                                                                          <  2:828044/150‑1 (MQ=255)
              cAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGAGTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                                                                                        >  1:1775767/1‑150 (MQ=255)
                 cGAGGCGGATACCCCCGGCAAGATCTGCAGGTCACCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgccgc                                                                                                                                     <  1:3039126/150‑1 (MQ=255)
                   aGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCTGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCa                                                                                                                                   <  1:957436/150‑1 (MQ=255)
                     gCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGTATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGt                                                                                                                                 <  1:5154561/150‑1 (MQ=255)
                       ggATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTc                                                                                                                               >  2:218604/1‑150 (MQ=255)
                             ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                         >  1:4386758/1‑150 (MQ=255)
                             ccccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGc                                                                                                                         >  2:275326/1‑150 (MQ=255)
                              cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGaa                                                                                                                            <  1:5326261/146‑1 (MQ=255)
                              cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                        >  2:2405425/1‑150 (MQ=255)
                              cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                        >  2:1826349/1‑150 (MQ=255)
                              cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                        >  1:1704259/1‑150 (MQ=255)
                              cccGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCa                                                                                                                        >  1:2374939/1‑150 (MQ=255)
                                  gCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTc                                                                                                                    <  2:1271342/150‑1 (MQ=255)
                                  gCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTc                                                                                                                    <  2:545379/150‑1 (MQ=255)
                                          tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                            >  1:2804978/1‑150 (MQ=255)
                                            cAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTg                                                                                                          >  1:1183014/1‑150 (MQ=255)
                                               gTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAt                                                                                                       <  2:1108244/150‑1 (MQ=255)
                                                 cAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtct                                                                                                     >  2:3856321/1‑150 (MQ=255)
                                                  aaCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCAtctc                                                                                                    >  1:3350903/1‑150 (MQ=255)
                                                   aCGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATctct                                                                                                   <  2:3753171/150‑1 (MQ=255)
                                                      cgATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTcg                                                                                                <  1:3153901/150‑1 (MQ=255)
                                                       gATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTcgc                                                                                               <  1:4801927/150‑1 (MQ=255)
                                                         tGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCgg                                                                                             <  2:2145854/150‑1 (MQ=255)
                                                          gCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                                                                                            >  1:4964148/1‑150 (MQ=255)
                                                               gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                       <  2:4890140/150‑1 (MQ=255)
                                                               gTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                       <  2:577064/150‑1 (MQ=255)
                                                                tGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCgg                                                                                      <  2:5031421/150‑1 (MQ=255)
                                                                 ggggTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgt                                                                                     <  1:2301854/150‑1 (MQ=255)
                                                                    gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGTCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                                  <  2:1181975/150‑1 (MQ=255)
                                                                    gTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                                                                                  >  2:2301273/1‑150 (MQ=255)
                                                                      cGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagt                                                                                <  2:4030809/150‑1 (MQ=255)
                                                                       gCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtg                                                                               >  1:2572546/1‑150 (MQ=255)
                                                                         tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                             <  2:4158177/150‑1 (MQ=255)
                                                                         tGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcagtgca                                                                             <  2:4812118/150‑1 (MQ=255)
                                                                           gATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTgcagc                                                                           >  2:1878873/1‑150 (MQ=255)
                                                                                gTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCAccga                                                                      <  2:27398/150‑3 (MQ=255)
                                                                                    aaTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcg                                                                  >  1:3956622/1‑150 (MQ=255)
                                                                                     aTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgc                                                                 >  2:5302910/1‑150 (MQ=255)
                                                                                     aTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgc                                                                 <  1:46260/150‑1 (MQ=255)
                                                                                       gCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg                                                               <  1:5137799/150‑1 (MQ=255)
                                                                                        cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                              >  2:2668780/1‑150 (MQ=255)
                                                                                        cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccgc                                                              >  2:4702637/1‑150 (MQ=255)
                                                                                        cTTCGAGCACCGCAGAGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                                                                                       <  2:3046958/125‑1 (MQ=255)
                                                                                         ttCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCcgcg                                                             >  2:4909461/1‑150 (MQ=255)
                                                                                            gAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGt                                                                                         <  1:4725159/119‑1 (MQ=255)
                                                                                              gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACg                                                        >  1:3990236/1‑150 (MQ=255)
                                                                                              gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCACCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACg                                                        >  2:3484468/1‑150 (MQ=255)
                                                                                               cACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACgg                                                       <  1:2546677/150‑1 (MQ=255)
                                                                                                 ccGCAGCGGCCATGGCGAATACGCCCGTGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCTTGTCGGGCGCGCTGCATCTCTTCTCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgc                                                     <  2:5334609/150‑1 (MQ=255)
                                                                                                   gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                                   >  2:3181006/1‑150 (MQ=255)
                                                                                                   gCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgcca                                                   <  1:3699410/150‑1 (MQ=255)
                                                                                                    cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                                  <  1:1874248/150‑1 (MQ=255)
                                                                                                    cAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGgccag                                                  >  1:5055794/1‑150 (MQ=255)
                                                                                                             tcccaaaaaCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCCAGTACAAATGCATGGCGTGCACGCTGCATTTCCTCGCGAGTCTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCa                                         >  2:148457/6‑150 (MQ=255)
                                                                                                             tGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCa                                         >  1:3018642/1‑150 (MQ=255)
                                                                                                              ggCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCAc                                        <  1:2689907/150‑1 (MQ=255)
                                                                                                              ggCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCAc                                        <  1:5302910/150‑1 (MQ=255)
                                                                                                              ggCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCAc                                        >  2:2072823/1‑150 (MQ=255)
                                                                                                                 gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                            >  1:4442048/1‑79 (MQ=255)
                                                                                                                 gAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgcgc                                                                                                            <  2:4442048/79‑1 (MQ=255)
                                                                                                                  aaaaCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTaa                                    <  2:4395782/150‑1 (MQ=255)
                                                                                                                     aCGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGtt                                 >  1:1095312/1‑150 (MQ=255)
                                                                                                                     aCGCCCGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGtt                                 <  1:1928505/150‑1 (MQ=255)
                                                                                                                      cGCCCGGATCGCCTGACGATACCATCACCACATCCCGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGttt                                <  2:3139159/150‑1 (MQ=255)
                                                                                                                       gCCCGGGTCGCCTGCCGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                               <  1:462099/150‑1 (MQ=255)
                                                                                                                         ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                                     <  2:525683/110‑1 (MQ=255)
                                                                                                                         ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTgg                                                                     >  1:525683/1‑110 (MQ=255)
                                                                                                                         ccGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGt                             >  1:564129/1‑150 (MQ=255)
                                                                                                                         ccGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTcat                             >  2:2997130/1‑148 (MQ=255)
                                                                                                                         ccGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTcat                             >  1:4104977/1‑148 (MQ=255)
                                                                                                                          cGGATCGCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATa                            <  1:3595887/150‑1 (MQ=255)
                                                                                                                                gCCTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAg                      <  1:2997130/150‑1 (MQ=255)
                                                                                                                                gCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAg                      >  1:4144341/1‑150 (MQ=255)
                                                                                                                                 ccTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGc                     <  1:5361011/150‑1 (MQ=255)
                                                                                                                                 ccTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGc                     <  2:2374939/150‑1 (MQ=255)
                                                                                                                                  cTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGca                    <  2:83706/150‑1 (MQ=255)
                                                                                                                                  cTGACGATACCATCACCACATCACGGCCTGAAGCCGCGAGTTCAAATGCATGGCGTGCACGCTGCATTTCCTCGCGATTGTCGGTCATATGGATGATCTGGTCGTCCCGGAACGGGCCAGCCATGCGCACATAGGTTTCATAGCCAAGca                    >  2:2673314/1‑150 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    >  2:4636378/1‑150 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    >  1:4671293/1‑150 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    <  2:1704259/150‑1 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    >  2:1760236/1‑150 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    >  1:2093177/1‑150 (MQ=255)
                                                                                                                                  cTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGca                    >  1:1270046/1‑150 (MQ=255)
                                                                                                                                      cgacACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATc                >  2:2254451/1‑150 (MQ=255)
                                                                                                                                        acacCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCtt              <  2:2489477/150‑1 (MQ=255)
                                                                                                                                         caccacCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTc                               <  2:364642/132‑1 (MQ=255)
                                                                                                                                          accaccGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCg            >  2:4736684/1‑150 (MQ=255)
                                                                                                                                          accaccGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCg            >  1:2147756/1‑150 (MQ=255)
                                                                                                                                                gccaccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTgc      >  2:2524707/1‑150 (MQ=255)
                                                                                                                                                    ccgAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCaa  >  1:3489857/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCGCGGAACGGGCCAGCCATGCGCACGTAAGTTTCGTAACCGAGCACATCTTCGGCCTGCGCAA  >  NC_002947/5494682‑5494979

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: