Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 856,078 G→A 8.0% A9A (GCG→GCA pagL‑I → lipid A 3‑O‑deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947856,0780GA 8.0% 173.1 / 4.4 62A9A (GCG→GCApagL‑Ilipid A 3‑O‑deacylase
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAGCGTTGG  >  NC_002947/855941‑856224
                                                                                                                                                |                                                                                                                                                      
gTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCAGCCGCTG‑‑‑‑ccaa                                                                                                                                         <  8:114400/149‑3 (MQ=255)
gTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCAGCCGCTG‑‑‑‑cc                                                                                                                                           >  4:23946/1‑147 (MQ=255)
   gtgAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCAGCCGCTG‑‑‑‑ccaagc                                                                                                                                       <  3:23946/148‑5 (MQ=255)
    tgAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGtt                                                                                                                                     >  5:244879/1‑149 (MQ=255)
    tgAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGt                                                                                                                                      <  6:83066/148‑1 (MQ=255)
     gAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGttt                                                                                                                                    >  5:118170/1‑149 (MQ=255)
     gAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGttt                                                                                                                                    >  6:210445/1‑149 (MQ=255)
      aaaaaaaCGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGACTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGtttt                                                                                                                                   >  1:198773/1‑149 (MQ=255)
           aaCGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGc                                                                                                                              <  6:296260/149‑1 (MQ=255)
            aCGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGc                                                                                                                              >  5:113346/1‑148 (MQ=255)
              gTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGCGAGGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGc                                                                                                                              >  1:242317/1‑148 (MQ=255)
               tGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCAGCCGCTG‑‑‑‑ccaagc                                                                                                                                       >  1:180868/1‑132 (MQ=255)
               tGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCAGCCGCTG‑‑‑‑ccaagc                                                                                                                                       <  2:180868/136‑5 (MQ=255)
                  aaaaaGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGAt                                                                                                                       <  3:239984/149‑1 (MQ=255)
                    agaagtggAACGATAGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGt                                                                                                                     >  7:216040/8‑149 (MQ=255)
                           gAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGca                                                                                                                >  6:242535/1‑147 (MQ=255)
                                 tCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGc                                                                                                        <  6:244879/149‑1 (MQ=255)
                                      ccACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTgg                                                                                                                                 <  3:265265/119‑1 (MQ=255)
                                      ccACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTgg                                                                                                                                 >  4:265265/1‑119 (MQ=255)
                                               cTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGtt                                                                                            <  6:113346/147‑1 (MQ=255)
                                                tttCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTc                                                                                         >  6:41420/1‑149 (MQ=255)
                                                 ttCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGc                                                                                                                 <  3:109918/124‑1 (MQ=255)
                                                 ttCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGc                                                                                                                 >  4:109918/1‑124 (MQ=255)
                                                  tCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCg                                                                                        <  8:200752/148‑1 (MQ=255)
                                                      ccGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTgg                                                                                    <  3:46529/148‑1 (MQ=255)
                                                       cGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGGAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCg                                                                                        <  8:85486/143‑1 (MQ=255)
                                                       cGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGGAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCg                                                                                        >  7:85486/1‑143 (MQ=255)
                                                       cGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTgggg                                                                                  <  5:210445/149‑1 (MQ=255)
                                                         cTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTgggg                                                                                  >  1:286052/1‑147 (MQ=255)
                                                             gCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACt                                                                            >  7:263343/1‑149 (MQ=255)
                                                              ccTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGCAGCAAAGCGGCTGCGGAAAAGATTTTTTTTGGGGAAATCGATGTCGGCACAGGCTGCCGATGTTTGGTTTTCGGTgg                                                                                    >  1:95145/1‑149 (MQ=255)
                                                              ccTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTg                                                                           >  4:240397/1‑149 (MQ=255)
                                                                tGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGc                                                                         >  7:89424/1‑149 (MQ=255)
                                                                  agaTGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGa                                                                       >  7:111485/1‑149 (MQ=255)
                                                                     tGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTc                                                                    <  8:9100/149‑1 (MQ=255)
                                                                       aGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTACGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGa                                                                  >  4:205211/1‑149 (MQ=255)
                                                                             ggACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGa                                                                  >  6:77043/1‑145 (MQ=255)
                                                                             ggACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGa                                                                  <  5:77043/145‑1 (MQ=255)
                                                                               aCAAATAACAGCATCAGAAAGAAGGATGTTTCTTACGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATgg                                                            <  3:205211/149‑1 (MQ=255)
                                                                                cAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGt                                                           <  4:264360/149‑1 (MQ=255)
                                                                                      aCAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTAcc                                                      <  8:111485/148‑1 (MQ=255)
                                                                                         gCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTgg                                                                          <  5:227916/125‑1 (MQ=255)
                                                                                         gCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTgg                                                                          >  6:227916/1‑125 (MQ=255)
                                                                                                aaagaaGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTg                                           >  5:269502/1‑149 (MQ=255)
                                                                                                    aagGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAAt                                       >  2:8361/1‑149 (MQ=255)
                                                                                                     agGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATc                                      <  6:118170/149‑1 (MQ=255)
                                                                                                        aTGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGaa                                   <  8:263343/149‑1 (MQ=255)
                                                                                                          gTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGtttt                                                                                                                                   <  4:93294/51‑1 (MQ=255)
                                                                                                          gTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGtttt                                                                                                                                   >  3:93294/1‑51 (MQ=255)
                                                                                                           tttCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTg                                <  2:198773/149‑1 (MQ=255)
                                                                                                           tttCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTg                                <  8:216040/149‑1 (MQ=255)
                                                                                                                 attaGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGAc                            <  5:41420/144‑1 (MQ=255)
                                                                                                                  tGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGAcgc                          <  8:89424/148‑1 (MQ=255)
                                                                                                                    aGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGa                        <  2:242317/148‑1 (MQ=255)
                                                                                                                       aaaCTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTtg                    <  2:286052/149‑1 (MQ=255)
                                                                                                                       aaaCTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTtg                    <  3:240397/149‑1 (MQ=255)
                                                                                                                              tctgCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGATTTCGGTGCGGAAGACTGGCGATACGACAATGGTCTACCGGATAGGGCTGCAATCGAATTGGGAAGCGAGATCGTCGCa              >  8:251482/4‑148 (MQ=255)
                                                                                                                              tCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTtggtgg                >  2:95734/1‑146 (MQ=255)
                                                                                                                                   ttGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACAGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAg        >  5:196965/1‑149 (MQ=255)
                                                                                                                                   ttGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTAGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTACGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAg        >  4:68721/1‑149 (MQ=255)
                                                                                                                                             aGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAGCGTTg   <  5:242535/149‑1 (MQ=255)
                                                                                                                                              gCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAGCGTTgg  >  1:265808/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GTTGTGAAAAAAACGTGAAAAAAGGCTGAACTATCGTTCCACTTCTGCTTTCTACCGCTTTGCCTGAGATGAGGC‑‑GGACAAATAACAGCATCAGAAAGAAGGATGTTTCTTATGAGGAAACTGCTCGGCTTGGC‑‑‑‑‑AGCGGCTGCTG‑‑‑‑CCTTTGTTTTGGGGCAAGCGATGTCGGCACAGGCTGCCGATGTTTCGTTTTCGGTGGGGCAGACTGGCGATTCGACAATGGTCTACCGGCTAGGGCTGCAATCGAATTGGGACGCGAGTTGGTGGCAGACCAGCGTTGG  >  NC_002947/855941‑856224

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