Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,273,948 T→A 8.5% intergenic (+19/+17) PP_4643 → / ← radA xanthine/uracil permease family protein/DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,273,9480TA 8.5% 916.3 / 58.5 303intergenic (+19/+17)PP_4643/radAxanthine/uracil permease family protein/DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base T (142/132);  new base A (10/16);  total (155/148)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

Only 100 of 309 total aligned reads displayed.

CCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCTTC  >  NC_002947/5273799‑5274097
                                                                                                                                                     |                                                                                                                                                     
 cgcgCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCtt                                                                                                                                                      <  1:3886610/150‑1 (MQ=255)
   cgcTGCGTGATGGGGACTACTTCCAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgc                                                                                                                                                    <  2:6760096/150‑1 (MQ=255)
   cgcTGCGTGATGGGGACTACTTCCAGGAGGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGCCGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgc                                                                                                                                                    <  1:6493614/150‑1 (MQ=255)
     cTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctc                                                                                                                                                  <  1:3660340/150‑1 (MQ=255)
         gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                              <  1:2376357/150‑1 (MQ=255)
         gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                              <  2:2985403/150‑1 (MQ=255)
         gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                              <  2:1562237/150‑1 (MQ=255)
         gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                              <  2:6307174/150‑1 (MQ=255)
         gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                              <  2:1639714/150‑1 (MQ=255)
          tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTcc                                                                                                                                             <  1:6138166/150‑1 (MQ=255)
          tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTcc                                                                                                                                             <  2:508662/150‑1 (MQ=255)
          tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTcc                                                                                                                                             <  1:4544345/150‑1 (MQ=255)
           gATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCt                                                                                                                                            <  1:1829659/150‑1 (MQ=255)
           gATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCt                                                                                                                                            <  2:221257/150‑1 (MQ=255)
            aTGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTCCTg                                                                                                                                           <  2:5920337/150‑1 (MQ=255)
               gggACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCAc                                                                                                                                        >  2:5802322/1‑150 (MQ=255)
               gggACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCAc                                                                                                                                        <  1:2546049/150‑1 (MQ=255)
                gggctactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAAGGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACg                                                                                                                                       >  1:6389456/4‑150 (MQ=255)
                 gACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGc                                                                                                                                      >  2:6847584/1‑150 (MQ=255)
                   ctactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTa                                                                                                                                    <  2:3958605/150‑1 (MQ=255)
                    tactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTAt                                                                                                                                   >  2:3199271/1‑150 (MQ=255)
                       ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCg                                                                                                                                >  2:6526223/1‑150 (MQ=255)
                       ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCg                                                                                                                                >  2:3891616/1‑150 (MQ=255)
                       ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCg                                                                                                                                >  2:3802143/1‑150 (MQ=255)
                       ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCg                                                                                                                                >  1:1913721/1‑150 (MQ=255)
                       ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGATGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCGAGCTCGCTCCTGCACCCTATTCg                                                                                                                                >  1:6633043/1‑150 (MQ=255)
                         cAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACCGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGaa                                                                                                                              <  1:6744275/150‑1 (MQ=255)
                                ggCAAGCTGGTCGATGCCGAGGGTAACGAAGTGCCGATGCTGTAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTGGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCa                                                                                                                       <  1:1449925/150‑1 (MQ=255)
                                  cAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATc                                                                                                                     <  2:4871719/150‑1 (MQ=255)
                                  cAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATc                                                                                                                     >  2:7184978/1‑150 (MQ=255)
                                   aaGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATcc                                                                                                                    >  1:4384415/1‑150 (MQ=255)
                                    aGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCa                                                                                                                   <  2:2684471/150‑1 (MQ=255)
                                        gTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgcgc                                                                                                               >  1:3563128/1‑150 (MQ=255)
                                            gATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGt                                                                                                           >  1:843357/1‑150 (MQ=255)
                                             aTGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGtt                                                                                                          <  1:3455858/150‑1 (MQ=255)
                                                ccGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAATCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCa                                                                                                       <  2:6270925/150‑1 (MQ=255)
                                                ccGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCa                                                                                                       <  2:3256063/150‑1 (MQ=255)
                                                        aaCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTa                                                                                               >  2:3053282/1‑150 (MQ=255)
                                                        aaCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTa                                                                                               >  1:7051480/1‑150 (MQ=255)
                                                         aCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTaa                                                                                              >  2:5372028/1‑150 (MQ=255)
                                                           gacgTGCCGTTGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACg                                                                                            <  1:3108424/147‑1 (MQ=255)
                                                           gAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACg                                                                                            >  2:2016073/1‑150 (MQ=255)
                                                               tGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCaa                                                                                        <  1:2147081/150‑1 (MQ=255)
                                                                 ccGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATc                                                                                      >  1:5361779/1‑150 (MQ=255)
                                                                   gATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCAc                                                                                    >  1:5419114/1‑150 (MQ=255)
                                                                   gATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCAc                                                                                    >  2:4081831/1‑150 (MQ=255)
                                                                    aTGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCAcc                                                                                   <  1:4002150/150‑1 (MQ=255)
                                                                    aTGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCAcc                                                                                   <  2:7051480/150‑1 (MQ=255)
                                                                     tGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCt                                                                                  <  1:711516/150‑1 (MQ=255)
                                                                     tGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCt                                                                                  <  1:2169310/150‑1 (MQ=255)
                                                                         ggAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCaa                                                                              >  1:2923311/1‑150 (MQ=255)
                                                                          gAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAc                                                                             >  1:1091429/1‑150 (MQ=255)
                                                                           aaCTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAcc                                                                            >  2:3645478/1‑150 (MQ=255)
                                                                               ggACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCGTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGc                                                                        >  2:4333415/1‑150 (MQ=255)
                                                                                gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCa                                                                       <  2:550181/150‑1 (MQ=255)
                                                                                gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCa                                                                       >  2:5944326/1‑150 (MQ=255)
                                                                                gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCa                                                                       >  2:1853337/1‑150 (MQ=255)
                                                                                gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCa                                                                       >  1:1381658/1‑150 (MQ=255)
                                                                                 aCGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAg                                                                      <  2:864449/150‑1 (MQ=255)
                                                                                 aCGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAg                                                                      >  2:6685619/1‑150 (MQ=255)
                                                                                     gCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGAGTAACGGCAATCACCTGCAACCCCGCAGgcgc                                                                  <  1:6847826/150‑1 (MQ=255)
                                                                                       aCGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGctt                                                                >  1:6786940/1‑150 (MQ=255)
                                                                                         ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttct                                                              <  1:3351854/150‑1 (MQ=255)
                                                                                            aaTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                           >  2:4805927/1‑150 (MQ=255)
                                                                                             aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCg                                                          >  2:2905214/1‑150 (MQ=255)
                                                                                             aTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCg                                                          <  2:5705238/150‑1 (MQ=255)
                                                                                               gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                        >  1:6314240/1‑150 (MQ=255)
                                                                                               gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                        <  2:1916880/150‑1 (MQ=255)
                                                                                               gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                        <  2:6959462/150‑1 (MQ=255)
                                                                                                 aaCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgcgc                                                      <  1:4223550/150‑1 (MQ=255)
                                                                                                          gcgcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCtt                                             >  1:783348/1‑150 (MQ=255)
                                                                                                          gcgcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCtt                                             <  1:4867633/150‑1 (MQ=255)
                                                                                                           cgcACGGCCGTTGCCGAGTACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                            <  2:7260097/150‑1 (MQ=255)
                                                                                                           cgcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                            <  1:1023273/150‑1 (MQ=255)
                                                                                                            gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                           >  1:974660/1‑150 (MQ=255)
                                                                                                                  ccGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGa                                     <  1:5959549/150‑1 (MQ=255)
                                                                                                                    gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                                   <  1:6935611/150‑1 (MQ=255)
                                                                                                                        ccGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGccc                               >  1:2065631/1‑148 (MQ=255)
                                                                                                                        ccGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCAc                               >  2:1304509/1‑150 (MQ=255)
                                                                                                                         cGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGaa                                                                                                                              <  2:5997658/54‑1 (MQ=255)
                                                                                                                         cGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGaa                                                                                                                              >  1:5997658/1‑54 (MQ=255)
                                                                                                                          naGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGc                             >  1:2243913/2‑150 (MQ=255)
                                                                                                                                tGACGGGTTGCAGGATCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCAGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAg                       >  2:3285869/1‑150 (MQ=255)
                                                                                                                                tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAg                       >  2:7024704/1‑150 (MQ=255)
                                                                                                                                  aCGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGcc                     >  2:4563980/1‑150 (MQ=255)
                                                                                                                                    gggTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCGGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGt                   <  2:2283913/150‑1 (MQ=255)
                                                                                                                                      gTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAgtagtcc                 <  2:6389456/150‑7 (MQ=255)
                                                                                                                                       gtGCAGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCt                >  1:5628801/2‑150 (MQ=255)
                                                                                                                                          cAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTgg             >  2:7374924/1‑150 (MQ=255)
                                                                                                                                           aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                                   <  2:4642898/127‑1 (MQ=255)
                                                                                                                                           aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                                   >  1:4642898/1‑127 (MQ=255)
                                                                                                                                             gagcgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCg          >  2:33852/1‑150 (MQ=255)
                                                                                                                                             gagcgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCg          >  1:2774866/1‑150 (MQ=255)
                                                                                                                                             gagcgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCg          >  1:7448553/1‑150 (MQ=255)
                                                                                                                                                cgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGctt       >  1:5036383/1‑150 (MQ=255)
                                                                                                                                                 gagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttc      >  2:2434355/1‑150 (MQ=255)
                                                                                                                                                 gagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttc      <  1:1362510/150‑1 (MQ=255)
                                                                                                                                                 gagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttc      <  1:7256194/150‑1 (MQ=255)
                                                                                                                                                 gagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCGTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttc      >  2:5605679/1‑150 (MQ=255)
                                                                                                                                                    cTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttcctt   <  1:3199271/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCTTC  >  NC_002947/5273799‑5274097

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: