Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,870,471 G→C 5.2% intergenic (‑97/+35) gllP ← / ← glllA porin‑like protein/gallate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,870,4710GC 5.2% 873.3 / 22.9 276intergenic (‑97/+35)gllP/glllAporin‑like protein/gallate dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (150/110);  new base C (2/13);  total (152/124)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 280 total aligned reads displayed.

TGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGACACACCC  >  NC_002947/2870323‑2870620
                                                                                                                                                      |                                                                                                                                                     
tggctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGa                                                                                                                                                      >  1:123164/1‑150 (MQ=255)
  gctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcg                                                                                                                                                    >  2:8765705/1‑150 (MQ=255)
  gctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcg                                                                                                                                                    >  1:3822858/1‑150 (MQ=255)
   ctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcgc                                                                                                                                                   >  1:8812620/1‑150 (MQ=255)
    tgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcgcg                                                                                                                                                  >  1:7916511/1‑150 (MQ=255)
    tgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcgcg                                                                                                                                                  >  1:4664826/1‑150 (MQ=255)
    tgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcgcg                                                                                                                                                  >  1:3614797/1‑150 (MQ=255)
      cctgTGCGATTGTGGGGACGCGCGGGTACCGGGTGCGGCATTTCATTATTTTTATTCGCCCGGGAGGGGCAAC‑GC‑GCGGCCGGCCAGCGGCGCATCGGCGCTGGGCCGGTTGAAGCCGGGGGCGGTGGGGGTGGCGGCGCGTCCCGCGcc                                                                                                                                                <  1:1184044/147‑1 (MQ=255)
         gTGCGAATGTGGGGACGAGCGGGTACCGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCGGCCCGCCAGCGGCGCATCGGCGCTGGGCAGTTTGAATCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCgcc                                                                                                                                             <  1:7173683/150‑3 (MQ=255)
          tGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTca                                                                                                                                            >  2:3078066/1‑150 (MQ=255)
          tGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCGCca                                                                                                                                            <  2:2417960/150‑1 (MQ=255)
           gcgcATGTGGGGCCGCGCGGGTACCGGGGGCGGCATTTCATTATTTTGATTCTCCCGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAAGCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCGCcac                                                                                                                                           <  2:6906893/150‑1 (MQ=255)
             gAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCacac                                                                                                                                         >  1:3366105/1‑150 (MQ=255)
                      gACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGGTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTGACACGATCAAGCg                                                                                                                                >  1:6028198/1‑150 (MQ=255)
                      gACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCg                                                                                                                                >  2:3365376/1‑150 (MQ=255)
                       aggcgggggTACTGGGTGGGGCTTTACTTTTTTTTTTTGCCCCCGGGTTGGCCAAC‑GC‑GCTCCCCGCCTGCGGCGAACCGGCGCTGGGCTGGTTCAATCGGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCgg                                                                                                                               <  2:1214210/144‑1 (MQ=255)
                        ggggCGGGTACGGGGGGGGGATTTACTTTATTTTTTTGCTCCGGGGATGGGCAAC‑GC‑GCGGCCCGCCGGGGGCGCAGGGGCGCGGGGCTGGTTGAATCGGGGGGGGGGGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCggg                                                                                                                              <  2:8309749/147‑1 (MQ=255)
                          cggggttACTGGGTGGGGATTTTATTTTTTTTTTTTCTCCTGGGTTGGGAAAA‑GC‑GGTGCCCGCCTGCGGGGCAGGGGCGGTGGGCTGTTTGATTCTGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGtt                                                                                                                            <  2:287106/144‑1 (MQ=255)
                              ggTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCAt                                                                                                                        >  2:8480934/1‑150 (MQ=255)
                                    gggTGGGGCATTTATTTTTTTTTATTCGCCTGGGTTGGCAAAC‑GC‑GGGGCCCGCCAGCGGGGCATCGGCGGTGGGCTGGTTGAAGCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCa                                                                                                                  <  2:6309804/150‑1 (MQ=255)
                                     ggTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAg                                                                                                                 >  2:6818176/1‑150 (MQ=255)
                                         cGGCTTTACTTTATTTTTATTCCCCCGGGGTGGGCAAC‑GC‑GCGCCCCGCCAGCGGCGCTGGGGCGCTGGGCTGTTTGATGCTGGGGGGGGTGGGGGTGGCGGCGGGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTgg                                                                                                             <  2:7620671/150‑1 (MQ=255)
                                          ggCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTggg                                                                                                            >  2:2790139/1‑150 (MQ=255)
                                          ggCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTggg                                                                                                            >  1:6713020/1‑150 (MQ=255)
                                           ggattttattttTTTTTTTTCCCCTGGTTGGGGAAACGC‑GCTGCCCGCCCGGGGCGAATGGGGGGTGGGCTGTTTGAACCCGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTggg                                                                                                            <  2:3951921/148‑1 (MQ=255)
                                             aTTACATTATTTTTTTTCTCCTGGGTTGGCCAAC‑GC‑GCGGCCCGCCCGGGGCGCTTGGGCGCGGGGCTGGTTGAATCCGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGgcgc                                                                                                         <  2:4995962/150‑1 (MQ=255)
                                              tttCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCt                                                                                                        >  2:5732815/1‑150 (MQ=255)
                                              tttCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCt                                                                                                        >  1:322356/1‑150 (MQ=255)
                                                 ccttttttttCTCCCCCGGGGAT‑GGCAAA‑CCGGCGGCCGGCCCGCGGCGCAGGGGGGCGGGGCTGGTTGAACCGGGGGGGGGGGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTg                                                                                                     <  1:5808190/145‑1 (MQ=255)
                                                  atttttTTTTTTCCCCGGGGTTGGCAAAC‑GC‑GGGGCCCCCCCGCGGGGCCGGGGCGGGGGGCTGGTTGAAGCCGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgc                                                                                                    <  2:8216226/150‑1 (MQ=255)
                                                  attattTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgc                                                                                                    >  1:2922697/1‑150 (MQ=255)
                                                   ttattTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgcc                                                                                                   >  1:8568922/1‑150 (MQ=255)
                                                    ttttttttttCCCCGGGGTGGGGCAAC‑GG‑GCGGCCCCCCCGCGGCGCATGGGGGCGGGGCTGGTTGAACCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgcca                                                                                                  <  1:5359965/148‑1 (MQ=255)
                                                    tattTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgcca                                                                                                  >  2:406096/1‑150 (MQ=255)
                                                      tttttATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgccagc                                                                                                >  1:5785610/1‑150 (MQ=255)
                                                      tttttATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTCGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCAGCGGTCAGTTGGGCGCTTTgccagc                                                                                                >  2:708443/1‑150 (MQ=255)
                                                            tcccccAGGCATGGCCACC‑GC‑GCGGCCCGCCCGCGGCGAACCGGCGCGGGGCCGGTTGAATCCGGGGGCGGGGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGa                                                                                          <  2:1614887/146‑1 (MQ=255)
                                                               tcCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCACCCACCGAGTa                                                                                       >  1:8163940/1‑150 (MQ=255)
                                                                cccgggTTGGGCAAC‑GC‑GGGGCCCGCCTGCGGCGCAGCGGCGGGGGGCCGGTTCACTCTGGGGGCGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAg                                                                                      <  2:2205460/147‑1 (MQ=255)
                                                                cccgggTTGGGCAAACGC‑GCGCCCCCCCCGGGGCGCCGCGGCGGGGGGCTGGTTCAACCGGGGGGGGGGGGGGGGGGGGGCGGGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTa                                                                                       <  1:5293561/147‑1 (MQ=255)
                                                                ccAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCACCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAg                                                                                      >  2:2744176/1‑150 (MQ=255)
                                                                 cAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGc                                                                                     >  1:1652246/1‑150 (MQ=255)
                                                                  aGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGTAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCg                                                                                    >  2:8030376/1‑150 (MQ=255)
                                                                   gggTTGGGCAAC‑GC‑GCGGCCCGCCCGGGGCGCTTGGGCGGTGGCCTGGTGGAATCCGGGGGCGTGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGa                                                                                   <  2:5528420/150‑1 (MQ=255)
                                                                    ggttGGGCCAC‑GC‑GGGGCCCGCCCGGGGCGAAGCGGCGGGGGGCCGGTTGAACCGGGGGGGGGTGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGaa                                                                                  <  2:2542794/147‑1 (MQ=255)
                                                                    ggATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGaa                                                                                  >  2:5945607/1‑150 (MQ=255)
                                                                    ggATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTCGCCAGCCACCGAGTAGCGaa                                                                                  >  1:8032232/1‑150 (MQ=255)
                                                                         ggCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgc                                                                             >  1:894087/1‑150 (MQ=255)
                                                                           cAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgctg                                                                           >  2:6224673/1‑150 (MQ=255)
                                                                            aaacgc‑gcGGCCCGCCGGGGGCGCAGGGGCGCGGGGCTGGTTGAATCGGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgctg                                                                           <  2:1068339/146‑1 (MQ=255)
                                                                                gc‑gcTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGGCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTa                                                                       >  2:5954066/1‑150 (MQ=255)
                                                                                  ggcggCCCGCCGGGGGCGCTGCGGGGGTGGGCTGTTTCATGCCGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTAc                                                                      <  1:6696966/146‑1 (MQ=255)
                                                                                   gcTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACg                                                                     >  2:833249/1‑150 (MQ=255)
                                                                                    gtGCCGGCCCGCGGGGATTCGGGGCTGGGCTGGTTGATTCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACgg                                                                    <  2:7514620/149‑1 (MQ=255)
                                                                                    cggcccgcccGCGGCGCAGGGGGGGGGGGCCGTTTTAATCCGGGGGGGGGGGGGGGGGCGGGGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACgg                                                                    <  1:7764219/148‑1 (MQ=255)
                                                                                        ccGCCCGGGGCGAATGGGGGGTGGGCTGGTTTACGCGGGGGGGGGGGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCt                                                                <  1:4620104/150‑1 (MQ=255)
                                                                                               gggcgCATGGGGGGTGGGCTGTTTTAATCGGGGGGGGGGGGGGGTGGCGGGGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGaa                                                         <  1:1321761/149‑1 (MQ=255)
                                                                                                   gcATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTcc                                                     >  1:1820802/1‑150 (MQ=255)
                                                                                                       ggggggtGGGCTGGTTTAATCGGGGGGGGTTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGa                                                 <  2:8514840/144‑1 (MQ=255)
                                                                                                       cggcggGGGGCTGGTTTATTCGGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGa                                                 <  1:838491/150‑1 (MQ=255)
                                                                                                          cgcTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGc                                              >  2:446589/1‑150 (MQ=255)
                                                                                                            cggggCTGGTTTAATCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGctc                                            <  2:5288523/148‑1 (MQ=255)
                                                                                                              gggCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTg                                          >  1:3646598/1‑150 (MQ=255)
                                                                                                               ggCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGc                                         >  1:5444241/1‑150 (MQ=255)
                                                                                                                gCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGcc                                        >  2:220137/1‑150 (MQ=255)
                                                                                                                gCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGcc                                        >  1:3498393/1‑150 (MQ=255)
                                                                                                                    gTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCcgcg                                    >  2:3991370/1‑150 (MQ=255)
                                                                                                                          tCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATggc                              >  2:7791165/1‑150 (MQ=255)
                                                                                                                          tCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATggc                              >  1:7318611/1‑150 (MQ=255)
                                                                                                                          tCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCAGCTGGTTACGGGTCTTCATGAATTCCTCGATGCTCTGCCCCCGCATggc                              >  1:248396/1‑150 (MQ=255)
                                                                                                                             gggggCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATggcggc                           <  2:7728122/150‑1 (MQ=255)
                                                                                                                             gggggCGGTGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATggcggc                           <  1:6518991/150‑1 (MQ=255)
                                                                                                                             gggggCGGGGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATggcggc                           <  1:8425632/150‑1 (MQ=255)
                                                                                                                              ggggCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGgcggcg                          >  2:3206583/1‑150 (MQ=255)
                                                                                                                              ggggCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGgcggcg                          >  2:7547191/1‑150 (MQ=255)
                                                                                                                              ggggCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATAAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCCAGTAGCGAACCAGCTGGTTACGGGTCTTCAAGAAGTCCTCGATTCTCTGCCCGCGCAAGgcggcg                          >  1:3927031/1‑150 (MQ=255)
                                                                                                                                ggCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTa                        >  2:7222729/1‑150 (MQ=255)
                                                                                                                                 gCGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTaa                       <  2:6387810/150‑1 (MQ=255)
                                                                                                                                        gggTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCa                <  2:1660270/150‑1 (MQ=255)
                                                                                                                                        gggTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCa                <  2:791942/150‑1 (MQ=255)
                                                                                                                                        gggTGGCGGCGCGGGACGCGCCGTAACACTTTCCAGAGGCCTGCATCGGTCAGTTGGGCGCGGTGACAGCCACCGAGTAGCGAAACTGCTGGTTACGGGACTTCATGACGACCTCGATGCTATGCCCGCGCATGGCGGCGTAAATGTCCa                >  2:5487928/1‑150 (MQ=255)
                                                                                                                                         ggTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAg               <  2:4601914/150‑1 (MQ=255)
                                                                                                                                         ggTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAg               <  2:4151403/150‑1 (MQ=255)
                                                                                                                                          gTGACGGCGCGTCCCGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCCCCGAGTAGCGACCCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCCTGGCGGCGTAAATGTCCAgg              >  2:6212497/1‑150 (MQ=255)
                                                                                                                                             gcggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTg           <  2:7322332/150‑1 (MQ=255)
                                                                                                                                             gcggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTg           <  1:487180/150‑1 (MQ=255)
                                                                                                                                             gcggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTg           >  1:1856120/1‑150 (MQ=255)
                                                                                                                                             gcggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTg           <  1:5667880/150‑1 (MQ=255)
                                                                                                                                             acggcgCGTCCCGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTg           >  2:2103670/2‑150 (MQ=255)
                                                                                                                                              gggcgGGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCaacctgg          <  1:7653431/149‑7 (MQ=255)
                                                                                                                                               ggcgGGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         <  1:406206/150‑1 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCGAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         >  1:8279179/1‑150 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         >  1:2201460/1‑150 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         <  1:1029580/150‑1 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         <  2:5175626/150‑1 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         >  1:5174350/1‑150 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         >  2:4728121/1‑150 (MQ=255)
                                                                                                                                               ggcgCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGa         >  2:4373593/1‑150 (MQ=255)
                                                                                                                                                     ggACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGACACAcc   <  2:6054415/150‑1 (MQ=255)
                                                                                                                                                      gACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAACCTGGacatttac  <  2:5224690/150‑6 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAA‑GC‑GCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCTGCCCGCGCATGGCGGCGTAAATGTCCAGGTTGGACACACCC  >  NC_002947/2870323‑2870620

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