Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,315,249 C→G 7.4% intergenic (‑598/+116) PP_16SG ← / ← msrQ 16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,315,2490CG 7.4% 496.7 / 22.1 148intergenic (‑598/+116)PP_16SG/msrQ16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (63/74);  new base G (8/3);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCGCAACA  >  NC_002947/5315100‑5315398
                                                                                                                                                     |                                                                                                                                                     
gTCGCGGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGccgcccgc                                                                                                                                                       <  2:4253854/150‑1 (MQ=255)
   gCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGACACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCcgcgcg                                                                                                                                                    >  1:2679386/1‑150 (MQ=255)
    cTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                   <  1:2888606/150‑1 (MQ=255)
    cTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                                   <  1:476908/150‑1 (MQ=255)
    cTGGAAGAGGGGCGCATTATAAGGGGATTCGAGACGGCGTCAACCTTTAATTACAAGAAAGTCAAATATAAGGGGAAATACCGTAACGGGCCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCCCgcgg                                                                                                                                                   >  2:365273/1‑150 (MQ=14)
     tGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                  >  2:310166/1‑150 (MQ=255)
       gAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                                <  2:243315/150‑1 (MQ=255)
         aGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                              <  1:304552/150‑1 (MQ=255)
          gTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                             <  1:1453677/150‑1 (MQ=255)
           tGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            >  2:1341930/1‑150 (MQ=255)
           tGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            >  2:786014/1‑150 (MQ=255)
             gggCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                          <  1:162789/150‑1 (MQ=255)
               gcgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                                                        <  2:821008/150‑1 (MQ=255)
               gcgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtctc                                                                                                                                        >  2:4481493/1‑150 (MQ=255)
                cgcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                       <  2:1546421/150‑1 (MQ=255)
                     tataAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGAc                                                                                                                                  <  2:4157584/150‑1 (MQ=255)
                     tataAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGAc                                                                                                                                  <  2:1888160/150‑1 (MQ=255)
                       taAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                                >  2:979185/1‑150 (MQ=255)
                       taAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                                <  1:3153645/150‑1 (MQ=255)
                         aGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                              <  2:3332723/150‑1 (MQ=255)
                          ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCag                                                                                                                             >  2:4056450/1‑150 (MQ=255)
                          ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCag                                                                                                                             >  2:1864687/1‑150 (MQ=255)
                           gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCaga                                                                                                                            >  1:2598136/1‑150 (MQ=255)
                           gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCaga                                                                                                                            <  2:2832919/150‑1 (MQ=255)
                           cggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCaga                                                                                                                            <  2:4099483/149‑1 (MQ=255)
                             gATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAgaga                                                                                                                          <  2:3227598/150‑1 (MQ=255)
                               ttCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGaaa                                                                                                                        <  2:2762703/150‑1 (MQ=255)
                                   agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGCGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                                    >  1:3016246/1‑150 (MQ=255)
                                   agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                                    <  1:1938111/150‑1 (MQ=255)
                                   agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                                    <  1:81246/150‑1 (MQ=255)
                                   agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                                    <  2:3573354/150‑1 (MQ=255)
                                   agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                                    <  2:32680/150‑1 (MQ=255)
                                         cGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCggg                                                                                                              <  1:1744146/150‑1 (MQ=255)
                                          gTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGc                                                                                                             <  2:4004552/150‑1 (MQ=255)
                                          gTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGc                                                                                                             >  1:2361553/1‑150 (MQ=255)
                                                       ttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGcc                                                                                                >  1:4038969/1‑150 (MQ=255)
                                                         cAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCAt                                                                                              >  1:2831811/1‑150 (MQ=255)
                                                           aGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTg                                                                                            <  1:2727714/150‑1 (MQ=255)
                                                             aaaGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCt                                                                                          <  1:2004315/150‑1 (MQ=255)
                                                                gTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtc                                                                                       <  2:74201/150‑1 (MQ=255)
                                                                 tGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                                      <  1:3105021/150‑1 (MQ=255)
                                                                    aaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                                   <  2:1071081/150‑1 (MQ=255)
                                                                    aaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                                   <  2:1071078/150‑1 (MQ=255)
                                                                    aaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                                   <  1:2320758/150‑1 (MQ=255)
                                                                       ataAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGAc                                                                                >  1:3128477/1‑150 (MQ=255)
                                                                         aaGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAg                                                                              <  1:1025020/150‑1 (MQ=255)
                                                                          aGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAgg                                                                             <  1:2799615/150‑1 (MQ=255)
                                                                          aGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAgg                                                                             >  2:1513739/1‑150 (MQ=255)
                                                                           ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGc                                                                            >  1:1683416/1‑150 (MQ=255)
                                                                           ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGc                                                                            >  2:1473737/1‑150 (MQ=255)
                                                                               aaaTACGGTTTGGGGGGTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGTACGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGTCGCCAGAGAATCAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCtt                                                                        <  1:365273/150‑1 (MQ=37)
                                                                                  tACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGc                                                                     >  1:1524099/1‑150 (MQ=255)
                                                                                   aCCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCa                                                                    <  2:3534324/150‑1 (MQ=255)
                                                                                   aCCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCa                                                                    <  1:2826668/150‑1 (MQ=255)
                                                                                     cGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                                  >  2:1784087/1‑150 (MQ=255)
                                                                                      gTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGcc                                                                 >  1:1754862/1‑150 (MQ=255)
                                                                                       gaTCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccc                                                                <  2:1270823/149‑1 (MQ=255)
                                                                                            gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccaccc                                                           <  1:3975380/150‑1 (MQ=255)
                                                                                            gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccaccc                                                           >  1:1590870/1‑150 (MQ=255)
                                                                                            gggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccaccc                                                           >  1:3237853/1‑150 (MQ=255)
                                                                                                 gCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCActc                                                      <  1:802746/150‑1 (MQ=255)
                                                                                                   cAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACtctc                                                    >  2:486881/1‑150 (MQ=255)
                                                                                                    aaGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCa                                                   <  2:4091606/150‑1 (MQ=255)
                                                                                                      gtgtTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa                                                 <  2:4060380/150‑1 (MQ=255)
                                                                                                         ttCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCg                                              <  2:629354/150‑1 (MQ=255)
                                                                                                          tCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGc                                             <  2:157810/150‑1 (MQ=255)
                                                                                                            gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGccc                                           <  1:3687944/150‑1 (MQ=255)
                                                                                                              cacaACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTa                                         >  1:4507564/1‑150 (MQ=255)
                                                                                                              cacaACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACGCGCAAACCGCCCTa                                         >  1:2256503/1‑150 (MQ=255)
                                                                                                                 aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTAccc                                      >  2:2687108/1‑150 (MQ=255)
                                                                                                                 aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTAccc                                      >  2:2813206/1‑150 (MQ=255)
                                                                                                                 aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTAccc                                      >  1:891869/1‑150 (MQ=255)
                                                                                                                  acaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCt                                     >  1:4325981/1‑150 (MQ=255)
                                                                                                                   caGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTg                                    <  1:3234556/150‑1 (MQ=255)
                                                                                                                    ttctctggcgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGc                                   >  2:2442246/8‑150 (MQ=255)
                                                                                                                    aGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGc                                   >  2:2963382/1‑150 (MQ=255)
                                                                                                                        tCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAcc                               >  2:4284976/1‑150 (MQ=255)
                                                                                                                         tggcgtCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCATAGGCGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCt                              >  2:924065/3‑150 (MQ=37)
                                                                                                                         cAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCt                              <  1:2687593/150‑1 (MQ=255)
                                                                                                                          aGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTg                             <  2:4421884/150‑1 (MQ=255)
                                                                                                                           gcgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAACCAAAGGCGGGCCCTCTCCACAGCCCTTGCTATCTCACGACAGGCCCTTGGCCGCCCCCCCCCACTCTCAAACCGCCCGACCCTGCTACCTgg                            >  1:1390184/1‑150 (MQ=38)
                                                                                                                            cgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGACACAAAGGCGGGCCCTCTCCACAGCCATTGCCATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                           >  2:3021169/1‑150 (MQ=255)
                                                                                                                                   gagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTg                    <  1:3498042/150‑1 (MQ=255)
                                                                                                                                    agctCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg                   <  2:2882905/147‑1 (MQ=255)
                                                                                                                                    agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTgg                   <  1:1095663/150‑1 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa                                                 <  2:1359573/118‑1 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa                                                 >  1:1359573/1‑118 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGaa                <  2:1759235/150‑1 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGaa                <  1:1174909/150‑1 (MQ=255)
                                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCCCCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAg               >  1:2911222/1‑150 (MQ=255)
                                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCACGTGGCAGCCCCACCCCACTCTCAAACCGCCCTAACCTGCTACCTGGCAACCTTTGGAAg               >  1:674176/1‑150 (MQ=39)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgc              <  2:2598136/150‑1 (MQ=255)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgc              >  1:697349/1‑150 (MQ=255)
                                                                                                                                          aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgcg             >  1:416015/1‑150 (MQ=255)
                                                                                                                                            cgccgcCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGcc           >  1:3152747/1‑150 (MQ=255)
                                                                                                                                             gccgcgCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCACGACAGGCCCTTGGCAGCCCCACCCCACGCGCACACCGCCCTACCCGGCTACCTGGCCACCGGGGGAAGCGccc          >  2:2156182/1‑150 (MQ=255)
                                                                                                                                                 cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGACCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCgcgc      <  1:786014/150‑1 (MQ=255)
                                                                                                                                                 cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGACCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCgcgc      <  1:1300703/150‑1 (MQ=255)
                                                                                                                                                   cgcgcgGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCGCaa    <  1:4087999/150‑1 (MQ=255)
                                                                                                                                                     cgcgGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCGCAAca  <  1:870343/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCGCAACA  >  NC_002947/5315100‑5315398

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