Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,520 T→A 5.8% 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 5.8% 266.6 / 4.3 86Q149L (CAG→CTG) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base T (38/43);  new base A (3/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGC  >  NC_002947/5688372‑5688666
                                                                                                                                                     |                                                                                                                                                    
cGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                                      >  5:204432/1‑149 (MQ=255)
  aGCCAAGGATGTGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCATGGTCGCTGCCACCTgg                                                                                                                                                    <  4:27772/149‑1 (MQ=255)
  aGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCCACGGTCTCGGCCAGGGTCGCTGCCACCTg                                                                                                                                                     <  1:192473/149‑1 (MQ=255)
  aGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                                      >  4:208708/1‑147 (MQ=255)
  aGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                                      <  3:208708/147‑1 (MQ=255)
     cAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑g                                                                                                                                                 >  7:235689/1‑148 (MQ=255)
      aaGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑g                                                                                                                                                 >  5:240661/1‑147 (MQ=255)
      aaGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCa                                                                                                                                               <  3:160063/149‑1 (MQ=255)
        ggATGCGTTCGGCCAGCAACTGCTGGCGCTGGGCCACAGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGc                                                                                                                                             <  4:10611/149‑1 (MQ=255)
          aTGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGa                                                                                                                                           <  2:182102/149‑1 (MQ=255)
           tGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGa                                                                                                                                           >  8:88124/1‑148 (MQ=255)
             cGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACa                                                                                                                                         <  1:10916/148‑1 (MQ=255)
             cGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACa                                                                                                                                         <  6:128579/148‑1 (MQ=255)
                  ggCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACg                                                                                                                                    <  2:37433/148‑1 (MQ=255)
                   gCCAGCAACAGCTGGCGCTGGGCGACCGCAACCTGATTGGCCGGGGCACGGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                                  <  6:283147/149‑1 (MQ=255)
                   gCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                                  >  5:84706/1‑149 (MQ=255)
                     cAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCt                                                                                                                                <  7:88124/149‑1 (MQ=255)
                     cAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTG‑CACCAGGT‑ggcag                                                                                                                                             >  4:207820/1‑131 (MQ=255)
                     cAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTG‑CACCAGGT‑ggcag                                                                                                                                             <  3:207820/135‑5 (MQ=255)
                       gCAACAGCTGGCGCTGGTCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGc                                                                                                                              <  8:179588/149‑1 (MQ=255)
                          aCAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCg                                                                                                                           >  3:2596/1‑149 (MQ=255)
                           cAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCg                                                                                                                           >  8:59908/1‑148 (MQ=255)
                            aGCTGGCGCTGGGCCATCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCt                                                                                                                         <  8:265593/149‑1 (MQ=255)
                                 gcgcTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGc                                                                                                                      >  7:160629/1‑147 (MQ=255)
                                    cTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGga                                                                                                                 <  8:160629/149‑1 (MQ=255)
                                     tGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCACCGCCAGCACGGTCGGCTCGCGGGCCAGgat                                                                                                                >  6:111748/1‑149 (MQ=255)
                                       ggCCACCGCAACCTGAGTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTgg                                                                                                                                                    <  4:168970/112‑1 (MQ=255)
                                        gCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGgatgat                                                                                                             >  8:23312/1‑149 (MQ=255)
                                         ccACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGatgatg                                                                                                            >  2:311017/1‑149 (MQ=255)
                                          cACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGt                                                                                                           <  6:240661/149‑1 (MQ=255)
                                           aCCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                                  >  5:285344/1‑125 (MQ=255)
                                           aCCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGt                                                                                                                                  <  6:285344/125‑1 (MQ=255)
                                           aCCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTc                                                                                                          >  4:261154/1‑149 (MQ=255)
                                            ccGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCg                                                                                                         <  8:235689/149‑1 (MQ=255)
                                            ccGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCg                                                                                                         <  3:261154/149‑1 (MQ=255)
                                              gCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgt                                                                                                       <  1:311017/149‑1 (MQ=255)
                                                aaCCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                     >  5:211590/1‑149 (MQ=255)
                                                       ttGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTg‑cg                                                                                             >  7:27664/1‑149 (MQ=255)
                                                          gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGcc                                                                                                                     >  2:180208/1‑123 (MQ=255)
                                                          gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAg                                                                                          <  5:111748/149‑1 (MQ=255)
                                                          gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑CCAGCGACAGCACGGTCTGCTCGCTGGcc                                                                                                                     <  1:180208/123‑1 (MQ=255)
                                                              gggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCTCTGTCCAGGATGATGTAGGTTTTCTTG‑CGCAGGttt                                                                                      >  4:193735/1‑149 (MQ=255)
                                                               ggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGtttt                                                                                     >  5:73160/1‑149 (MQ=255)
                                                                gCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTc                                                                                    <  6:204432/149‑1 (MQ=255)
                                                                  aCCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTcc                                                                                   >  6:205217/1‑148 (MQ=255)
                                                                      ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTc                                                                              <  8:27664/149‑1 (MQ=255)
                                                                           cagcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCgcggc                                                                         <  6:211590/149‑1 (MQ=255)
                                                                              cagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTa                                                                      <  7:59908/149‑1 (MQ=255)
                                                                               agGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTAc                                                                     <  3:234298/149‑1 (MQ=255)
                                                                                 gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GTAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTAcc                                                                    >  7:237461/1‑148 (MQ=255)
                                                                                 gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTAcc                                                                    >  8:301643/1‑148 (MQ=255)
                                                                                      cGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGc                                                                                                                          >  2:315803/1‑90 (MQ=255)
                                                                                      cGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGc                                                                                                                          <  1:315803/90‑1 (MQ=255)
                                                                                       gACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                     >  7:163075/1‑110 (MQ=255)
                                                                                       gACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtgt                                                                                                     <  8:163075/110‑1 (MQ=255)
                                                                                       gACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                 <  1:191774/146‑1 (MQ=255)
                                                                                       gACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                 <  3:260324/145‑1 (MQ=255)
                                                                                        aCCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTcc                                                              <  2:161879/147‑1 (MQ=255)
                                                                                        aCCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTcc                                                              >  1:161879/1‑147 (MQ=255)
                                                                                          cACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtct                                                          >  2:26750/1‑149 (MQ=255)
                                                                                             tttttCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGt                                                       >  1:90543/1‑149 (MQ=255)
                                                                                             tttttCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGt                                                       >  5:18025/1‑149 (MQ=255)
                                                                                               tttCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTc                                                      <  7:23312/148‑1 (MQ=255)
                                                                                                 tCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCAc                                                    <  5:205217/148‑1 (MQ=255)
                                                                                                        ccACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCg                                              >  1:81684/1‑147 (MQ=255)
                                                                                                         cACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTc                                                      >  1:87102/1‑138 (MQ=255)
                                                                                                         cACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTc                                                      <  2:87102/138‑1 (MQ=255)
                                                                                                            cTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                                        <  2:81684/149‑1 (MQ=255)
                                                                                                              gCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                                        >  1:114220/1‑147 (MQ=255)
                                                                                                                aGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGgcg                                     >  3:102209/1‑148 (MQ=255)
                                                                                                                  cTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAg                                  >  6:135767/1‑149 (MQ=255)
                                                                                                                   tGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGCGGCGCAgg                                 >  1:225514/1‑149 (MQ=255)
                                                                                                                     cGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggc                                <  2:114220/148‑1 (MQ=255)
                                                                                                                     cGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAgg                                 <  6:84706/147‑1 (MQ=255)
                                                                                                                     cGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGCGGCGCAggc                                <  2:225514/148‑1 (MQ=255)
                                                                                                                      gg‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggcgg                              <  2:90543/149‑1 (MQ=255)
                                                                                                                       g‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCTCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggcgg                              >  8:285066/1‑148 (MQ=255)
                                                                                                                          aCGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGATCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAg                           >  1:67946/1‑149 (MQ=255)
                                                                                                                          aCGGTCTCGGCCAGGGTCGCTGCCACCTGGG‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAg                           <  5:287452/149‑1 (MQ=255)
                                                                                                                              tctcGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGc                         >  3:245801/1‑147 (MQ=255)
                                                                                                                              tctcGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCa                       <  4:2596/149‑1 (MQ=255)
                                                                                                                                 cGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTg                    >  5:278959/1‑149 (MQ=255)
                                                                                                                                   cgtgctgTCGCTG‑CACCAGGTGGCAGCGGCCGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGa                  >  1:179207/7‑149 (MQ=255)
                                                                                                                                      aGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt               <  3:193735/149‑1 (MQ=255)
                                                                                                                                       ctgTCGCTG‑CACCAGGTGGCAGCGACCGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt               <  2:179207/146‑1 (MQ=255)
                                                                                                                                         gTCGCTG‑CACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGt             >  8:153765/1‑148 (MQ=255)
                                                                                                                                              tGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCt       >  5:51768/1‑149 (MQ=255)
                                                                                                                                              tGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCt       >  8:168404/1‑149 (MQ=255)
                                                                                                                                                   ccTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgcgc  <  5:135767/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGG‑CACGGTCTCGGCCAGGGTCGCTGCCACCTGGT‑GCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTG‑CGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGC  >  NC_002947/5688372‑5688666

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