Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,523:1 +G 8.3% coding (443/2061 nt) pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,5231.G 8.3% 227.5 / 13.2 60coding (443/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base . (28/27);  new base G (3/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAGG  >  NC_002947/5688376‑5688669
                                                                                                                                                      |                                                                                                                                                   
ccAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑g                                                                                                                                                    <  8:99868/149‑1 (MQ=255)
      attcGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACATGAT‑‑TGGCCGGGGCACCGCTATGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCAAGGGTCGCTGCCACCTGGT‑GCA‑gag                                                                                                                                              >  4:201821/4‑146 (MQ=255)
      aTGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGa                                                                                                                                             <  5:161414/149‑1 (MQ=255)
         cgctCGGCCAGACACAGCGGGCGCTGGGACACCTAAACCTGATCCTCAAAGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTG‑CACCAGGTGGCA‑GCGa                                                                                                                                             <  1:308373/145‑1 (MQ=255)
         cGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACa                                                                                                                                           <  6:150572/148‑1 (MQ=255)
          gTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAg                                                                                                                                          >  7:86296/1‑148 (MQ=255)
           ttCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCa                                                                                                                                        >  2:121254/1‑149 (MQ=255)
             cGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACg                                                                                                                                      <  4:306746/149‑1 (MQ=255)
              ggCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACgg                                                                                                                                     >  4:175982/1‑149 (MQ=255)
                 cAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACTGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCt                                                                                                                                  <  3:201821/149‑1 (MQ=255)
                   gCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGc                                                                                                                                <  4:166679/149‑1 (MQ=255)
                      aCAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCACTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTG‑CACCAGGTGGCA‑GCGa                                                                                                                                             <  6:202622/133‑1 (MQ=255)
                      aCAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCACTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTG‑CACCAGGTGGCA‑GCGa                                                                                                                                             >  5:202622/1‑133 (MQ=255)
                          cTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTgg                                                                                                                         >  5:301224/1‑149 (MQ=255)
                              cgcTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAg                                                                                                                     <  6:301224/149‑1 (MQ=255)
                                 tGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGgat                                                                                                                  >  6:114375/1‑149 (MQ=255)
                                         cGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCgg                                                                                                          >  4:342412/1‑149 (MQ=255)
                                              ccTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTc                                                                                                     >  1:185045/1‑149 (MQ=255)
                                               cTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCt                                                                                                    >  5:130798/1‑149 (MQ=255)
                                                 gAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTc                                                                                                            <  3:108853/139‑1 (MQ=255)
                                                        gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgt                                                                                                         <  2:132708/137‑1 (MQ=255)
                                                         ccGGGGCACCGCCCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                            >  2:17363/1‑149 (MQ=255)
                                                            gggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                        >  3:256598/1‑134 (MQ=255)
                                                            gggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                        <  4:256598/134‑1 (MQ=255)
                                                                  cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCg                                                                                                                             >  7:23663/1‑107 (MQ=255)
                                                                  cGCTCTGCAGAAGGATCTCGACCACTTTTTCGTATTCTACCTTCTGTTGAGGCACGGTCTCGGCGAGTGTCGCTTCCTCCTGGT‑GCA‑GCGAGAGCACTTTCTGCTCg                                                                                                                             <  8:23663/107‑1 (MQ=255)
                                                                    ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTc                                                                                 >  1:125437/1‑149 (MQ=255)
                                                                       tGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgcg                                                                              >  3:248465/1‑149 (MQ=255)
                                                                         cagcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCgcgg                                                                             <  3:175982/148‑1 (MQ=255)
                                                                            cagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTa                                                                         <  1:17363/149‑1 (MQ=255)
                                                                             agGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTAc                                                                        <  1:121254/149‑1 (MQ=255)
                                                                                 gctcGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                    <  4:11815/148‑1 (MQ=255)
                                                                                     gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGa                                                                     <  4:248465/144‑1 (MQ=255)
                                                                                      aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                               <  8:86296/149‑1 (MQ=255)
                                                                                       ccACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtc                                                              <  5:114375/149‑1 (MQ=255)
                                                                                       ccACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtc                                                              >  5:189443/1‑149 (MQ=255)
                                                                                          cTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑Gc                                                                                                                                               >  5:2917/1‑65 (MQ=255)
                                                                                          cTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑Gc                                                                                                                                               <  6:2917/65‑1 (MQ=255)
                                                                                           tttttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGga                                                                                                                   >  4:169747/1‑92 (MQ=255)
                                                                                           tttttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGga                                                                                                                   <  3:169747/92‑1 (MQ=255)
                                                                                             tttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGt                                                          <  6:130798/147‑1 (MQ=255)
                                                                                               tCGTATTCCACCTGCAGCTGCGGCAAGGTCTCGGCCAGGGTAGCTTCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTAGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTGCATAACGTCCCg                                                      >  7:234346/1‑149 (MQ=255)
                                                                                                  tATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCAc                                                   >  5:132771/1‑149 (MQ=255)
                                                                                                   aTTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCAcc                                                  >  3:317240/1‑149 (MQ=255)
                                                                                                    ttCCACCTGCAGCTGCGGCACGGTCTCGGCGAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                       >  1:30746/1‑95 (MQ=255)
                                                                                                    ttCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                       <  2:30746/95‑1 (MQ=255)
                                                                                                    ttCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCg                                                 <  3:342412/149‑1 (MQ=255)
                                                                                                        aCCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAgg                                             >  4:35031/1‑149 (MQ=255)
                                                                                                         ccTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                            <  3:35031/149‑1 (MQ=255)
                                                                                                            gCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                               <  1:295557/111‑1 (MQ=255)
                                                                                                            gCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                               >  2:295557/1‑111 (MQ=255)
                                                                                                             cAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTgg                                          <  6:189443/147‑1 (MQ=255)
                                                                                                              aGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGgcgc                                       <  2:185045/149‑1 (MQ=255)
                                                                                                                     gCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggcggc                                >  8:239951/1‑149 (MQ=255)
                                                                                                                           tctcGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCAGGCCAGGATGATGACGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGc                            >  2:82594/1‑147 (MQ=255)
                                                                                                                                   aGGGTCGGGGCCAACTGGG‑GCATGCTA‑AGCACGTGATGCTCGCTTGCCAGGTTGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTACATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt                  <  8:234346/149‑1 (MQ=255)
                                                                                                                                     tgTCGCTG‑CACCAGGTGGCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGt                >  1:238380/2‑149 (MQ=255)
                                                                                                                                            tgcACCAGGTGGCA‑GCGACAGCACGGTCTGCTCCCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCt          >  7:58482/3‑149 (MQ=255)
                                                                                                                                              cACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgc       >  6:156123/1‑149 (MQ=255)
                                                                                                                                                 cTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgcg      >  7:324022/1‑147 (MQ=255)
                                                                                                                                                   ggT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAgg  >  4:118107/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAGG  >  NC_002947/5688376‑5688669

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