Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,523:1 +G 6.5% 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 6.5% 364.1 / 19.9 92coding (443/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base . (45/41);  new base G (4/2);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

AGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAGG  >  NC_002947/5688379‑5688669
                                                                                                                                                  |                                                                                                                                                  
aGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAg                                                                                                                                                <  1:306110/149‑1 (MQ=255)
 ggATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGc                                                                                                                                               <  2:47962/149‑1 (MQ=255)
      cGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAg                                                                                                                                          <  4:362371/149‑1 (MQ=255)
         tCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCAc                                                                                                                                       >  7:349750/1‑149 (MQ=255)
          cGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGc                                                                                                                                         <  5:12705/146‑1 (MQ=255)
          cGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGc                                                                                                                                         >  6:12705/1‑146 (MQ=255)
           ggCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACg                                                                                                                                      <  2:235323/148‑1 (MQ=255)
            gCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGt                                                                                                                                    >  4:20205/1‑149 (MQ=255)
            gCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGt                                                                                                                                    >  5:148494/1‑149 (MQ=255)
             ccAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTc                                                                                                                                   >  6:236724/1‑149 (MQ=255)
                     aGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 >  7:29936/1‑143 (MQ=255)
                     aGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 <  8:29936/143‑1 (MQ=255)
                     aGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCt                                                                                                                           <  8:349750/149‑1 (MQ=255)
                      gCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGCCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTg                                                                                                                          >  8:265741/1‑149 (MQ=255)
                                  ccACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 <  1:322254/130‑1 (MQ=255)
                                  ccACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 >  2:322254/1‑130 (MQ=255)
                                      cGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCgg                                                                                                          >  5:105606/1‑149 (MQ=255)
                                               aTTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTTGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCtt                                                                                                   <  3:20205/147‑1 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTc                                                                                                                              <  5:75154/116‑1 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTc                                                                                                                              >  6:75154/1‑116 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAg                                                                                             >  5:289423/1‑149 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAg                                                                                             <  3:265030/149‑1 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCACCAGGTGGCAGCGAc                                                                                                                                            >  4:250396/1‑102 (MQ=255)
                                                   gCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCACCAGGTGGCAGCGAc                                                                                                                                            <  3:250396/102‑1 (MQ=255)
                                                    ccGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                            >  7:168014/1‑149 (MQ=255)
                                                     cGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGt                                                                                           >  6:61070/1‑149 (MQ=255)
                                                        ggCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGtttt                                                                                        <  5:236724/149‑1 (MQ=255)
                                                         gCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTc                                                                                       <  6:148494/149‑1 (MQ=255)
                                                           aCCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGCTGTTCTTGCGCAGGTATTcc                                                                                      >  1:112832/1‑149 (MQ=255)
                                                             cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCCCCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAg                                                                                   >  4:37974/1‑149 (MQ=255)
                                                             cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGGACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTg                                                                                                  >  6:150811/1‑134 (MQ=255)
                                                             cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTg                                                                                                  <  5:150811/134‑1 (MQ=255)
                                                             cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAg                                                                                   <  6:289423/149‑1 (MQ=255)
                                                             cGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTACACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAg                                                                                   >  1:101105/1‑149 (MQ=255)
                                                               ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 <  5:320393/101‑1 (MQ=255)
                                                               ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                                 >  6:320393/1‑101 (MQ=255)
                                                                tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                            <  3:231899/137‑1 (MQ=255)
                                                                tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                            >  4:231899/1‑137 (MQ=255)
                                                                tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcg                                                                                >  8:68099/1‑149 (MQ=255)
                                                                     agcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCgcggcg                                                                           >  1:171609/1‑149 (MQ=255)
                                                                     agcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCgcggcg                                                                           >  3:298796/1‑149 (MQ=255)
                                                                          gATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTAcc                                                                       >  8:347733/1‑148 (MQ=255)
                                                                           aTCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGa                                                                     >  5:302260/1‑149 (MQ=255)
                                                                            tctcGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTCTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                    >  1:369584/1‑149 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGgatga                                                                                                                <  8:350240/101‑1 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGgatga                                                                                                                >  7:350240/1‑101 (MQ=255)
                                                                                gACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGAt                                                                    <  7:90798/145‑1 (MQ=255)
                                                                                 acgtcTTTTTCGTATTCCACCTGCAGCTGCGGCTCGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                               <  2:112832/145‑1 (MQ=255)
                                                                                 aCCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAt                                                               <  5:61070/149‑1 (MQ=255)
                                                                                  ccACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtc                                                              >  7:155433/1‑149 (MQ=255)
                                                                                  ccACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtc                                                              <  7:347733/149‑1 (MQ=255)
                                                                                   cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                               <  7:167978/131‑1 (MQ=255)
                                                                                   cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                               >  8:167978/1‑131 (MQ=255)
                                                                                      tttttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGt                                                          >  1:374136/1‑149 (MQ=255)
                                                                                         ttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                        >  5:42610/1‑100 (MQ=255)
                                                                                         ttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                        <  6:42610/100‑1 (MQ=255)
                                                                                         ttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCa                                                        <  3:319287/148‑1 (MQ=255)
                                                                                            gTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGc                                                     <  2:152242/148‑1 (MQ=255)
                                                                                             tATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCAc                                                   >  7:256267/1‑149 (MQ=255)
                                                                                                  cACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAg                                              >  1:142983/1‑149 (MQ=255)
                                                                                                  cACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAg                                              <  2:171609/149‑1 (MQ=255)
                                                                                                    ccTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCa                                                                                    >  5:111628/1‑109 (MQ=255)
                                                                                                    ccTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCa                                                                                    <  6:111628/109‑1 (MQ=255)
                                                                                                    ccTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGt                                            <  2:101105/149‑1 (MQ=255)
                                                                                                          gCTGCGGCACGGTCTCGGCCAGGGCCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCa                                      <  7:265741/149‑1 (MQ=255)
                                                                                                                gCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggcgg                                 >  8:335669/1‑148 (MQ=255)
                                                                                                                gCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggcgg                                 >  7:276768/1‑148 (MQ=255)
                                                                                                                 cACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTTTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCGATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCa                               <  2:369584/149‑1 (MQ=255)
                                                                                                                    ggTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGc                            <  6:105606/149‑1 (MQ=255)
                                                                                                                      tctcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGc                            >  1:340918/1‑147 (MQ=255)
                                                                                                                       ctcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAg                         >  3:150148/1‑149 (MQ=255)
                                                                                                                        tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                                           <  7:83673/130‑1 (MQ=255)
                                                                                                                        tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                                           >  8:83673/1‑130 (MQ=255)
                                                                                                                        tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGt                        >  1:128210/1‑149 (MQ=255)
                                                                                                                            gCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGa                     <  3:249263/149‑1 (MQ=255)
                                                                                                                             ccAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCa                               >  2:85154/1‑138 (MQ=255)
                                                                                                                             ccAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCa                               <  1:85154/138‑1 (MQ=255)
                                                                                                                               aGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt                  <  2:374136/149‑1 (MQ=255)
                                                                                                                               aGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt                  <  2:128210/149‑1 (MQ=255)
                                                                                                                               aGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt                  <  8:155433/149‑1 (MQ=255)
                                                                                                                               aGGGTCGCTGCCACCTGGGGCAGCGACAGCACCGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt                  <  8:168014/149‑1 (MQ=255)
                                                                                                                                 tgTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGt                >  8:313987/2‑149 (MQ=255)
                                                                                                                                  gTCGCTGCACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGt                >  2:260616/1‑148 (MQ=255)
                                                                                                                                   tCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCa              >  2:353165/1‑149 (MQ=255)
                                                                                                                                     gCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCAcc            >  2:268141/1‑149 (MQ=255)
                                                                                                                                          cACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgc       >  8:32737/1‑149 (MQ=255)
                                                                                                                                          cACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTAGTCACCCTGgc       >  3:159876/1‑149 (MQ=255)
                                                                                                                                          cACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggc                                   <  5:121434/122‑1 (MQ=255)
                                                                                                                                          cACCAGGTGGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAggc                                   >  6:121434/1‑122 (MQ=255)
                                                                                                                                           aCCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCAGgcg      <  1:353165/149‑1 (MQ=255)
                                                                                                                                             cTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGgcg      >  6:311831/1‑147 (MQ=255)
                                                                                                                                                gTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAgg  <  7:335669/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCGCAGG  >  NC_002947/5688379‑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: