Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,520 T→A 6.5% 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 6.5% 280.1 / 6.2 91Q149L (CAG→CTG) pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base A (4/2);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688372‑5688665
                                                                                                                                                     |                                                                                                                                                  
cGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCt                                                                                                                                                    >  8:90798/1‑149 (MQ=255)
       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                                                                                                                                     >  6:12705/1‑146 (MQ=255)
                 cGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGc                                                                                                                                     <  5:12705/146‑1 (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                                                                                                                             <  8:29936/143‑1 (MQ=255)
                            aGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                             >  7:29936/1‑143 (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                                                                                                                             >  6:320393/1‑101 (MQ=255)
                                                                      ctctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTg                                                                                                                             <  5:320393/101‑1 (MQ=255)
                                                                       tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                        <  3:231899/137‑1 (MQ=255)
                                                                       tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAgg                                                                                        >  4:231899/1‑137 (MQ=255)
                                                                       tctGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcg                                                                            >  8:68099/1‑149 (MQ=255)
                                                                            agcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCgcggcg                                                                       >  3:298796/1‑149 (MQ=255)
                                                                            agcagGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCgcggcg                                                                       >  1:171609/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:347733/149‑1 (MQ=255)
                                                                                         ccACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCAtc                                                          >  7:155433/1‑149 (MQ=255)
                                                                                          cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                           <  7:167978/131‑1 (MQ=255)
                                                                                          cACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTcgc                                                                           >  8:167978/1‑131 (MQ=255)
                                                                                             tttttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGt                                                      >  1:374136/1‑149 (MQ=255)
                                                                                                ttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                    <  6:42610/100‑1 (MQ=255)
                                                                                                ttCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGgtg                                                                                                    >  5:42610/1‑100 (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                                       >  8:83673/1‑130 (MQ=255)
                                                                                                                               tcGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTg                                       <  7:83673/130‑1 (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              <  8:155433/149‑1 (MQ=255)
                                                                                                                                      aGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACtt              <  2:128210/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)
                                                                                                                                                     |                                                                                                                                                  
CGAGCCAAGGATGCGTTCGGCCAGCAACAGCTGGCGCTGGGCCACCGCAACCTGATTGGCCGGGGCACCGCTCTGCAGCAGGATCTCGACCACTTTTTCGTATTCCACCTGCAGCTGCGGCACGGTCTCGGCCAGGGTCGCTGCCACCTGGTGCAGCGACAGCACGGTCTGCTCGCTGGCCAGGATGATGTCGGTGTTCTTGCGCAGGTTTTCCCAGTCGCGGCGTACCGATTCCATCTCGTCACGCACCGCAGGTGGCGCAGGCGGCAGGCCAGTGGACTTGTCACCCTGGCG  >  NC_002947/5688372‑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: