Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 78,508 G→C 8.1% intergenic (‑36/+25) fmt ← / ← def‑I methionyl‑tRNA formyltransferase/peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294778,5080GC 8.1% 604.8 / 33.2 185intergenic (‑36/+25)fmt/def‑Imethionyl‑tRNA formyltransferase/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base G (74/96);  new base C (9/6);  total (83/102)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 188 total aligned reads displayed.

GCCGGCCGGGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGTACGCACACCGCCA  >  NC_002947/78359‑78657
                                                                                                                                                     |                                                                                                                                                     
gccggccgGGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCgg                                                                                                                                                       >  1:4359600/1‑150 (MQ=255)
    gccggGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCa                                                                                                                                                   >  2:4034381/1‑150 (MQ=255)
     ccggGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCTCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCaa                                                                                                                                                  <  2:3328983/150‑1 (MQ=255)
     ccggGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCaa                                                                                                                                                  <  1:5112915/150‑1 (MQ=255)
     ccggGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCaa                                                                                                                                                  <  2:3073011/150‑1 (MQ=255)
          cGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCTCCGGCAAGCCtt                                                                                                                                             <  1:3636967/150‑1 (MQ=255)
           ggTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCTCCGGCAAGCCttt                                                                                                                                            <  1:2251396/150‑2 (MQ=255)
           ggTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCTCCGGCAAGCCttt                                                                                                                                            <  1:1211212/150‑2 (MQ=255)
            gTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCTCCGGCAAGCCtttt                                                                                                                                           <  1:1754240/150‑3 (MQ=255)
            gTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCa                                                                                                                                           <  1:3838079/150‑1 (MQ=255)
            gTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCa                                                                                                                                           <  1:3782315/150‑1 (MQ=255)
            gTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAACAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCa                                                                                                                                           <  1:2162297/150‑1 (MQ=255)
                    gggTGTAGACGGCCACGATCTCGTAGGTGCTGTCGAGCAGGGACTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTACTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTg                                                                                                                                   <  1:109855/150‑1 (MQ=255)
                    gggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTg                                                                                                                                   <  2:4691061/150‑1 (MQ=255)
                    gggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTg                                                                                                                                   <  1:2849917/150‑1 (MQ=255)
                    gggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTg                                                                                                                                   <  1:2872091/150‑1 (MQ=255)
                     ggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTgg                                                                                                                                  >  1:2067878/1‑150 (MQ=255)
                     ggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTgg                                                                                                                                  >  1:3049579/1‑150 (MQ=255)
                     ggTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTgg                                                                                                                                  <  2:1857166/150‑1 (MQ=255)
                       tgtAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTgggg                                                                                                                                <  2:1096855/150‑1 (MQ=255)
                          aGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATc                                                                                                                             <  2:2499072/150‑1 (MQ=255)
                           gACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCa                                                                                                                            >  2:3406016/1‑150 (MQ=255)
                            aCGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAg                                                                                                                           <  1:1935316/150‑1 (MQ=255)
                             cGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAgg                                                                                                                          <  1:103368/150‑1 (MQ=255)
                              ggCCACGATCTCGGAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGc                                                                                                                         <  1:2874925/150‑1 (MQ=255)
                                ccACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCtt                                                                                                                       >  1:4858508/1‑150 (MQ=255)
                                 cACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTg                                                                                                                      >  1:711635/1‑150 (MQ=255)
                                   cGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCt                                                                                                                    <  2:4467090/150‑1 (MQ=255)
                                    gATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTg                                                                                                                   >  1:4681205/1‑150 (MQ=255)
                                    gATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTg                                                                                                                   <  1:3536259/150‑1 (MQ=255)
                                     aTCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTgg                                                                                                                  <  2:2095897/150‑1 (MQ=255)
                                      tctcGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGc                                                                                                                 <  2:1817532/150‑1 (MQ=255)
                                          gTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGAtg                                                                                                             >  1:834458/1‑150 (MQ=255)
                                              gggCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGAtgctgc                                                                                                         <  1:2863728/150‑1 (MQ=255)
                                              gggCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGAtgctgc                                                                                                         <  2:3100872/150‑1 (MQ=255)
                                                  tGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCtttt                                                                                                     <  1:639057/150‑1 (MQ=255)
                                                  tGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCtttt                                                                                                     >  2:4197507/1‑150 (MQ=255)
                                                  tGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCtttt                                                                                                     >  2:4197352/1‑150 (MQ=255)
                                                     cGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCa                                                                                                  <  1:4376858/150‑1 (MQ=255)
                                                      gAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAg                                                                                                 >  1:2316181/1‑150 (MQ=255)
                                                        gCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGct                                                                                               <  1:1634660/150‑1 (MQ=255)
                                                             gCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGTTCAGGCTTGCTGGCGATGCTGCTTTTCCAGcttcttc                                                                                          <  2:4782696/150‑1 (MQ=255)
                                                              ccTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGcttcttct                                                                                         <  2:541463/150‑1 (MQ=255)
                                                                    ggTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATcc                                                                                   >  1:1505817/1‑150 (MQ=255)
                                                                      tgtTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCgg                                                                                 <  2:3437032/150‑1 (MQ=255)
                                                                        ttCGGCGGCAAACTCTGGAGTGCCTGCATCGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGTTGCTGCTTTTCCAGCTTCTTCTTGATCCGGTc                                                                               <  1:2906761/150‑1 (MQ=255)
                                                                        ttCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTc                                                                               <  2:498200/150‑1 (MQ=255)
                                                                          cggcggCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTcgc                                                                             <  2:1401842/150‑1 (MQ=255)
                                                                            gcggcAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCgcgc                                                                           <  2:4330727/150‑1 (MQ=255)
                                                                               gcAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTg                                                                        >  1:1057519/1‑150 (MQ=255)
                                                                                cAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGa                                                                       >  2:2240448/1‑150 (MQ=255)
                                                                                cAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGa                                                                       <  2:3847784/150‑1 (MQ=255)
                                                                                 aaaCTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAg                                                                      >  2:2155560/1‑150 (MQ=255)
                                                                                 aaaCTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAg                                                                      >  2:4698311/1‑150 (MQ=255)
                                                                                 aaaCTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAg                                                                      >  2:2155461/1‑150 (MQ=255)
                                                                                     tctGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGc                                                                  >  1:1358662/1‑150 (MQ=255)
                                                                                     tctGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGc                                                                  <  2:4681167/150‑1 (MQ=255)
                                                                                     tctGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGc                                                                  <  2:1496677/150‑1 (MQ=255)
                                                                                     tctGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGc                                                                  >  1:3938129/1‑150 (MQ=255)
                                                                                       tGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGa                                                                <  1:1151434/150‑1 (MQ=255)
                                                                                        ggAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGAc                                                               <  1:3110449/150‑1 (MQ=255)
                                                                                          aGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAg                                                             <  2:834327/150‑1 (MQ=255)
                                                                                          aGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAg                                                             >  2:3403586/1‑150 (MQ=255)
                                                                                           gTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAgg                                                            <  1:2538909/150‑1 (MQ=255)
                                                                                               cTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAGGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGt                                                        >  2:570251/1‑150 (MQ=255)
                                                                                                tGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTc                                                       <  1:3424703/150‑1 (MQ=255)
                                                                                                  cAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCTGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTcga                                                     >  2:2644515/1‑150 (MQ=255)
                                                                                                  cAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTcga                                                     <  2:923630/150‑1 (MQ=255)
                                                                                                   aaaGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTcgac                                                    >  1:792667/1‑150 (MQ=255)
                                                                                                          aTGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGc                                             <  1:1876740/150‑1 (MQ=255)
                                                                                                            gcgcATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCtt                                           <  1:4123196/150‑1 (MQ=255)
                                                                                                              gcATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTc                                                                                                    <  2:279921/91‑1 (MQ=255)
                                                                                                              gcATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTc                                                                                                    >  1:279921/1‑91 (MQ=255)
                                                                                                              gcATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGc                                         >  2:3584369/1‑150 (MQ=255)
                                                                                                               cATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGcc                                        >  1:2261849/1‑150 (MQ=255)
                                                                                                                aTGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCg                                       <  1:1877868/150‑1 (MQ=255)
                                                                                                                  ggAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGACGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGtt                                     <  2:1472735/150‑1 (MQ=255)
                                                                                                                  ggAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGtt                                     >  2:3174709/1‑150 (MQ=255)
                                                                                                                     ggtctcGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAg                                  <  1:2474184/148‑1 (MQ=255)
                                                                                                                     gTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAg                                  <  2:1922890/150‑1 (MQ=255)
                                                                                                                     gTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAg                                  <  2:1285485/150‑1 (MQ=255)
                                                                                                                      ttCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAgg                                 >  2:1531750/1‑150 (MQ=255)
                                                                                                                        ctcGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTg                               >  2:2478444/1‑150 (MQ=255)
                                                                                                                         tcGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGa                              <  1:463971/150‑1 (MQ=255)
                                                                                                                            cTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCg                           >  1:4216925/1‑150 (MQ=255)
                                                                                                                             ttCAAAAAAGAAAATGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGa                          >  2:5048934/1‑150 (MQ=255)
                                                                                                                             ttCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGc                                                                                                                     >  1:986991/1‑59 (MQ=255)
                                                                                                                                aaaaaaGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGa                                                                                      <  2:5011219/87‑1 (MQ=255)
                                                                                                                                aaaaaaGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGa                                                                                      >  1:5011219/1‑87 (MQ=255)
                                                                                                                                  aaaagaaaaagGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCg                     >  1:1737908/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaagGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCCTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGt                    >  1:2646819/1‑150 (MQ=255)
                                                                                                                                     agaaaaagGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGc                  <  1:3227718/150‑1 (MQ=255)
                                                                                                                                     agaaaaagGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCGTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGc                  <  2:1774837/150‑1 (MQ=255)
                                                                                                                                        tgaagGCTTGCTCCGGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTgg                                                                                                                  >  2:2901606/3‑51 (MQ=18)
                                                                                                                                        tgaagGCTTGCTCCGGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGt               >  2:4036937/3‑150 (MQ=255)
                                                                                                                                        tgaagGCTTGCTCCGGCAAGCCCTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGt               >  1:3620819/3‑150 (MQ=255)
                                                                                                                                        aaaagGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGt               <  1:1410116/150‑1 (MQ=255)
                                                                                                                                               ttGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGTACGcaca        >  2:1962461/1‑150 (MQ=255)
                                                                                                                                                 gCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGTACGCACAcc      <  1:2738759/150‑1 (MQ=255)
                                                                                                                                                     gAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGTACGCACACCGCca  <  1:2641701/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGCCGGGCGGTCGGGCTGGGTGTAGACGGCCACGATCTCGTAGGGGCTGTCGAGCAGGGCCTTGAGGTGTTCGGCGGCAAACTCTGGAGTGCCTGCAAAGACGATGCGCATGGAGTTCTCGCTTCAAAAAAGAAAAAGGCTTGCCGGAGCAAGCCTTCAGGAAGGTGGGGATCAGGCTTGCTGGCGATGCTGCTTTTCCAGCTTCTTCTTGATCCGGTCGCGCTTGAGCTGCGACAGGTAGTCGACGAACAGCTTGCCGTTGAGGTGATCGAATTCGTGCTGTACGCACACCGCCA  >  NC_002947/78359‑78657

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: