Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,520 T→A 8.3% Q149L (CAG→CTG)  pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,5200TA 8.3% 174.2 / 4.9 60Q149L (CAG→CTG) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base A (3/2);  total (30/30)
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‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688376‑5688665
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                                         >  5:202622/1‑133 (MQ=255)
                      aCAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCACTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTG‑CACCAGGTGGCA‑GCGa                                                                                                                                         <  6:202622/133‑1 (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                                                                                                    <  4:256598/134‑1 (MQ=255)
                                                            gggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                    >  3:256598/1‑134 (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                                                                                                               <  3:169747/92‑1 (MQ=255)
                                                                                           tttttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGga                                                                                                               >  4:169747/1‑92 (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)
                                                                                                                                                  |                                                                                                                                                   
CCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGAT‑‑TGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCA‑GCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688376‑5688665

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: ‑